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PREFACE 


NATA is a ccxnputer program for calculating steady, quasi- 
one-dimensional flow of a reacting gas mixture in ^ nczzle or 
rectangular channel. It also computes stagnation-point condi- 
tions on axisynimetric or two-dimensional models and the cond.i- 
tions on the flat surface of a blunt wedge. The code's primary 
purpose is the prediction and interpretation of test conditions 
in arc-heated wind tunnels used for laboratory evaluation of 
thermal protection materials for reentry vehicles such as the 
Space Shuttle Orbiter. The theory and analysis underlying the 
operation of NATA have been documented in Volume I of this re- 
port.* The present volume is a user's manual for the code. 
defines the inputs and outputs, documents the precoded data or 
gas species and nozzle geometries, explains the diagnostic out- 
puts, and includes illustrative results from test problems. In 
addition, this volume contains a user's manual for an auxiliary 
progreim (NOZFIT) which can be used to set up nozzle profile 
curvefits of the form used in NATA. The programming of NATA and 
NOZFIT is documented in Volume III.** 


*V7, L. Bade and J. M. Yos, The NATA Code - Theory and Analysis, 
NASA CR-2547. 


**VJ. L. Bade and J. M. Yos, The NATA Code - Programmer's Manvial, 
N7\SA CR- 14 174 4. 
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TI-IE NATA CODE - JMANUAL 

By W. L, Bade and J. M, Yos 
Avco Systems Division 
VJilmington, Massachusetts 


1. INTRODUCTION 

The N/vTA code is a computer program for solving proljlcms 
of steady, quasi-onc>-dimcns.i onal gas flov; in nossles. The code's 
capabilities, and the theory and analysis underlying its opera- 
tion, have been documented in the first volume of this report 
(ref. 1). The present volume is a user's m^inual for the code. 
Section 2 is a comprehensive discussion of NATA' s inputs. Sec- 
tion 3 defines the normal outputs which present the results of 
the flow calculations and the calculations of test conditions on 
models. Section 4 documents the preceded data on gas models, 
transport properties, and nozzle geometries. 7'ppendix A discusses 
the reaction data assiamcd in the electronic nonequ.ilibrium models 
for argon and helixim. Appendix B lists and cznplains the diagnos- 
tic outputs v.’hich NATA produces to aid the user in identifying 
the. causes of abnormal conditions and code f;iilu?:cs. Apjjer.dix C 
presents the ii^puts and outputs of test problems which illustrate 
the code's use and capabilities. Finally Apj.)andix D is a user's 
manual for the NOZFIT code, relatively small au.'ziliary program 
for setting up nozzle profile curvefits of the tom required by 
NATA. 
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2, INPUTS 

NATA employs a flexible, user— oriented input system 
called "Namelist", which is a standard feature of the Fortram 
IV programming language. The format requirements for the in- 
put card deck are summarized in Section 2,1. Section 2.2 dis- 
cusses the few inputs that are required for running most prob- 
lems of interest to NASA/JSC. Section 2.3 is a complete list 
of the definitions of all input variables accepted by NATA, 
except those used to read in the properties of chem'.cal species 
and the rates of reactions. Input of these types of gas data 
is discussed in Section 2.4. In this discussion, the reader is 
assumed to be generally familiar with NATA' s capabilities as 
summarized in the Introduction of Volume I (ref. 1) . 

Excunples of sets of NATA input data for various types of 
problem, together with portions of the output produced when the 
code was run with these data, are presented in Appendix C. 

2.1 Format Requirements for Input Cards 

The input data for a NATA code run are punched on com- 
puter cards. The data for each case in tlie run require a deck 
of at least four cards, as explained below. Tlie decks for the 
cases are stacked to obtain the input deck for the entire run 
or job. The cases are run in tl^.e order in v;hich they appear 
in the job deck, from the top dov;n. Figure 1 illustrates the 
input data for a job consisting of three cases. In this fig- 
ure, the data are written* on an 80-column coding form. Each 
line of the form corresponds to a card in tlie input deck. 

Each column corresponds to one of the 80 columns in which 
data can be punched on a computer card. 


*In figure 1, and in the text of this report, the letter 0 
is written 0 to distinguish it from zero (0), and the letter 
I is written v/ith serifs to distinguish it from the numeral 1. 
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The first card in the input deck for each case may con- 
tain any descriptive information desired by the user. This 
information is reproduced at the head of the printed output 
for the case. This card may be left blank, if desired; but 
it may not be left out. 

The second card in the deck for each case must contain 
the follcv/ing characters: 


$ INPUT 


in columns 2-7.* 

The nxime'^ical input data for the case begin on the third 
card of the deck for the case. As many additional cards may 
be used as necessary. The data are ptinched in the form of 
equations : 


or 


variable name = value (la) 

array name = list of values (Ib) 


These equations are separated by commas. Also, the individ- 
ual values of a list being read into an array are separated 
by commas. The admissible names for input variables and 
arrays are listed and defined in Sections 2.2 and 2.3, be- 
low. 


The follov/ing is a condensed summary of the format re- 
quirements for the data cards. A more ccroplete discussion 
may be found in the UNIVAC and IBM Fortran manuals (refs. 2, 
3), under "Nciraelist". 


*V\’lien the program is run on an IBM 360 system, the $ is re- 
placed by an EBCDIC ampersand or a BCD + sign. 
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(1) T'ne data must be punched in columns 2-80 of 
the cards; column 1 should be blank. 

(2) Tlierc may be no embedded bl^lnks v/ithin the 
field occupied by a variable name or a num- 
erical value (including the sign, if ciny) . 

With this exception, blanks may be inserted 
freely to improve legibility cind thus facili- 
tate checking the ca^is. 

(3) The last item on each data card must be a 
constant (i.e., a value) followed by a comma. 

On the last data card for a case, the final 
comma is optional. 

(4) The nximber of values listed for an array must 
be less than or equal to the "number of entries" 
or the product of dimensions given in Sections 
2.2 to 2.4. 

(5) The value given for any variable (except a log- 
ical variable) may be a mamber with or without 
a decimal point, or a number with a decimal 
follov/ed by an exponent of 10 expressed in "E" 
notation. For example, 1.23 x 10”^ could be 
punched as 1.23E— 5, and 9.8 x IC^^ as 9.8E+23 
or 9.8E23. For example, see the input for 
TPRINTI in figure 1. 

(6) Tlie typing of variables as integer or real 
(see Section 2.2) is determined by the pro- 
gram. If a nxunber without decimal is pro- 
vided as input to a real variable, the name- 
list input system converts it to a real value 
before storing it. For example, in figure 1, 
the specification ANGLE = 10 has the same 
effect as 7\NGLE = 10.0. If a decimal number 
is provided as input to em integer variable, 
the s^'stem rounds it down to the next smailest 
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integer. For example, NANGLE = 1.3 would be 
« juivalent to NANGLE = 1; however, NANGLE =5.0 
ight have the effect of NANGLE = 4 on some com- 
puter systems where the 5.0 is represented as 
4.999999... Thus, the decimal should be omitted 
in input values to integer variables. 

Note: Tlie remaining features (7-10) in the present list are 

not needed for setting up standard-type cases using the 
inp its of Section 2.2. 

(7) The input value to a variable typed "logical" (see 
Section 2.3) may be T for "true" or F for "false". 
For example, see the input for AAMS in figure 1. 

(8) xn an input of array data in the form (lb), if sev- 
eral successive values in the list arc equal to the 
same value v, these values can be given in the form 
n'^v, where n is the number of values equal to v. 

(9’ In the case of a multiply dimensioned array, the 

order in which the values must be listed is deter- 
:nined by the rule that the left-most index varies 
no3t rapi-^ry and the right-most index least rap- 
idly. I or example, in a doubly dimensioned array 
such as ISHAPE(J, M) , which is dimensioned (12,2), 
the values must be listed in the order 

ISHAPEd, 1) 

ISHAPE(2,1) 

ISILAPE(3,1) 


• 

ISIIAPE(12, 1) 

ISHAPE(1,2) 

ISI1APE(2,2) 


ISILAPE(12,2) 
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It is not necessary to set all of the elements 
of an array in the injput list. However, the 
list in equation (lb) must begin with the first 
element and must include values for all elements 
up to the last one to which a value is assigned. 

For example, in a channel problem v/ith four sec- 
tions in each profile, the input for ISHAPE might 
be 

ISHAPE = 1, 2,2,1, 8*0, 1,2, 2,1, 

The entry 8*0 fills up the elements ISHAPE(5,1), 
(ISHAPE(6, 1) , . . . , ISHAPE(12, 1) , v;hich are not 
actually to be used. This entry is required to 
place the remaining data (1,2, 2,1) into the loca- 
tions ISHAPE (1,2;, ISHAPE (4, 2). If the 8*0 

entry were omitted, the second set of data (1, 2,2,1) 
v;ould be loaded into ISHAPE (5,1),..., ISH/'I>E (8,1), 
instead. Note that the elements ISHAPE(5,2) , . . . , 
ISHAPE(12,2) are not referenced in the above list 
and are not required. 

(10) A single element of an array can be set in the form 
(la), if the array name is given v;ith its numerical 
subscripts; e.g., JSHAPE(3,2) = 1. 

Tlie last card in the input deck for each case, following the 
cards containing the data, must be pxinched* 

$END 


in columns 2-5. 


The Namelist input system processes the ini>uts of the form 
(1) one at a time, as they are encountered in the input deck. 
Thus, the order in which the input variables arc referenced is 
arbitrary. If a varia]Dle is set more than once, the last v^’.lue 
read is used by the program. For example, in figure 1, the pro- 
greuTi v;ov:ld run with CXILA/CI = 45. Inpiit variables whi.ch are 
not referenced in the input deck to a case arc not ciianged; the 
program I'uns with the values already in storage in these loca- 
tions. For example, in the second case of figure 1, only FL0W 


*If NATA is run on an IbM 360, the $ is rcpjaccd by an EBCDIC 
ampersand or a BCD + sign. 
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and TPRNTI are different from the values in the first case. Most 
of the input variables in NATA are preset to values which are 
either usually satisfactory or frequently desired, as indicated 
in Sections 2.2 and 2.3. If these variables are not referenced 
at all in the input to the job, the program runs with the preset 
values. This feature reduces the amount of input data required 
in most NATA runs by orders of magnitude. However, those vari- 
ables which are not preset (such as the reservoir pressure, 

PRESAI) must be set in the input to the first case in every job; 
otherwise, the program would tiry to run with garbage data left 
in the conputer by the preceding job, or with zero in the case of 
a computer in which core is cleared before each job. 

A few of the input variables can be reset by operation of 
the code using certain options. Such exceptions to the rule 
that variables not referenced in the input do not change from 
case to case v/ill be pointed out in Section 2.3. Examples in- 
clude N0ZZLE, ATPI, ISHAPE, NPRFLS, and IGAS. Apart from changes 
in IGAS due to automatic air model selection, these exceptions do 
not arise in jobs containing cases with only a single type of geo- 
metry, i.e., channels or nozzles, 

2.2 Input for Air Cases with a Standard Geometry 

The NATA code contains compiled-in data on the thermochem- 
istry and reaction kinetics of certain gas mixtures (including 
air) and on the geometries of standard nozzles and channels in 
use at NASA/JSC. These precoded gas model and geometry data 
allow the NATA user to rvin certain standard t^^es of problems 
by providing input data for just a few varicibles. The present 
section lists and defines these key inputs, under the following 
assumptions: 

(1) the gas is air; 

(2) the flov; is confined by one of the available precoded 
standard nozzles or channels; and 

(3) the reservoir conditions are to be determined from 
data on the reservoir pressure and the total mass flow. 

Sections 2.3 and 2.4 present a more comprehensive discussion of 
NATA inputs for users desiring to run nonstandard type problems 
or to use some of the special options which give the code its 
flexibility. 
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It is recommended that NATA users employ the input lists 
given below and in Section 2.3 as checklists in setting up 
problems. 

The table below lists the names and definitions of the key 
input variables. For each variable, the "nximbcr of entries" is 
the number of numerical values which can be punched on the right 
hand side of the input equation (1). The preset values are com- 
piled int'v the program and will be used unless a different value 
is supplied in the input deck. The eibbreviations under "Type" 
have the following meanings: 

R - Real (number containing a decimal point) 

I - Integer (number without a decimal point) 


Variable 

Name 

Number 
of Entries 

Preset 

Values 

Type 

Definition 

PRESAI 

1 

— 

R 

Reservoir pressure (atm)» 

FL0W 

1 

— 

R 

Total Mass flow (Ib/sec). 

N0ZZLE 

1 

0 

I 

Index of standard nozzles: 


1 - DCA nozzle with 1.905 cm 

(0.75 inch) throat dia- 
meter 

2 - DCA nozzle with 3.81 cm 

(1.5 inch) throat dia- 
meter 

3 - MRA nozzle with 5.715 cm 

(2.25 inch) throat dia- 
meter 

4 — MRA nozzle with 2.54 cm 

(1 inch) throat diameter 

5 - EOS nozzle with 0.813 cm 

(0.32 inch) tliroat dia- 
meter 



f 


Variable Number Preset 

Name of Entries Value Type Definition 

6 - EOS nozzle with 1.968 cm 

(0.775 inch) throat dia- 
meter 

7 - MRA nozzle with 1.905 cm 

(0.75 inch) throat dia- 
meter 

8 - MRA nozzle with 3.81 cm 

(1.5 inch) throat dia- 
meter 

9-10 MW (Aerothenti) nozzle 

with 5.715 cm (2.25 inch) 
throat diameter 

10 - EOS nozzle with 2.764 cm 
(1.088 i "^h) throat dia- 
meter 

Note: DCA, EOS, MRA, aid 10 MW are designations for electric-arc 

gas heaters in use at the NASA/Johnson Space Center Arc 
Tunnel Facility: 

DCA - Dual-Constrictor Arc 

EOS - Electro-Optical Systems Heater 

MRA - Modified Ring Arc 

10 - Aerotherm Heater 

ICHAN 1 01 Index of standard rectangular 

channels: 

1 - channel with 2.54 x 5.08 cm 

(1x2 inch) throat for DCA 
(use CXMAXI =57.) 

2 - channel with 5.08 x 5.08 cm 

(2x2 inch) throat (nominal 
geometry) for 10 MV-7 heater 
(use CXMAXI = 100.) 
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Note ; Only one of the inputs, N0ZZLE and ICHAN, is used in a 
given case. For ICHAN = 0, the flow geometry is determined by 
N0ZZLE. For ICHAN >0, the input data for N0ZZLE are both ignored 
and overwritten, and the ^low geometiry is determined by ICHAN. 

If an axisymmetric flow problem follows one or more channel prob- 
lems in the same job, it is necessary to input ICHAN = 0 and 
NPRFLS = 1. 


Variable 

Name 


Nuinber 
of Entries 


Preset 

Values Type 


Definition 


CXMAXI 1 1.E5 R Maximum distance beyond throat 

at which free-stream and model 
calculations will be done 
(inches) . if CXMAXI is omit- 
ted, the calculations will 
continue until the free-stream 
temperature drops to 0.4 per- 
cent of the reservoir tem- 
perature . 


TSDIAM(I) 20 1.E20 R For nozzle flow problems, spec- 

ified nozzle diameters at which 
stagnation point of model or 
leading edge of wedge will be 
placed for calculations of mod- 
el test conditions. For chan- 
nel problems, specified chan- 
nel widths at which the free- 
stream flow and conditions on 
the channel wall will be cal- 
culated. Values assumed to 
be in inches. 


KDIM 1 1 I 0 - two-djmensional model 

geometry 

1 - axisymmetric model geo- 
metry 

Note ; The following 5 inputs are needed only if wedge calculations 
are desired. Wedge calculations cannot be obtained in channel flow 
problems. 
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Variable 

Number 

Preset 



Name 

of Entries 

Values 

Tvoe 

Definition 

NANGLE 

1 

0 

I 

NunOaer of angles of attack 
for wedge. 

ANGLE (I) 

10 

0. 

R 

Wedge angles of attack (deg- 
rees) in ascending order. 

NRADLE 

1 

0 

I 

Number of leading edge radii 
for wedge. 

RADLE(I) 

5 

0. 

R 

Wedge leading-edge radil 
( inches) , 

WXI(I) 

20 

1.E30 

R 

Distances from leading edge 
at which conditions on wedge 
will be calculated (inches) . 

ISWIA 

1 

1 

I 

0 suppresses frozen flow sol- 
ution , 

ISW2A 

1 

1 

I 

0 suppresses non equilibrium 
flow solution. 

ISW3A 

1 

1 

I 

0 suppresses equilibrium flow 
solution. 

1SW4A 

1 

0 

I 

Must be nonzero if another 
case follows in the job. Must 
be 0 for last case. 


2.3 General Inputs 

The main inputs to the NATA code are read in under the Name- 
list name INPUT. All of the input variables in this group are 
defined in the present section. Data for new species and reactions 
and transport cross section data are read in under other namelist 
names, as explained in Section 2.4 Problems in which the precoded 
data for species and reactions are used require only the inputs 
discussed in the present section. 
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In the table of definitions below, the ir ut variables are 


arranged 

in ten groups, as follows: 

(1) 

General control variables 

(2) 

Output controls 

(3) 

Reservoir conditions 

(4) 

Geometry 

(5) 

Gas model 

(6) 

Test model conditions 

(7) 

Wedge conditions 

(8) 

Controls for the flow solution 

(9) 

Electronic nonequilibrium 

(10) 

Controls for diagnostic dxmps. 


The table format is the same as in the preceding section, 
except that the array dimensions are listed in place of the "num- 
ber of values" . The type designations R for real and I for integer 
are defined as before; L denotes logical variables, whose admis- 
sible values are T for "true" and F for "false". 


Grov^p 1; General Control Variables 


Variable 

Nan\e 

ISWIA 


ISW2A 


ISW3A 


ISW4A 


Preset 

Dimensions Values 


Type Definition 

I Value 0 suppresses frozen 

solution, 

I Value 0 suppresses nonequi- 

librium solution. 

I Value 0 suppresses equilibrium 

solution, 

I Must be nonzero if another 

case follows in the job. Must 
be 0 in last case. 
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General Control Variables fCont*d) 

Variable Preset 

Name Dimensions Va lues 


ISW6A 1 0 

ditions and transport proper- 
ties in the reservoir are cal- 
culated. For the value 2, the 
reservoir transport property 
calculations are omitted. If 
< 0, tables of speci>^^s rnermal 
properties at temporatn-'* ' i up 
to the reservoir temp are 

are produced and no t ':nl- 

culations are done. 


Type Definition 

I If >0, only the reservoir con- 


Note ; ISW6A<0 should be used together with ISW2B = 1 and input of 
CTAPI and PRESAI (see Group 3 below) . 


ISWIB 1 0 I lf>0, an edit of the steps 

in rhe transport property 
cross section calculations is 
produced before the flow sol- 
utions. If <0, averaged trans- 
port cross sections are also 
printed for temperatures up to 
CTAPI and the flow solutions 
are not computed. If = -1, 
these cross section data are 
also punched on cards. 


Note : For 

(see Group 

ISWIB < 0, 
3 below) 

also set 

ISW2B = 

1 and read in CTAPI and PRESAI 

ISW3B 

1 

1 

I 

If 0, boundary layer on nozzle 
wall is omitted. 

TWALL 

1 

300 

R 

Nozzle wall temperature (°K) 

N0TRAN 

1 

.FALSE. 

L 

If .TRUE,, all transport prop- 
erty, bcjndary layer, heat 
flux, and wedge calculations 
are suppressed. 


r 


( 

i 


Variable 

Name 

TST0PI 


CXMAXI 


General Control Variables (Concl'd ) 


Dimensions 


Preset 

Values 

0 . 


1.E5 


Tyng J3e£iiutiQn 

R Pree~stream temperature at 

which the flow solutions will 
be termirated (°K) . For value 
0., the case is stopped at 0.004 
times the reservoir temperature. 

R Distance beyond the throat at 

which tlie solutions will be 
stopped (inches) . 


Note ; The solutions are stopped by either TST0PI or CXMAXI, which- 
ever condition ^s reached first. 


READG 


READXS 


.FALSE. 


.FALSE. 


If .TRUE,, data on elements, 
species, and/or reactions will 
be read in xinder the namelist 
name E INPUT. 

If .TRUE., cross section data 
for transport property calcu- 
lations will be read in under 
namelist naune TINPUT. 


Group 2: Out - ut Controls 


The inputs in this group allow some user control of the types 
and amount of printed output produced by the code. 


Variable 

Name 

Dimensions 

Preset 

Values 

Tvpe 


Definition 


ISW6B 

1 

1 

I 

Value 0 
species 

suppresses output 
mole fractions in 

of 

free 


stream and mode] point output. 
Positive value gives mole frac- 
tion output every ISWGBth printed 
step. Negative value also gives 
output of reaction 3.*atc data ev- 
ery tlSV76Bl th printed step. 
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Output 

: Controls 

(Concl'd) 

VaxiaCale 

x'lcune 

Dimensions 

Preset 

Valxies 

Type 

Definite, on 

ISW7B 

1 

0 

I 

Value > 0 gives output of he 
boundar*/ layer parameters r=n 
and XSK = n’ . 

TPRNTI 

1 

0.01 

R 

The free-stream nonequilibrium 
solution is printed out at tem- 
perature intervals greater than 
or equal to TPRNTI times the 
reservoir temperature. For 
TPRNTI = 0., evejry step is 
printed. 

DATAPE 

1 

.FALSE. 

L 

If .TRUE., data are vnritten onto 
tape 3 for subsequent platting. 

NRECO 

1 

0 

I 

Nuitber of records already on 
data tape at beginning of run. 

IRUN 

1 

0 

I 

Rvm Number (for identification) . 


Groyi.n> 3: Rcserv^oir Conditions 


Tlie variables in this group control the calculation of the gas 
state in the upstream reservoir. The methods used are explained in 
Section 6.5 of Volume I (ref. 1), 


Variable 

Kar.ie 

ISV72B 


PRrsAI 


Dimensions 

1 


Preset 

Values 

0 


Type Definition 

I If 0, reservoir temperature is 

computed from reservoir pressure 
(PRESAI) and mass flow (FL0W) . 

If positive, reservoir tempera- 
ture (CT2VPI) and pressure (FRES/iI) 
ax'e road in. If negative, reser- 
voir temperature and pressure 
are computed from mass flov; (FL^^V7) 
and stagnation enthalpy (HSTAG) . 

R Reservoir pressure (atm). Re- 

quired in input if ISW2B > 0. 
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Reservoir Conditions fCr>ncl*d) 


Variable 

Name 

Dimensions 

Preset 

Values 

iVpe 

Definition 

FI^W 

1 


R 

Total mass flow (Ib/sec) if 
JDIM = 1 (see under "Geometry"); 
mass flow per inch (Ib/in-sec) if 
JDIM - 0. (FL0W is required in 
the input if ISW2B<0.) 

CTAPI 

1 

— 

R 

Reservoir tenperature (°K) . (Re- 
quired in input if ISW2B > 0.) 

HSTAG 

1 

— 

R 

Stagnation enthalpy (Btu/lb) . 
(Required in input if ISW2B< 0.) 

MFITER 1 

ORIGINAIS PA6B IS 
OF POOR QUAUTB 

1 

I 

0 value suppresses iteration to 
take displacement thickness into 
account in reservoir condition 
calculations based on mass flo^v 
(ISW2B < 0) . 



Grout) 4: 

Geometrv 


The geometry of a nozzle is specified, in NATA, by describing 
the nozzle profile. The geometry of a rectangular channel is spec- 
ified by giving tv;o profiles. There are four optional methods for 
defining the flov/ geometry in the input: 

(1) Standard nozzle - ICUAN raust be 0, NPRFLS raut .. oe 1, 

N0ZZL^ must be an integer in the range from 1 to 10, 
inclusive. 

(2) Standard channel - ICHAl^ must be 1 or 2. 

(3) Nonstandard nozzle - must be 0, NPRFLS must be 1, 

N0ZZLE must be 0, and input data must be provided for 
DIAM(l), NSECTS(L,1), ISHAPE { J, 1) , PARAMI (K, J, 1) , ATP.T(J,1), 
and XZER^ZJl. 

(4) Nonstandard channel - must be 0, NPRFLS must be 2, 

N0ZZLE, i;pR0FL(l), and NPR0FL(2) must be 0, and input data 
must be provided for DI.AM(M) , i.:<=:£CTS (L,M) , ISHAPE (J,M), 
PAIt^M::(K, J,M) , and ATPI(J,H) for M = 1 and 2, and for 
XZER01 and NBL. 

The description of nozzle and channel geometries in NATA is discussed 
in Sections 4.2 and 4.3 of Volume I (ref. 1). 


- 17 - 




Gecxnotrv (Cont'dl 


Variable Preset 

Name Dimensions Values Type Definition 

N0ZZLB 1 01 Index of standard nozzles: 

0 - nonstandard nozzle 

1 - DCA noi.zle with 1.905 can 

{0.75 inch) throat dia- 
meter 

2 - DCA nozzle with 3.81 cm 

(1.5 inch) throat dia- 
meter 

3 - MRA nozzle with 5.715 cm 

(2.25 inch) throat dia- 
meter 

4 - HRA nozzle with 2.54 can 

(1 inch) throat diameter 

5 - EOS nozzle with 0.813 cm 

(0.32 inch) throat diameter 

6 - EOS nozzle with 1.968 cm 

(0.775 inch) throat c7iameter 

7 - MRA nozzle v/ith 1.905 cm 

(0.75 inch) throat diameter 

8 - MRA nozzle with 3.81 cm 

(1.5 inch) throat diameter 

9-10 MW (Aerotherra) nozzle 
with 5.715 cm (2.25 inch) 
throat da meter 

10 - EOS nozzle with 2.764 cm 

(1.0S8 inch) throat diametei. 

Notes ; N0ZZLE is altered by input of NPR0Fi, or by execution of a 
case involving ICI1AN>0. In cases with M0ZZLE = 0, the firsu 4 char- 
acters on the description card at the head of the input data are used 
as a facility name. 
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Geometry (Cont*d) 

Variable Preset 

Name Dimensions Values Types Definition 

ICHAN 1 01 Index of standard rectangular 

chaimels: 

0 - not a channel, or nonstan- 

dard channel 

1 - channel with 2.54 x 5.08 cm 

(1x2 inch) throat for 
DCA. (Use CXI4AXI =57.) 

2 - channel with 5.08 x 5.08 cm 

(2x2 inch) throat (nomi- 
nal gecxaetry) far 10 MW 
heater. (Use CXMAXI = 100.) 

Note ; In channel cases with ICHAN = 0, the second 4 characters on 
the description card at the head of the input data are used as a 
channel name. 


JDIM 1 1 I 0 - two-dimensional nozzle 

1 - axisymmetric nozzle 

Not e; A tv7o-dimensional noi^-zle may be considered as the limiting 
case of a rectangular channel when one of the channel profiles is 
at an infinite distance from the channel axis. Txiis limit is of 
little practical interest. However, JDIM = 0 gives a convenient 
way of treating the flov? in a channel in which two of the walls 
are straight and parallel, when the boundary 'ayer is neglected 
(ISW3B = 0) . 


NPRFLS 1 II Number of profiles: 

1 - nozzle 

2 - rectangular channel 

Note ; The program sets NPRFLS =2 if ICHAN >0. 
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Geometry (Cont’dl 


Variable Preset 

Name Dimensions Values 

NPR0FL(I) 2 0 


Definition 

Indices of profiles in a chan- 
nel; NPP0FL(1) is equivalent 
to N0ZZLE. 

0 - nonstandard profile 


1 to 10 - profiles for standard noz- 
zles (see above) 

11 - profile 1 for T12 and T22 
channe Is 


DIAM(M) 


NSECTS(L,M) 2x2 


12 - profile 2 for T12 channel 

(ICHAN = 1) 

13 - profile 2 for T22 channel 

(ICHAN = 2) 

Index (1 or 2) of the profile 
which diverges from the channel 
axis least rapidly downstream 
of the throat. 

For a nozzle, DIAM(l) is the 
throat dieimeter (inches). For 
a channel, DIAM(M) is the throat 
diameter of the Mth profile for 
M = 1 and 2 (inches). 

NSECTS(1,M) = number of upstream 
sections in curvefit for Mth pro- 
file; NSECTS(2,M) = number of 
downstream sections in curvefit 
for Mth profile. For a nozzle, 

M = 1; for a channel, M = 1,2. 


- 20 - 



r 


Geometry (Cont'd ) 

Preset 

Values Typ e Definition 

— I Shape index for Jth section of 

Mth profile: 

ISHAPE = 1 straight section 

ISHAPE = 2 circular section 
convex toward axis 

irHAPE = 3 circular section 
concave toward axis 

Note ; ISHAPE is altered by use of N0ZZLE > 0 or ICHAN>0. 

PARAMI(K. J,M) 3 x 12 X 2 R Parameters for profile curvefit 

sections (lengths in centimeter 
units) : 

For ISHAPE (J,M) = 1, equation of 
straight profile is r(x) = 
PARAMI(1,J,M) + PARAMI(2, J,M) *X 
For ISHAPE (J,M) = 2 or 3, 

PARAMI (1, J,M) = distance of cir- 
cle center from axis 
PARAMI (2, J,M) = X coordinate of 
circle center 

PARAI>'II(3, J,M) = circle radius 

ATPI(J,M) 11 X 2 R Do\\mstreara boundaries of pro- 

file curvefit sections, measured 
from throat (cm) . 

Note ; ATPI is altered by use of N0ZZLG > 0 or ICH/\N >0. 


Variable 

Ntune Dimensions 


ISHAPE (J,M) 12 X 2 


XZER0I 1 R Nozzle or channel inlet position 

at v;hich boundary layer is as- 
sumed to begin (negative value, 
measured in inches upstream from 
th . th-. oat) . 
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Geomolr y^ fConcl'd) 

Note : XZER0I is altered by use of N0ZZLE > 0 or ICIIAN>0. 

Note : A separate progreun (NOZFIT) is available for computing 

the inputs PARAMI and ATPI from nozzle or channel design data 
such as dimensions, angles, and radii of curvature; see Appen- 
dix b. 


Group 5: Gas Model 

NATA provides three methods for input specifications of the 
con.positi an , thermochemistry cuid kinetics of the gas mixture: 

(1) Standard gas models • invoked simply by setting IGAS 
to 1, 2, 3, 4, 5, or 6. 

(2) Standard gas models \ ith altered elemental composition - 

obtained by setting IGAS = -1, -2, -5, or -6 and specify- 
ing the mole fractions of the cold species (QPJ; see 
Section 4.5). Using this option, the standard air models 
can be modified to obtain models for (nearly pure) oxygen 
or nitrogen by setting the rn^le fraction for the other 
cold species (nitrogen or oxygen) to a small value. The 
proportions of CO 2 , N2» planetary atmos- 

phere models can also be changed in this way. 

(3) Nonstandard gas model s - specified by setting IGAS - 0 
and reading in NCS, JCS, OpJ, ISCI, ISSI, ISRI, ICI, IE, 
IS, IR, ISAT0M, and ISM0L. If sp'^cies or reactions other 
than those compiled inco the code are desired, they can 
also be read in as explained in Section 2.4. 

Tlie compiled-in species, reactions, and gas models availab].e in NATA 
are fully described in Section 4. 


Variable 

Name 


Preset 

Dimensions Values Tyre 


Definition 


IGAS 


1 


1 


I Gas model index; 


0 - nonstandard gas mixture; 

NCS, JCS, QPJ, ISCI, ISSI, 
ISRI, ICI, IE, IS, IR, ISAT0;>>. 
and ISM0L must all be speci- 
fied to the input. 

1 - high- temperature air model 


- moderate-temperature air 
model 




ORIGINAD PAGB H 
OF POOR QUAIOT 
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Gas Model 

(Conf d) 



Variable 

Preset 




Name 

Dimensions Values 

Type 

Definition 



3 - argon model including elec- 

tronic nonequilibrixim 

4 - helium model including el- 

ectronic nonequilibrium 

5 - planetary atmosphere model 

(75% CO^, 20% hr, 5% N 2 ) 
for use at reservoir tem- 
peratures above 7000°K. 

6 - planetary atmosphere model 

(75% CO 2 , 20% hr, 5% N 2 ) 
for use at reservoir tem- 
perature below 7000°K 

Note ; If a negative value of IGAS is specified, then llGASi is 
the index of a standard gas mixture for which the mole fractions 
of cold species (QPJ) are to be provided in the input. 


AAMS 

1 

.TRUE. 

L 

Control for automatic air model 
selection. If IGAS = 1 or 2 and 
AAMS = .TRUE., NATA resets IGAS 
to 1 or 2 based on an enthalpy 
or temperature criterion; for 
AAMS = .FALSE., the IGAS value 
specified in the input is used. 

NCS 

1 

— 

I 

Number of cold species in mix- 
ture ( < 10) , 

JCS(I) 

10 

— 

I 

Indices of cold species in the 
master list of species*. 

QPJ (I) 

10 

— — 

R 

Mole fractions of cold species 
in the seme order as JCS (must 
be provided if IGAS < 0) . 


*Seo Section 4,2. 
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Gets Model fCont:‘cn 


Variable 

Name 

ISCI 


ISSI 


Preset 
Dimensions Values 


ISRI 


ICI 


IE(I) 


Definition 

I Number of chemical elements in 
mixture, including e“ if model 
contains ion species (< 10), 

I Number of chemical species in 
mixture, including e“ if model 
contains ion species (<20), 

I Number of reactions included 
in gas model ( <64). 

I Ntimber of ions in gas model 
excluding e“. 


10 — - I Indices of elements present in 

mixture, in master Ij.st of ele- 
ments;* if electrons are present 
they should be the first element 

20 — I Indices of species present in 

mixture, in master list of spec- 
ies. **These species must be 
listed in the following order; 
e“ (if present) 

Neutral species which are 
stable at low temperatures 

Other neutral species 

Ion species 


Note ; The first ISCI species in this list must be linearly inde- 
pendent combinations of the ISCI chemical elements. 


IS(J) 


IR 


64 


Indices of reactions included, 
in master list of reactions.*** 


*Sco Section 4.1. 

**Sec Section 4.2. 

***Scc Section 4.3. 
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Gas Model fConcl«dl 

Note ; If a nonequilibrium flov; solution is to be run with the 
read-in gas model, there must be (ISSI-ISCI) linearly independent 
reactions in the chemical kinetic model; see Section 7.3.4 of 
Volume I (ref. 1). The standard gas models all satisfy this re- 
quirement. 


Variable 

Name 

Dimensions 

Preset 

Values 

Type 

Definition 

ISAT0M 

1 

— 

I 

Index of atom used for Lewis 
number calculations, in master 
list of species. 

ISM0L 

1 


I 

Index of molecule used for 
Lewis number calculation, in 
master list of species. 

CTMXXI 

1 

5000. 

R 

Temperature (°K) above which 
species thermal properties are 
computed from the thermo fit for 
those species for which thermo 
fits are supplied. 

BZER0I 

1 

0.0 

R 

Constant in imperfect gas cor- 
rection; the 0 value suppresses 
the correction, which is negli- 
gible for the conditions in 
v;hich NATA is normally applied. 

INEQVI 

1 

0 

I 

0 - equilibrium molecular vib- 


ration 


1 - molecular vibration frozen 
at the reservoir tempera- 
ture 


Group 6: Test Model 


NATA provides calculations of test conditions on two types of 
models: blunt bodies (stagncition point only) and wedges. A sin- 

gle set of inputs (XM0DP1, NM0DPT, TSDliVM) controls the positions 
in the flov; at which test conditions are calculated for both types 
of model. These inputs and the parameters controlling options in 
the calculations for blunt models are in the present group. The 
wedge model inputs arc in group 7. 
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Test Model (Cont'd) 

In the case of blunt models, the test model position determined 
by the inputs is assumed to be the location of the model stagnation 
point. In the case of wedge models, it is assvimed to be the loca- 
tion of the leadjjig edge. There are two options for specifying 
the test model positions: 

(1) A geometric sequence of distances x downstream of 
the throat from x = XM0DP1 to x = CXMAXI. 

(2) The positions at which the nozzle diameter is equal 
to the input values TSDIAM(I) . 

These two options operate independently. In channel flow sol- 
utions, model condition ca:i culations are not done but extra points 
in the free streami solutions are inserted at the locations speci- 
fied by XM0DP1, NM0DPT, and TSDIAM, to provide results for compar- 
ison with experimental data from pressure taps and heat transfer 
gages located at known positions on the channel wall. 

Regardless of the model-position inputs, no model or wedge 
calculations are done at positions where the flov/ Mach number 
is less than 1.5. 


Variable 

Name 

Dimensions 

Preset 

Values 

Tvpa 

Definition 

XM0DP1 

1 

1.E20 

R 

Initial x for model condition 
calculations, measured in inchei 
dov/nstreani of the throat. 

NM0DPT 

1 

20 

I 

Number of model points to be 
placed in a geometric progres- 
sion from XM0DP1 to CXMAXI; for 
NM0DPT - 1, the model calcula- 
tion is done at x = X’I0DP1. 

TSDI7^i(l) 

20 

1.E20 

R 

In nozzle flow problems, nozzle 
diameters specifying model pos- 
itions; in channel problems, 
specified channel widths at 
which extra flow calculations 
are done (inches). 

AXIM0D 

1 

. TRUE . 

L 

Value .FALSE, suppresses stag- 
nation j3oint model condition 
calculations, if only wedge 
conditions are desired. 
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Test Modol (Cont:*d) 


Varia33le Preset 

Name Dimensions Values Type Definition 


KDIM 


FSTAG 


1 


1 


1 I Indicator for blvmt model 

geometry : 

0 -• t\i70“dimensional 

1 - axi symmetric 

1.0 R Control for normal shock cal- 

culations: 


CATFAC 


LEI^IS 


1 


1 


Negative - frozen shock 

Zero - equilibrium shock 

Positive - both frozen and 
equilibrium shock computed 

1,0 Catalytic wall parameter § for 

stagnation point heat flux cal- 
culations;* 0 for noncatalytic 
modeJ. surface; tiieoretical up- 
per limit is <1. 

1 I Indicator for use of Fay-Riddell 

Lewis nur.fr-er factor in stagna- 
tion point heat flux calcula- 
tions: 

1 - include Lewis no. factor 

2 - omiu Lewis no. factor 


Note ; LEWIS is not vmder input control v;hon a standard gas model 
is used. 


TM0DEL 1 300. R Surface temperature at stagna- 

tion point on model (*^K). 


*See Section 8.1.2 of Volume I (ref. 1) . 
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Test Model 

(Concl'd.) 

Variable 

Preset 


Name 

Dimensions Values 

Type 


Definition 


TPLATE 1 300, R Flat plate temperature for cal- 

culations of heat flux to a flat 
plate 1 ft from loading edge (°K), 


Group 7; Wedge Models 

Wedge model calculations are done only if the following condi- 
tions are satisfied: 

(a) Model positions have been specified by input of 

Xfl^JDPl and or of TSDTT'i'I(I) ; 

(b) The flow Mai'h nuniber at the specified model posi- 
tions is greater than 1.5; 

(c) Positive values have been specified for both NTfNGLE 
and NRADLE ; and 

(d) Either N1'K>0, or a value htis been specified for 
^7X1(1) . 

The positions along the surface of the wedge at v/hich the con- 
ditions are calculated can be specified in two ways: 

(1) Uniform secfiience - The inputs V7X1, D^\!X, and NTiJX deter- 
mine a uniformly spaced sequence of distances from the 
leading edqe. 

(2) -r '-'i fieri di stances - up to 20 arbitrary distances 
from the leading edge can be specified using the input 
array IfXI, 

Both options may Vjg vised together, if desired. In all cases, the 
specified distances from the leading edge are measured along the 
surface of the wedge (rather than parallel to the direction of the 
incident fiov;.) 
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o Models (Cont'd) 


Variable 

Name 

Dimensions 

Preset 

Values 

. Tvp 

c Definition 

NANGLE 

1 

0 

I 

Nun^er of wedge angles of attuc]. . 

ANGLE (I) 

10 

o 

• 

o 

R 

Angles of attack of wedge surface 
relative to the direction of in~ 
cident flow (degrees) , in asccijo.- 
ing order. 

NRADLE 

1 

0 

I 

Nvuiibcr of leading edge radii . 

RADLE(J) 

5 

o 

• 

o 

R 

Radii of leading edge finches) . 


1 

1.0 

R 

Distance of the first co^aput^’.t;i o’i 
point from the leading edge (in- 
ches) . 

DWX 

1 

1.0 

R 

Distance betv;een computati.on 
points (inches) . 

hWA 

1 

0 

I 

Number of compi . points in 

uniform soquen«.. . 

v;xi(i) 

20 

1.E30 

R 

Specified distances o.c couipii'-a- 
tion points froiti 3oadini.r edge 
(inches) 

T^ffiDGE 

1 

300. 

R 

Wedge surface temperaturu (®J'.) . 

WK 

1 

1.333 

R 

Nose drag coo.rf iciont for Ci’cnc;-- 
Kemp wedge theory; the preot.'t 
value is for a cylindrical lead- 
ing cxige. 

ISV79B 

1 

0 

I 

Control for \,'odgc modfjl c'r IcuLe- 
tions aiH'i out]')Ut*. 





•to - inc.ludo caleulatio.n;; 

unmodi.'i.cd Cjicng-lt ' p tirecr: 





iO - omit uni'.’odi.[icd tioo!-' * 

*Scu Soc'r..i 

.on V.2.A C)Z 

Vo].uiuo I 

(rci!. 

1). 
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Wedge Moclols fConcl*d) 


Variable Preset 

Nome Dimensions Value s Type Definition 

With IS9 = IISW9bI : 

IS9 - 1 — print shock ordinate 
Ys 

IS9 = 2 — print nondimensional 
distance C from leac 
ing edge 

IS9 = 3 - print both Yg and X 


Group 8: Controls for the Flow Solutions 

The inputs in this group are control parameters for the frozen, 
equilibriuni and nonequilibrixuti flow solutions. They are all preset 
to values v.’hich have proven satisfactory in practice, and need be 
varied only rarely, to treat cases in which the code has failed to 
produce a satisfactory solution when run v;ith the standard values. 

Variable Preset 

Name D5-ir.cn si on s Value s T\T>e Definition 


WSAVE 1 3.0 R 


DELTII 1 0.01 


Parameter controlling the averag' 
ing distance for the boundary 
layer correlation pariimeter, n. 
Instability due to coupling of 
the inviscid flow with the boun- 
dary layer can be suppressed by 
reducing V7SAVE.* 

R Nondimensional teirperature dec- 
rement used in fronen and equi- 
librivun calculations and in 
starting the nonequilibrirua sol- 
ution. 


DELTXI 1 0.01 R Initial step size in X for non- 

cqi’d librium integration (cm) 
(may be overruled by code) . 


*Gcc SecLion 5.11 of Voljmc I (ref. 1) . 


ORIGINAL PAGE IS 
OF POOR QUAIOT 
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Variable 

Name 


CCHI 


NQSI 


TTEST 


GTEST 


HTEST 


TETEST 


QTEST 


GAMIN 


Controls for the Flow Solutions (Cont'd) 


Preset 

Dimensions Values Type 


Definition 


1 


1 


1 


1 


1 


1 


1 


1 


0.1 R Criterion value for switch 

from pertxirbation technique to 
numerical integration in non- 
equilibrium solution; increase 
to switch farther dovmstream, 
decrease to switch farther up- 
stream,* 


4 


0.05 


I Number of successful integration 
steps before increasing step 
size in the nonequilibri\jm cal- 
culation. 

R Maximvun |aT/t1 in one step of 
the nonequilibrium integration; 
decrease to force a smaller step 
size. 


0,1 


0.01 


0.05 


R Maxiraxm re lative species concen- 
tration change in one integra- 
tion step; decrease to force a 
smaller step size, 

R Maxinnm relative change in the 

total enthalpy (due to radiative 
losses) in an integration step. 

R Maximum relative change in the 
electron temperature in an in- 
tegration step. 


0.1 R Criterion value for maximum 

allov70}Dle change in the energy 
transfer to the electron gas 
during one integration step, 

10~^^ R Concentration (moles/g) belov; 

which a species will be frozen 
if it decreases so rapidly that 
it controls the integration step 
size. 


*Sce Section 7.3.6 of Volume I (ref. 1) . 
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Controls for the Floy; Solutions (Co -I'd) 


Variable 


Preset 



Name 

Dimensions 

Values 

Type 

Definition 

DCHLL 

1 

i 

o 

R 

Parameter limiting the initial 

integration step size to 

0.01 . /DCHLL. 

X raxn 

DCHRAT 

1 

10-4 

R 

Parameter controlling the arti- 


ficial increase in reaction rates 
in the perturbation solution to 
avoid premature startup of the 
nvimerical integration; minimum 
allowable l^^ilmax 

value. 


Group 9: Electronic Nonequilibrium 

The standard models for helixam and argon include electronic 
nonequi librium effects such as inequality of the electron tempera- 
ture and gas temperature and nonequilibrivuti population of electronic 
excited states. KATA allov/s nonstandard gas models containing t lese 
features to be set up by the user. The inputs in the present group 
provide the extra gas model data required to specify these effects. 

Variable Preset 

Name Dimensions Values Type Definit ion 

INT 1 01 Indicator for electronic non- 

equilibrium: 

0 - conventional one-temperature 
gas model 

Nonzero - tv;o- temperature (elec- 
tronic nonequilibrium) 
model 

KTF(IR) 25 I Indicator for forward rate con- 

stant kf for IRth reaction in 
gas model 
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Variable 

Name 


Electronic Nonequilibrium (Cont'd) 


Preset 

Dimensions Values Type 


Definition 


1 - kf = kf(T) 

2 - kf = kf (Tg) 


where T = gas temperature, 

Tg - electron temperature, 
anil the functional dependence 
of kf is as given by equation 
(69) in Section 2.3 of Volume 
I (ref. 1) 


3 

4 


. . A (l-e-=aAoT) 

^ 10^ 

- kf = A (2^)^ /max (1, T ) , 
^ 104 


where t = b UpR/Ng 


5 


- kf = A//R 


Note ; In the standard gas models (IGAS = 1 to 6), the rate formu- 
las indicated by KIF - 3, 4, and 5 are used only in the argon model 
(IGAS =3). Note that for KTF =3, kf depends on both and T. In 
the formulas for KTF =4, 5, R denotes the local nozzle radius (or 
a corresponding effective value in the case of a channel) . Also, 

Up is the number density of the atomic species appearing on the pro- 
duct side of the reaction. See Appendix A for a discussion of these 
rate formulas. 


I-?rR(IR) 25 T Indicator for reverse rate con 

stant kj- for IRth reaction in 
gas model: 

0 - kj- = 0 

1 - kj. = kj.(T) 

2 - kj. = kr(Te) 
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Electronic Monecrui librium (Cont ' dj 


Variable 

Name 


Dimensions 


Preset 
Va lues 


Definition 


ITR(IR) 


Indicator of rule for partition- 
ing the reaction energy of the 
IRth reaction; values may be 1 
to 6. 


Note ; In the definitions belov/, Sf and - 6^, denote the energies 
gained by the electron gas in Nq reactions in the forward and re- 
verse directions, respectively, and q^, -q^ denote the correspond- 
ing energies lost by radiation. Also, Nq = Avogadro's number. The 
admissible values of ITR correspond to the following relations:* 


ITR = 1 


- a R^T^, 


qf = 


ITR = 2 


ITR = 3 


= < 1 ^ = 0 
= - f 






= e = q = q = 0 


ITR = 4 


ITR = 5 


ITR = 6 


Sf = = - t 


qf =•- qr = 0 

■If “ <Jr = ° 

q = e 
f 0 


Tlie application of these reaction energy partition rules to reactions 
in argon and helium is discussed in Appendix A. 


*The formulas for ITR = 2 arc a special case of those for ITR - 1, 
The reasons for tlxl s formulation are historical rcatlior than logical. 
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Electronic Nonequilibrium (Concl'cO 


Variable 

Name 

Diniensiona 

Preset 

Values 

Type 

Definition 

EPAR(I,IR) 

2 X 25 


R 

EPAR(1,IR) = parameter for 


the IRth reaction in cal per Nq 
reactions; EPAR(2,IR) = para~ 
meter a for the IRth reaction if 
ITR(IR) = 1. 


BPAR 

1 


R 

Parameter b for all reactions 
v;ith KTF = 4. 

TLIST(J) 

30 


R 

Temperature values for teible of 
elastic collision cross section 

P0M(J) 

30 

— 

R 

Elastic collision cross section 
values q( 1|1) table (cm-) 


Group 10: Controls for Diagnostic Dumps 

NATA contains a nturiber of coded-in provisions for special out- 
put to facilitate tracing the operation of certain sections of the 
program These diagnostic dumps are intended for use by programmers 
in analyzing causes of code failure. Ordinary users of NATA will 
rarely find occasion to invoke them. The input variables control- 
ling these diagnostic outputs are defined below. 


Variable 

Name 

Dimensions 

Preset 

Values 

Tvpe 

Definiticn 

ISW5A 

1 

0 

I 

If nonzero, the execn.tion of 
subroiitine RESTMP is traced by 
dumps. 

ISW4B 

1 

0 

I 

If >0, a large dump is written 
each time the boundary layer 
routine BLAYER is called; if 
<0, a one- line dump is v/ritten. 

ISVJ5B 

1 

0 

I 

If >0, largo dumps arc v;rittcn 
each time the subroutines C0MII, 
EXACT, RNICT, PRTA are called and 
at a point in subroutine N0irEQ. 
If <0, these diimps arc written 
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Controls 

for Diacfnostic 

Dumps fCont'd) 

Variable 

N^u:ne 

Dimensions 

Preset 

Values 

Tme 

Definition 

ISW5B 

(Cont'd) 




every IlSWSBlth entry into 
C0MI'I, and a one- or two- line 
dump is v/ritten by NONEQ in every 
step. 

ISW8B 

1 

0 

I 

If nonzero, diagnostic dumps are 
printed in the transport prop- 
erty routines. For ISV78B> 0, 
PUTQIN dump is produced every 
ISW3B times. For ISW8B < 0, 
PUTQIN dvmip is suppressed. 


2.4 Input of Gas Species and Reactions 

Elements, chemical species, and reactions other than those avail- 
able in the precoded data can be defined and used in NATA flov/ cal- 
culations. These additional data are read in under the namelist 
neanc EINPUT. If such data are to be provided, the input variable 
READG in the main input must be set to .TRUE,. Then, immediately 
following the $EI>ID card of the main input, there must be a card 
containing 


$EINPUT 

in columns 2-3. This card is followed by the inpur data, discussed 
below, in the namelist format described in Section 2.1, and tlie 
input data cards must be follov/ed by another card with $END in 
columns 2-5. 

The input variables for defining chemical elements are NEELS, 
a ten-entry array lEEP(I), and ten two-entry arrays EEPl, EEP2,..., 
EEPlO. These variables are defined in the follov.dng table. None 
of them are preset. 
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Variable 




Name 

Dimensions 

Tvoe 

Definition 

NEELS 

1 

1 

Number of elements being defined ( ^ 10) . 

lEEP(I) 

10 

1 

Indices assigned to the defined ele- 
ments in the master list of elements. 

EEPn (J) 

2 

R 

Data for nth defined element 

J = 1 atomic number 
J - 2 atomic weight (g/mole) 


The elements available in the precoded data and their assigned pos- 
itions in the master list of elements are specified in Section 4.1. 

The input variables for defining chemical species are SPl, 

SP2, SP30, each of which is a 43-entry array. The nvimber n 

in the array name SPn is the index assigned to the species in the 
master list. The available species, their properties, and their 
locations in the master list are all specified in Section 4.2. 

Data for any of the standard species can be changed for a particu- 
lar NATA mm by reading in the SPn array used to store its proper- 
ties. 

The data in the SPn arrays are defined and discussed in detail 
in Section 4.2. The definitions of the SPn array entries are sum- 
marized briefly below for convenient reference. All entries are 
real, but those with integer values may be pvmched without decimal 
points, as the Namelist input system will supply the decimals and 
NATA provides for reliable rounding down to the correct integer val- 
ues in cases where this is required. 

SPn(l) Read in 0. (Contains species name in compiled-in 

data; the code supplies a name for identifying the 
species in the output.) 

SPn(2) Number of chemical elements in species (43). 

SPn (3-5) Indices of elements in the master list of elements 

(as modified by input data for elements) . 

SPn (6-8) Numbers of atoms of elements in a molecule of the 

species. 

SPn (9) Thermo-fit coefficient* n. 

*See Section 2.2 of Volume I (ref. 1). 
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SPn(lO) 

SPn(ll) 

SPn(12) 

SPn(13) 

SPn(14) 

SPn(15) 

SPn(16) 

SPn(17) 

SPn(18) 

SPn(19) 

SPn(20) 

SPn(21-30) 

SPn(31-40) 

SPn (41-43) 


Thermo- fit coefficient b. 

Thermo-fit coefficient c. 

Thermo-fit coefficient d. 

Tliermo-fit coefficient e. 

Thermo-fit coefficient k. 

Formation enthalpy at 0®K (cal/mole) . 

N\ainber of atcxns per molecule*. 

Chemical constant,** b. 

Characteristic vibrational temperature (°K) . 

Nxiraber of electronic levels ( < 10) . 

1 if t.ermo fit data are used for species, 0 if not. 

Degeneracies of the electronic levels. 

Energies of the electronic levels (cal/mole) , 

Vibrational temperatures for the second, third, and 
fourth vibrational modes (triatomic species only) 
(°K). 


When a new species model (as defined by an SPn array) is first 
set up, it is advisable to make a preliminary rvin with ISW6A =~1 
(Section 2.3, Group 1) to print a table of species properties as 
calculated from the model. Errors in the species inputs can be 
detected more readily in such a table than in the results of flow 
calculations. 


*Input of SPn (16) = 0 suppresses all use of the "physical model" 
for calculating the thermal properties of the species. In this 
case, the properties are calculated from the thermo fit at all 
temperatures, and SPn(17-19) and SPn(21-43) are not used; SPn(20) 
must be equal to 1 in this case, 

**See Section 2.2 of Volume I (ref. 1). 
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The input variables for defining reactions are RPl, RP2,.,., 
RP64, each of which is a 29-entry array. The entries are defined, 
and the available compiled-in reactions are specified, in Section 
4.3. The definitions are repeated here for ease of leterence. 

All of the entries eire real, but those with integer values may 
be punched without the decimal point. 

RPn(l) Constant A in rate equation (sec“^, cm^/mole-sec, 

or cm^/mole^-sec) , 

RPn(2) Ejqjonent in rate constant formula.* 

RPn(3) Activation energy E in rate constant formula.* 

a 

RPn(4) 1.0 if a list of third-body species is provided 

in RPn (20-29); 0.0 if not 

RPn(5) NxMtiber of reactant species ( < 3) . 

RPn(6) Nvunber of product species ( <3). 

RPn (7-9) Indices of reactant species in the master list of 
species, as modified by the input data for species 
(if any) . 

RPn (10- 12) Indices of product species in the master list of 
species. 

RPn (13 -15) Numbe3:s of molecules of reactants. 

RPn(16-l8) Nun±)ers of molecules of products. 

RPn(19) Nxmber of third bodies ( llO) . 

RPn (20-29) Indices of third body species in master list of 
species . 

If any transport property calculations are to be done for 
new species read in under EINPUT, it is also necessary to provide 
transport cross section data for the species. However, if only 

*Equation (69) in Section 2.3 of Volume I (ref, 1); see also the 
definition of KTF(IR) under Group 9 in Section 2.3 above. 
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inviscid flow calculations are desired, the code can be run with- 
out cross section data by setting the control variable N0TRAN 
to .TRUE, in the main input (Section 2.3, Group 1); this suppres- 
ses all transport property, boundary layer, heat flux, and wedge 
calculations everywhere in the code. 


If transport property calculations involving a new spe<’:ies 
are required, the v-iriabT-e READXS in the main input must be set 
to .TRUE,, and the cross section data for the ^=pecies are then 
read in under the namelist name TINPUT, The input cards contain- 
ing these data immediately follow the deck of cards read vmder 
the name EINPUT. They begin with a card containing 

$TINPUT 

in columns 2-8, This card is followed by the cards containing 
the cross section input data in the namelist format described in 
Section 2.1. The final card of this group must contain $END in 
colximns 2-5. 

The transport property. inputs are as follows: 

Variable 

Name Dimensions Type Definition 


NNKQ 


KKQ(N) 


NNQ(N) 


1 


100 


100 


I Nximber of steps in the cross section 
calculation for which data are speci- 
fied (including compiled-in data) . 

I Index specifying the option to be used 
in the Nth step of the cross section 
calculations (allowed values, 2 through 
14). NNKQ values are required. (The 
meaning of each of the allowed KKQ val- 
ues is given in Section 4.6). 

I Nximber of pairs of species to which the 
cross sections calculated in the Nth 
step are to be applied. NNKQ values 
are required. (NNQ(N) <5). 
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Variable 

Name Dimensions Type 


Definitions 


In(K) 5 

Jn(K) 5 


Vn(K) 5 


ISEU(L) 100 


TL 

1000 

0MEGA1 

1000 

ASTAR 

1000 

BSTAR 

1000 


I Indices of the species to which the 
cross sections calculated in the Nth 
I step are to be applied, referred to the 
master species list. In these variable 
names, n denotes an integer (equal to 
N) which is part of each name. There 
are 100 arrays of each type, namely 
I1(K), I2(K)‘..., 1100 (K), J1(K)..., 
JlOO (K) . There are NNQ (N) pairs of 
indices In, Jn for each step N, Only 
pairs with In(K) < Jn(K) are used in 
the transport property calculations. 

R List of input parameters for the Nth 

step of the cross section calculations. 
There are 100 of these arrays, V1(K), 
V2(K),.,,, VlOO(K). The number of par- 
ameters required for each option and 
their definitions are discussed in Sec- 
tion 4.6. 

I Sequencing array for specifying the or- 
der in which the defined steps will be 
carried out during the cross section 
calculations. The index N in the pre- 
ceding arrays is given by N = ISEQ(L) 
where L = 1, 2, 3, 4,,,., NNKQ is the 
order in which the steps are executed.* 

R Additional storage locations for cross 
R section data. The use of these arrays 
R is discussed in Section 4,6, 

R 


Whenever a new set of cross section data is used for the first 
time, it is advisable to check these inputs by making a preliminary 
run v;ith ISWlB set equal to -2 in the main input to invoke the com- 
plete edit of cross sections. Tliis edit consists of three parts. 


*Tlais input allows steps to be added to the cross section calcu- 
lation (e.g., for new species) without shifting any of the compilcd- 
in data in the li, Ji, and Vi arrays. 
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The first part lists all of the defined steps in the cross sec- 
tion calculation, incl\;ding these compiled in for computing 
the transport properties of the standard species and any steps 
which have been added by input. The steps are listed in th 
order in which they would be performed if the current gas model 
were to include all of the standard and defined species. The 
second part of the edit lists the steps selected by the transport 
routines for the current gas model. This list omits steps which 
are required only for calculating the cross sections of species 
which are not present in the gas model, and includes steps which 
have been added by the default options. This second part of the 
edit thus shows how the transport properties will actually be 
calculated in the current problem. The third part of the edit 
is a set of tables giving the cross sections for each pair of 
species in the current gas model as a function of temperature up 
to the input reservoir temperature, CTAPI, 

The simplest method for specifying the cross sections for 
pairs involving a new species is to rely, to the maximum extent 
possible, upon the default options in the NATA transport property 
routines, T!ie inputs required may then be summarized briefly as 
follows: 

(1) If the new species is an ion, no cross section in- 
puts for it are required, provided the species is 
assigned a previously unoccupied location in the 
master list of species, 

(2) If the new species is neutral, one step must be 
added to the transport cross section calculations 
to compute the cross sections for interaction of 
the species with itself. Alternatively, the like- 
like interaction for the new species can be added 
to an existing step in the cross section calcula- 
tion by increasing NNQ for the step by 1 and adding 
the new species to the corresponding In and Jn lists. 

In the absence of other specifications, all of the 
unlike pair exoss sections involving the new species 
will automatic ally be calculated using a mixing rule 
(option KKQ = 10) . 

(3) Adding a step to the calculations requires the fol- 
lowing changes in the transport inputs: 

a, NNKQ must be incx'cased by 1 (see transport 
block data routine for the original value). 
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b. The option to he used in calculating the 
cross sections for ihe like-like interac- 
tion of the new t^^ecies with itself must 
be specified in the form 

KKQ (n) = option nijuber 

where n is the numerical inder (n) of a 
step which is not already used by the 
compiled-in data. 

c. The number of pairs of species to which the 
new step is applied is set to unity; 

NNQ(n) = 1 


d. The indices of the species are set: 

In = index assigned to new species 
Jn = same index 


e. The parameters for the option are set 
Vn = list of values 


f. The ISEQ array must be modified to insert 
the new step ahead of the step in which the 
mixing rule (option 10) is exercised. The 
proper location can be determined by examin- 
ing the KKQ and ISEQ arrays in transport block 
data. This positioning of the nev; step is re- 
quired to allow’ application of the mixing rule 
to the unlike pairs involving the nev; species. 

g. If the option selected required addition of 
data to the TL, OMEGAl, ASTAR, and BSTAR 
arrays, these data must be sot, Tliis can 
be done most conveniently by reading in each 
array entry as a subscripted variable; e.g., 

TL(k) = value, TL(k + 1) = value, etc. 
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where k, k + 1, etc. are numerical index 
values of array elements which are not used 
by the compiled-in data. 

If information is available concerning the cross sections 
for some of the unlike interactions involving the new species, 
the acctiracy of the NATA tremsport property calculations can be 
im^' roved by adding additional steps to calculate these cross sec- 
tions. These steps, like the step for the like— like cross sec- 
tions, should be inserted ahead of the step in which option 10 is 
applied. 


2,5 Execution Time 

An estimate of execution time is normally required whenever 
a job is submitted for running on a computer system. The per-case 
execution time of NATA can vary greatly depending on the types of 
solutions requested and the gas models used. Typical times for 
the various parts of the calculation are listed below for runs 
using the standard air models (IGAS = 1, 2) . The times are roughly 
similar on the UNIVAC 1108 and the IBM 360/75. 

(1) Reservoir calculations 

Based on temperatvire and pressure: 1 sec 

Based on pressure and mass flow: 30 sec 

Based on enthalpy and mass flow: 60 sec 

(2) frozen or equilibrium solution: 1/3 to 1/2 min 

(3) Nonequilibrivim soxation: 1 to 3 min 

(4) Model calculations (per model point) : 2 to 3 sec 

The times are scHnewhat greater (by up to a factor or two) when 
the planetary atmosphere models (IGAS = 5, 6) are used. Non- 
equilibrium solutions based on the helium and argon models are 
in a class by themselves; ejq>erience is limited, but times of 
che order of an hour should be anticipated. 
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3 . OUTPUTS 

Normal outputs produced in various types of NATA runs are 
described and discussed in the present section. Diagnostic mes- 
sages cind dumps produced under abnomal circumstances are listed 
and es^lained in Appendix B. 

NATA outputs are discussed in the order in which they appear 
in a normal run in Sections 3.1 to 3.8. The special outputs mak- 
ing up the transport cross section edit and the species thermal 
property edit are described in Sections 3.9 cind 3.10. 

3.1 Listings of Inpat Variables 

The output for each NATA case begins with a con 5 >lete listing 
of the input variables read in vinder namelist INPUT (Section 2.3). 
The names and values of the variables are printed in the namelist 
format, equations (1). The variables are listed in the order in 
which they are defined in Section 2.3. The printed arrangement 
of these data differs between the UNIVAC 1108 and the IBM 360. 

The 1108 prints each single variable on a separate line, and each 
array in a separate block. The 360 runs the outputs of the form 
(1) together to make up lines rtmning the full width of the page. 
The arrangement produced by the 360 is more compact but more dif- 
ficult to read. 

To provide illustrations of the various types of NATA output, 
a job consisting of two cases dealing with air flow in an axisym- 
metric nozzle was rvin on an IBM 360/75. Figure 2 is a listing 
of the input cards for these cases. Figure 3 shows the listing 
of input variables produced by the code for the first case (to 
be referred to as Test Problem No. 1) . Apart from the differ- 
ence in arrangement described above, the output shown in figure 
3 differs from that which would be produced on a UNIVAC 1108 in 
three additional respects: 

(1) Variables which are not referenced in the input 
data (figure 2) , and which are not preset in the 
coding, contain meaningless values left in the 
computer core by a previous job. Since the HOC 
clears core before each job, such variables would 
be printed as zero in the output from an 7 108 run. 
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problem no. I - a;r flow xn an axisymmetric nozzle with model and wedge 

♦INPUT 

ISW4A=l. CXMAXIsSO. ISW2B=-1. FLOW=.03. HSTAG=10880. NOZZLE=l5. 

TSD1AM=20* 25. KOIM=0. CATFAC=0. NANGLE=1. ANGLE=15. NRA0LE=1 . RADLE=.375» 
NWX=8. WXI=2.S. 4,5. lSW9B=-2 
♦END 

TEST PROBLEM NO. lA - OUTPUT CONTROLS 
♦INPUT 

ISW4A=0. iswra^l. ISW6B=-4. ISWIA=0. ISW3A=0 
♦END 




FIGURE 2 - LISTING OF INPUT CARDS FOR TEST PROBLEMS NO. 1 AND lA 
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(2) NATA is run as a double-precision program on the 
IBM 360, Hence, the exponent in floating-point 
nximbers is printed as D in place of E, and the 

alues are printed to 16 digits. 

(3) Under the IBM 360 operating system used in this 
r\in,when a double precision variable is set to 
a single-precision value either by input or in 
a data statement in the program, only the first 
six significant figures are set correctly. The 
ten trailing digits are usually inaccurate or 
meaningless. 

If data describing elements, species, or reactions are 
read in, the listing of the variables in namelist INPUT is 
followed b_' a listing of those in namelist EINPUT (Section 
2.4). Also, if transport cross section data are read in, a 
listing of the variables in namelist TINPUT is produced. These 
listings are similar in format to the INPUT listing illustrated 
in figure 3. 


3 . 2 Problem Summary 

Immediately following the listings of program inputs, NATA 
produces a summary of the input specifications for the case. This 
output is headed by the title "NATA III Code Output" , The problem 
summary for Test Problem No. 1 is shown in figure 4. In general, 
the svimmary contains the following information: 

(1) A line giving the run number (IRUN) , the case 
number in the cxirrent job (set internally by the 
code) , and the contents of the description card 
for the case, 

(2) A line stating the input specifications for the 
reservoir conditions. 

(3) A summary of the nozzle or channel geometry, in- 
cluding a table giving the parameters of the noz- 
zle profile curvefit in centimeter units. This is 
omitted in subsequent cases with the same geometry. 

(4) A specification of the gas model as either a stan- 
dard gas or a nonstandard gas. 
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(5) A table listing the cold species with their 
indices in the master list of species, mole 
fractions in the cold gas, molecular weights 
and chemical formulas. The mean molecular 
weight of the cold gas is then^ven, 

(6) A table giving the elemental con^osition of 
the gas; the "atom fraction" in this table 
is the nvurber of atoms of a given element 
per molecule of the cold gas. 

(7) A table containing the rate data for all of 
the chemical reactions in the current gas 
model. This table includes the coefficient j 
A, the temperature exponent , and the ac- 
tivation energy (in ca 1/mole) in the for- 
mula (69) of Volume I (ref. 1) for the forward 
rate constant. This table also includes the 
criterion value Cy^ for the switch from the 
perturbation solution to the nonequilibrium 
integration. Also, for each reaction with a 
third-body list, the code prints a string of 
zeros and ones containing as many characters 

as there are species in the gas model. A 1 
indicates that the corresponding species is a 
third body in the reaction, a 0 that is is not. 

(8) A table listing the chemical species in the 
gas model. In this table, a 1 under "thermal 
fit indicator" means that a thermo fit has been 
provided for the species, a 0 that no thermo 
tit has been specified. Under "alpha matrix", 
the nvimber of atoms of each chemical element in 
a molecule of the species is given. 

(9) A table giving the matrix v ! . of the stoichi- 
ometric coefficients on the product side of the 
reaction, Tlic entry under each species is the 
nmiber of molecules of the species appearing 
as a product of the reaction. 

(10) A table giving the matrix V. . of the stoichi- 
ometric coefficients on the rSactant side of 
the reaction. 
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(11) A tat)le containing data for the "physical model" 
description of the thermal properties of the 
species. The coliomns of this table contain the 
number of nuclei in a molecule of the species, 
the chemical constant b, the vibrational tem- 
perature Qy in degrees Kelvin, the enthalpy of 
formation in cal/mole, and the number of elec- 
tronic energy levels. The degeneracies g and 
energy levels E relative to the state of foma- 
tion are given in the next table, listed in pairs 
(g, E) across the page. The energies are again 

in cal/mole. ^ 

(12) A table listing the thermo fit data for the species 
for whici. such a fit is provided. The table includes 
the coefficients a, b, c, d, e, and k for each spec- 
ies, and the enthalpy of formation in cal/mole (also 
given in the table (11) for the physical model in 
the case of species for which both methods are used) . 

All of the tables from (5) through (12) are omitted 
in subsequent cases with the same standard gas model. 

(13) A statement of the species pair for which the Lewis 
number is calculated. 

(14) A statement of whether or not boundary layer effects 
are to be included in the flow solutions. 

(15) A summary of input data for calculations of condi- 
tions on models. 

(16) A summary of input data for calculations of condi- 
tions on wedge models. 

(17) A timing message giving the elapsed times since the 
beginning of the run and since the last time mes- 
sage. Such timing messages appear at several points 
in the output and allow the user to determine how 
much computer time is consumed in each major part 

of the calculations, 

(18) The value of the specific heat of the cold gas mix- 
ture at the nozzle wall temperature T^'7ALL. 




1 


1 * 


1 


3.3 Definitions of Output Identifiers 

In the first case of each NATA run, the code next prints 
a list of definitions of the alphanumeric identifiers used t 
Itibel the outputs of the flow and model condition calculations. 
This list is shown in figure 5. 


3.4 Reservoir Conditions 

Next appears a page of output summarizing the calculated 
gas conditions in the upstream reservoir, as illustrated in 
figure 6. Since the same output format is also used for the 
throat conditions, the flow velocity and mass flxox are included 
even though they are always zero in the reservoir. In Test 
Problem No. 1, the reservoir conditions were determined from mass 
flow and stagnation enthalpy inputs. The double iteration required 
in this calculation consumed over a minute of computer time, as 
shown by the timing message. 


3.5 Flow Solutions 

The nozzle flow solutions are computed and printed in the 
order: frozen, equilibrium, nonequilibrium. The format is the 

same for all three types of solution, but varies with the type 
of problem being r\an. For example, if the boundary layer on the 
nozzle wall is neglected, only a single area ratio is printed and 
the boundary layer quantities such as the displacement thickness 
and the Stanton number are omitted. In a channel flow problem 
including the boundary layer, tv;o complete sets of boundary 
layer data are printed, one for each pair of channel faces. In 
a problem involving an electron r c nonequilibrium gas model, the 
electron temperature, the radiative power loss, the energy trans- 
fer to the electrons, and the local stagnation enthalpy are added 
to the output. 

The first page of output from the frozen solution for Test 
Problem No, 1 is shown in figure 7. The output identifiers X, 

T, etc,, are defined in figure 5, The species mole fractions are 
not included in the output for the frozen solution because they 
are constant at their previously printed (figure 6) values in the 
reservoir. 


- 56 - 





FIGURE 5 - NATA CODE OUTPUT - DEFINt-’'ONS OF OUTPUT IDENTIFIERS (Firilt Page) 
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Figure 8 shows the first page of output for the equilibrixm 
solution. Here, the mole fractions are included because they 
vary along the flow. When the gas includes free electrons, the 
electron density in cm“^ is printed in place of the electron mole 
fraction. 

In the nonequilibri\am solution, the step size is often quite 
small. To avoid excessive output, the flow conditions are printed 
at intervals in the free streeim temperature controlled by the in- 
put TPRNTI (Section 2,3, Group 2). The number of integration 
steps between successive printouts of the flow conditions is stated 
in tlie line of asterisks which separates the output for different 
flow points, as illustrated in figure 9, which shows the first page 
of output frOTi the nonequi librixra solution for Test Problem No. 1. 
The method of calculation used in obtaining the conditions at the 
current flow point is also given by printing the value of the in- 
dicator, INEQ, which is 0 in the perturbation solution and 1 in 
the nonequi librium integration. 

When NATA is run with the preset value ISW6B=1, or with 
ISW6B= -1, the output includes the species mole fractions in 
every printed step as illustrated in figure 9. If |isW6b 1 is 
greater them 1, the mole fractions are printed only every I ISW6BI th 
printed step. Also, if ISV^GB is negative, the output of mole 
fractions is followed by a summary of reaction rate data as shown 
in figure 10. In this rate output, PI denotes the quantity P^ in 
equation (288) of Volume I; CHI is of equation (289); PICHI 

is Pj^ • and DLG is d In y ./dix., where y . is the concentration 

of the jth species in moles/^. The quantifies and 

are listed for all of the reactions (i = 1 to r) , together with 
the reaction index i. The queintity d In Tj /dx is listed for all 
the species together with the species names. 


3.6 Model and Wedge Conditions 

The conditions on models at a specified model point are 
printed immediately lollowing the flow conditions at the model 
point, as illustrated in figure 11. If both equilibrium and 
frozen shock calculations have been requested (FSTAG >0,), the 
model conditions for the case of the equilibrium shock are printed 
first. The species mole fractions at the stagnation point in the 
external flow over the model are then tabulated, and arc followed 
by the model conditions for the case of a frozen shock. The out- 
put identifiers for the model conditions are defined in figure 5. 
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FIGURE 11 - NATA CODE OUTPUT - MODEL AND WEDGE CONDITIONS IN TEST PROBLEM NO 
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If calculations of conditions on wedge models have been speci- 
fied in the input, the wedge results follow those for the stag- 
nation point, as shown in figure 11. The message about merging 
effects refers to equation (510a) in Volume I (ref. 1) . The 
results given under "modified Cheng-Kemp theory" are based on the 
formulas of Section 8.2.4 in Volxime I. In particular, "capital 
gamma" is T as defined by equation (501) in Volume I, and "omega" 
is as given by (502). The output quantities XW, PWW, QWW, 
etc., in the table are defined in figure 5. The wedge results 
given under "unmodified Cheng-Kemp Theory" are based on the for- 
mulas of Section 8.2.3 in Volvime I. The message referring to the 
strong-interaction approximation is based on equation (510b) of 
Volume I. 


3.7 Throat Conditions 

Immediately following the frozen solution is a page of out- 
put siamraarizing the sonic conditions in the solution. The format 
is similar to that used in the output of reservoir conditions, 
except that the transport properties are omitted. A similar 
page giving the sonic conditions in the equilibrium flow solu- 
tion follows that solution (figure 12) . Since the equilibrium 
sonic conditions are needed for starting the nonequilibrixim sol- 
ution, they are computed and printed even when the equilibrium 
solution is suppressed h setting ISW3A =0, if the nonequilibrium 
solution has been requestea (ISW2A 0) . The page containing the 
equilibrium sonic conditions also gives the coefficients o< and C 
in the analytical area-density relation (equation (383) of Vol- 
ume I) used in the upstream nonequi librium solution by the in- 
verse method. 


3.8 Informative Messages 

The normal forms of output, illustrated and discussed above, 
are sometimes interrupted with messages containing information 
of possible interest to the code user. For example, during the 
solution by the perturbation method, one line of additional out- 
put is printed for each step, giving the smallest \SX.\ (DCIIMIN) 
the largest (DCHMAX) , and the index IMAX of the largest^ 

(see Section 7.3 of Volume I (ref. 1)). Tliis type of output is 
illustrated in figure 9. 

T'Jhen the switch from the inverse method to direct nun^erical 
integration is made (Section 7.4, Volume I (ref. 1)), subroutine 
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FIGURE 12 - NATA CODE OUTPUT - EQUILIBRIUM SONIC CONDITIONS FOR TEST PROBLEM NO 
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THREAT is called to rescale the nozzle geometry’ so as to pro- 
vide continuity of effective flow area between the two solu- 
tions. THR0AT prints a two- line message as illustrated in fig- 
ure 13. The point at which this output is produced is neither 
the geometric throat nor the sonic point, but a location dis- 
tinctly downstream of both. The variables given in this output 
are as follows: 

CX Axial coordinate in nozzle (cm) 

AMACH Mach niimber 

AFNX Effective area ratio, A^, based on the inverse 

method 

DLOGA d In Ag/dx 

51 Effective area ratio, Ag, based on the nozzle 

geometry and boundary layer displacement thick- 
ness 

52 d In Ag/dx 

RSA Area rescaling factor, R^ = A^/Ag 

DELBL Nondimensional displacement thickness, ^*/Rq 

(two values in the case of a channel) 

During the nonequilibrium integration, under certain cir- 
cumstances, the code may ''freeze" a species of very low concen- 
tration which is decreasing so rapidly that it controls the stop 
size, by switching off all of the reactions in which it appears 
(Section 7.5.3 of Volvane I). When this occurs, a message such 
as the one appearing at the bottom of figure 14 is produced. 

The switch from the perturbation method to the nonequilibrium 
integration normally occurs where the largest is about 0.1. 

If one of the 152^^1 is very much larger than some of the others, 
it can cause the integration to start at a point where some of 
the reactions are very close to equilibrium. Should this occur, 
the step size required for stability of the finite difference 
equations would be too small to allow significant progress in 
the solvition. Under certain circumstances (Section 7.3.7 in Vol- 
ume I) , NATA prevents such a premature startup of the numerical 
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FIGURE 14 - NATA CODE OUTPUT - FREEZING C MINOR SPECIES MESSAGE IN TEST PROBLEM NO 
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CONC'^NTRATION UF NJC FRnZFN AT S.A'JSD-U MOLE/GM 
rsFO RCACTfON'S 17 IR 22 



> 


I 


f 


1 


1 


integration by artificially increasing the rate constant for the 
reaction with the largest for the duration of the pertur- 

bation* solution. When this occurs, the code prints a message 
such as the one illustrated in figure 15 (for a problem involv- 
ing the IGAS = 5 planetary atmosphere model) . 


3.9 Transport dross Section Edit 

Figures 16, 17 and 18 illustrate the three parts of the 
transport ciross section edit, which is invoked by setting ISWIE 
equal to -2. Figure 16 shows the edit of precoded cross ^ectien 
data for all species. In this table, STEP is the counter for the 
steps of the cross section calculation in the order in nich th 2 y 
are carried out, while INDEX is the index of the steps in array 
storage. For the preceded data, these two indices are equal for 
all steps. However, if new steps were added by input, this equal- 
ity might no longer hold. OPTI0N is the index of the option used 
to calculate the cross section generated in each step (see section 
4), The quantities W(l) to W{5) are parameter values for the 
steps. Under UnTERACTION, the pairs of species to which each 
cross section applies are listed. 

Figure 17 shows the first page of the second part of the 
cross section edit, which suinraarizes the steps or the cross sec- 
tion calculation for the current gas model. The edit shown in 
the figure is for the large standard air model (IGAS = 1). The 
format is similar to that in figure 16. 

Figure 18 shows the first page of the third part of tlie 
cross section edit, the tabulation of average! pair cross sec- 
tions as functions of temperature. The temperature is given in 


degrees Kelvin 
tabulated cross 

at 1000 degree intervals up to CTAPI. 
sections are 

The three 

1 

Q(l) 

= a (i.« 
ij 

(2a) 


Q(2) 

= a 

ij 

(2b) 


0(3) 

= a fi'i) 

(2c) 


ij ij 


(where the notation is defined in Section 3.1 of Volume I), 
and arc given in square Angstroms (vinits of 10“ 1*" . Tlie 
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FIGURE 17 - NATA CODE OUTPUT - '""ANSPORT CROSS SECTION EDIT (EDITED) 
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FIGURE 18 - NATA CODE OUTPUT - TRANSPORT CROSS SECTION EDIT (AVERAGED PAIR CROSS SECTIONS) 
AVERAGFD PAIR CROSS-SECTIONS / JNCTIONS OF TEMPERATURE 
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indices printed refer to the list of species for the current 
gas model, as given in the problem summary (Section 3.2 above). 


3.10 Species Thermal Properties 

Output from the edit of species thermal properties (invoked 
by input of a negative value for ISW6A) is illustrated by figures 
in Section 4 below. In these figures HOO denotes the species for- 
mation enthalpy, MUO the chemical potential at standard pressure 
(1 atm), H the enthalpy, CP the specific heat, and SO the entropy 
at standard pressure. The left half of the output gives the prop- 
erties as calculated from the physical model, and the right half 
gives values based on the thermo fit (see Section 2.2 of Volume I). 
The right half is left blank for species for which no thermo fit 
has been supplied. In the case of a nonstandard species for which 
no physical model data are provided, the left half of the table 
is filled with zeros. 
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4. PRECODED DATA 

To reduce the quantity of input data required for use of the 
code, NATA contains precoded data on the properties of elements 
and chemical species, on reaction rates and electronic nonequi- 
librium processes, on standard gas models, and on the geometries 
of nozzles and channels in use at the NASA/Jolmson Space Center 
Arc Tunnel Facility. This information is compiled into the 
code by means of data statements. In the present section, these 
precoded data are both summarized and, V/here applicable, docu- 
mented as to source. Section 4.1 deals with the data for ele- 
ments, Section 4.2 v;ith the thermochemical properties of species, 
Section 4.3 with reaction rates, and Section 4.4 with electronic 
nonequilibrium processes. Section 4.5 svimmarizes the six stan- 
dard gas models available in NATA. Section 4.6 reviews the trans- 
port cross section data for all of the standard species. Finally, 
Section 4.7 discusses the standard nozzle and channel geometries. 

4.1 Elements 

The only data required for the chemical elements are the 
atomi«» v;eights (in normal units of g/mole) and the chemical sym- 
bols. These data are compiled into the present version of the 
code for six elements, as follows; 

Index (IE ) Element 

1 e” 

3 He 

4 C 

5 N 


7 Ar 

Tliis list v/ill be referred to as tne "master list of elements". 

Tlie data for elements are stored in an array EPRP(I,IE), 
contained in common block /ELEM/. Tliis array is dimensioned 
(2,10), The entries in this iirray are defined as follov;s, for 
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the element with index IE in the master list: 

EPRP(I,IE) 1=1 Name of element 

1=2 Atomic weight 


For convenience in adding to or altering the compiled-in data, 
EPRP is equivalenced tc ten singly dimensioned arrays of dimen- 
sion (2), as follows: 


EPl(I) 


Equivalent to EPRP (1,1) 


EPIO(I) 


Equivalent to EPRP (I, 10) 


The data provided for elements in the current code version are 
summarized in Table I. 


TABLE I 

DATA FOR ELEMENTS 



Name 

Atomic Weierht 

E- 

5.48597 X 10"'^ 

1 

HE 

4.0026 

C 

12.0111 

N 

14.007 

0 

16.000 

AR 

39.948 



4.2 Thermochemical Data for Species 

Data for the follo\ ing chemical species £ire compiled into 
the current version of NATh: 
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Species 


1 

e“ 

2 

N 

3 

O 

4 

Ar 

5 

N2 

6 

°2 

7 

NO 

8 

no'*' 

9 

N+ 

10 

0+ 

11 

N 

2 

12 

02'^ 

13 

CO 2 

14 

CO 

15 

CN 

16 

He 1 

17 

C 

18 

C'*' 

19 

He"' 

20 

Ar+ 
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He (ground state) 




M ^ 





Index (IS) Species 


21 

He 

(^S metastable state) 

22 

He 

(^S metastable state) 

23 

Hc2+ 



24 

He2 



25 

C0+ 



26 

Ar 

(-P 2 and 

3 

Pq metastable states) 

27 

Ar 

(^PjL and 

^Pjl^ resonant states) 

23 

Ar2+ 

0 

{r molecular ion) 


This is the "master list of species" . The helium and argon spec- 
ies in specified electronic states are used in the electronic 
nonequilibrium models. 

The thermochemical data for the above species are stored 
in an array SPRP(I,IS), contained in common block /SPEC/. This 
array is dimensioned (43,30). The entries in this array are de- 
fined as follows, for the species with index IS in the master 
list of species: 


SPRP(I,IS) 1=1 Name of species 

*I = Nvunber of elements in species 

( <3) 

*I = 3—5 Indices of elements in master 

list of elements 

I = 6-8 Nxuubcrs of atoms cf elements 


*A11 values in the array are real. The values indicated by 
asterisks are converted by the program into integers. To cn- 
'varc rounding dov;n to the correct value, tlie stored values have 
been increased by 0.1. 


1 



SPRP(I,IS) 


1=9 


a 

I = 10 b 
I = 11 c 
I = 12 d 
I = 13 e 
I = 14 k 


Thermo- fit 
parameters* 


I = 15 Formation enthalpy (cal/mole) 


I = 16 


I 

I 

*** j 


17 

16 

19 


A** *** i = 20 


r , number of atoms in a mole- 
cule of the species 

b, chemical constant** 

Oy, vibrational temperature (^K) 

IGM, number of electronic lev- 
els ( <10) 

IGJ, 1 if thermo fit is used 
0 if not 


I = 21-30 gjjj, degeneracies of electronic 
levels 

I = 31-40 energies of electronic 

levels (cal/mole above ground 
state) 

X ~ 41-43 Second, thira, and fourth vib- 
rational temperatures f-''r tri- 
atomic species 


*Sce equations (33), (34) in Volume I (ref. 1). 

**Sce equation (51) in Volume I. 

***7\11 vrlucs in the array are real. The values indicated by 
asterisks arc converted l^y the program into integers. To ensure 
rounding down to the correct varuc, tlie stored values hcivc keen 
increased 0.1. 
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For convenience in adding to or altering the compiled-in data, 
SPRP is equivalenced to 30 singly dimensioned arrays of dimen- 
sion (43), as follows: 

SPl(I) Equivalent to SPRP (1,1) 


SP30(I) Equivalent to SPRP (I, 30) 


The species thermochemical data compiled into the current 
version of NATA are summarized in Tables IT through VI. The data 
are listed, in this series of tab! as, in the order in which they 
appear in the SPRP array, except that the species name is re- 
peated in each table. Table II shows how positively charged 
ionic species are represented as containing a negative nximbar of 
electrons. The thermo fits given in Teible III are from refer- 
ence 4, p. 131-132, except those for N 2 '*', 02 ”*^' CO”^, 

which are based upon unpviblished Avco calculations. The forma- 
tion enthalpy values in Table IV are based on data in the 
JANAF tables (ref. 5), NBS Circular 467 (ref. 6), the Handbook 
of Chemistry and Physics (N 2 ”^ ml O 2 ’*’) (ref, 7), Herzberg (ref. 

8) and a paper by Ginter and Brovm (ref. ?) . The electronic 
states for monatomic species in Tables V and VI are from NBS 
Circular 467. Some of the states listed are combinations of 
two or more actual states with nearly the seime energy. The 
degeneracies g for the atoms and atomic ions were calculated 
from the total angular momentvira quantum number J, which is also 
given in -_"rcular 467, using the relation 

g = 2J + 1 (j) 

The electronic energy states, degeneracies, vibrational tempera- 
tures 9^, and rotational temperatures 0^. for the diatomic spec- 
ies v;ere obtained from the spectroscopic data summarized in Table 
VII. The sources of the data are indicated in the .last colum.n 
of this table. The states of NO’*' do not appear to be very v/cll 
knowi. Rotational and vibrational constants arc given in Tabic 
VII only for the ground state of each molecule, because the phy- 
sical model in NATA assumes (as an approximation) that these 
constants arc t]ic same for all of the electronic states of each 
molecvile. The chemical constants b and vibrational tcmporiiturcs 
9^ in Table IV wore calculated from the dcita of Table VII using 
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TABLE II 


COMPOr.ITICN DATA FOR SPECIES 


IS 

Name 

No. Of 
Elements 

Indices 

IE of El 

Croats 

Nos. Of Atoms 
Of Elements 

1 

E- 

1 

1 

0 

0 

1. 

0. 

0. 

2 

N 

1 

5 

0 

0 

1. 

0. 

0. 

3 

0 

1 

6 

0 

0 

1. 

0. 

0. 

4 

AR 

1 

7 

0 

0 

1. 

0. 

0. 

5 

N2 

1 

5 

0 

0 

2. 

0. 


6 

02 

1 

6 

0 

0 

2. 

0. 


7 

N0 

2 

5 

6 

0 

1. 

1. 

0. 

3 

N0+ 

3 

1 

5 

6 

-1. 

1. 

1. 

9 

N+ 

2 

1 

5 

0 

-1. 



1C 

0+ 

2 

1 

6 

0 

-1. 

1. 

0. 

11 

N2 + 

O 

/. 

1 

5 

0 

-1. 

2. 

0. 

12 

02 + 

2 

1 

6 

0 

-1. 

2. 

0. 

13 

C02 

2 

4 

6 

0 

1. 

2. 

0. 

14 

C0 

2 

4 

6 

0 

1 . 

19 

0. 

15 

CN 

2 

4 

5 

0 

1. 

n 

0. 

16 

HE 

1 

3 

0 

0 

1. 

0. 

0. 

17 

C 

1 

4 

0 

0 

1. 

0. 

0. 

18 

C+ 

2 

1 

4 

0 

-1. 

19 

0. 

19 

HE+ 

2 

1 

3 

0 

-1. 

19 

0. 

20 

AR+ 

2 

1 

/ 

0 


19 

0. 

21 

HE3S 

1 

3 

0 

0 

1. 

0. 

0. 

22 

HEIS 

1 

O 

0 

0 

1. 

0. 

0. 

23 

HE2 + 

2 

1 

o 

E 

-1. 

2. 

0. 

24 

HE2 

1 

3 

0 

0 

2. 


0. 

25 

C0+ 

3 

1 

4 

6 

-1. 

■W 

1. 

26 

AR*M 

1 

7 

0 

0 

1. 

0- 

0. 

27 

AR*R 

1 

t 

0 

0 

1. 


0. 

28 

AR2+ 

2 

1 

7 

0 

-1. 


0. 
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TABLE III 






THERMO FIT DATA* 


0 



10^ . TFB 

10® X TFC 

10^^ X TFD 

10^® X TFE 

TFK 

3.088332 

-4.251428 

2.729295 

-0.546832 

3.071269 

4.963449 

-6.?ni753 

4.443339 

-1.000281 

5.915022 

2.083961 

-2.639548 

1.690332 

-0.361152 

3.611167 

3.749384 

-6.062030 

4.637506 

-1.107704 

4.200563 

4.47257 

-3.95880 

1.52963 

-0.21145 

4.95160 

3.37873 

-5.20841 

4.16207 

-0.97275 

4.66750 

3.22824 

-3.94364 

2.17519 

-0.42966 

4.20400 

4.33773 

-3.93346 

1.59712 

-0.23789 

5.63340 

2.10083 

-1.11714 

0.58582 

-0.13605 

4.2967 


*The thermo- fit technique is not used for monatomic 
species, for the helium and argon molecular species, 
or for carbon dioxide. 
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TABLE IV 


DATA FOR PHySICAL MODEL 


IS 

Name 

Ho° 

cal/roole 

n 

b 


IGM 

IGJ 

1 

E- 

0. 

1. 

-14.9276 

0. 


0 

2 

N 

112520. 

1. 

0.2944 

0. 


0 

3 

0 

58989. 

1. 

0.4938 

0. 


0 

4 

AR 

0. 

1. 

1.8664 

0. 


0 

5 

N2 

0. 

2. 

- 0.4106 

3352. 

5 


6 

02 

0. 

2. 

0.1140 

2239. 

5 


7 

N0 

21456. 

2. 

0.5455 

2699. 



8 

N0+ 

236660. ^ 

2. 

0.3841 

3373. 



9 

N+ 

447600. 

1. 

0.2943 

0. 


^9 

10 

0t 

372940. 

1. 

0.4938 

0. 


^9 

11 

N2+ 

357680. 

2. 

- 0.3763 

3114. 



12 

02+ 

288800. 

2. 

- 0.0317 

2628. 

5 


13 

C02 

0.* 

3. 

1.8958 

1977.** 

1 

0 

14 

C0 

66770.* 

2. 

0.3169 

3083. 

5 


15 

CN 

197170.* 

2. 

0.2226 

2939. 

3 


16 

HE 

0. 

1. 

- 1.5846 

0. 

mm 

^^1 

17 

C 

263550.* 

1. 

0.0637 

0. 


0 

18 

c+ 

523310.* 

1. 

0.0636 

0. 

19 

0 

19 

HE+ 

566840. 

1. 

- 1.5846 

0. 

■9 

0 

20 

AR+ 

363318. 

1. 

1.8663 

0. 

2 

0 

21 

HE3S 

456910. 

1. 

- 1.5846 

0. 


0 

22 

HEIS 

475260. 

1. 

- 1.5846 

0. 


0 

23 

HE2+ 

511490. 

1. 

- 3.5619 

2343. 


0 

24 

HE2 

406170. 

1. 

- 3.6287 

2492. 

99 

0 

25 

C0+ 

389950.* 

2. 

0.2931 

3142. 

3 

1 

26 

AR*M 

266350. 


1.8663 

0. 


0 

27 

AR*R 

267970. 

19 

1.8663 

0. 

^9 

0 

28 

AR2+ 

337040. 

2. 

3.597 

115. 


0 


*These formation enthalpies for the carbon-containing species are 
based on CO2 as the reference state for carbon. They are 93970 
cal/mole higher than the JANAF values, which are based on graphite 
as the reference state. 


**The other three vibrational temperatures for CO2 are 960, 960, 
and 3380OK. 
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ENERGIES OF ELECTRONIC STATES 
(CAL/MOLE) 
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equations (38), (40b), (51b), and (45) of Volume I, i.e,, 

b = -3.66505 + 2. lnW-(n-l)ln ll. 43879 - % «e>] (4) 

= 1.43879 (ft)e - 2 ‘♦^e ^e^ <5) 

in which W denotes the molecular weight and n the number of atoms 
per molecule (for monatomic and diatomic species). The degen- 
eracies g for the molecular states in Table V were obtained from 
the state symbols using the following rule, based on Herzberg 
(ref. 8): 


g = n 

( 6 ) 

g = 2n 


S states 
^TT , , etc. 


Figures 19 through 46* are tables of thermal properties for the 
corapiled-in species, computed by NATA from the precoded data, 
under the option ISWGA = - 1. For each species, the properties 
calculated from the physical model are given for temperatures 
up to 30,000°K. For the molecular species for which the thermo 
fit is used, corresponding results based on the thermo fit are 
given for comparison. At the normal switchover temperature, 

CTMXX = 5000°K, the results from the two techniques are seen to 
be in reasonably good agreement. Above about 15,000°K, the prop- 
erties calculated from the thermo fits become very inaccurate in 
all cases. This behavior results from the inability of the poly- 
nomial form used in the thermo fit to represent the actual prop- 
erty variations over excessively wide temperature ranges. This 
is not considered to be a serious problem because, at tempera- 
tures of 15000°K and higher, the mole fractions of the molecular 
species (for which the thermo fit is used) are quite small, so 
that errors in the thermal properties of these species have only 
very minor effects upon the properties of the gas mixture. How- 
ever, at temperatures above about 20,000®K, the properties com- 
puted from the thermo fit show wild variations which could lead 


*The equipment used to produce these figures printed the BCD 
sign appearing in some of the species names as an ampersand. 
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FIGURE 43 

THERMAL PROPERTIES OP COC ( HOO* 380*950 KCAL/MOLE > 
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to code failure or grossly inaccurate results. Since the thermo 
fit is not used in the helium and argon models, these models can 
be used at higher temperatures, up to the onset of significant 
second ionization. 


4,3 Data for Reactions 


Tables VIII, IX, X, and XI list the chemical reactions for 
which data have been compiled into the current version of NATA. 
'Hiese tcibles define the "master list of reactions". They also 
give the parameters A, , and in the curvefit* 


kf 


A ( 


io,ooook 


) 


^-Ea/RoT 


(7) 


to the forward rate constant for each reaction. The reaction 
system for air (Table VIII) is one recommended recently by Cornell 
Aeronautical Laboratory (ref. 14) , Those for argon and helivun 
(Tables IX, X) are documented in Appendix A of the present report. 
The reactions for the carbon-containing species (Table XI) are 
from Dvinn (ref, 15) , except the CN reactions which are from McKenzie 
and Arnold (ref, 16) . 

The rate constants of interest in NATA applications are not 
accurately known. Experimental determinations of the rate for a 
given reaction at a given temperature typically differ by factors 
ranging from 2 up to an order of magnitude or more. In many cases, 
the temperature range over which the rate has been measured is 
considerably smaller than the range over which it is used in NATA. 
The extrapolation which is thus required is a further source of 
error. Hov/ever, in spite of these vincer taint ies, calculations 
based on the reaction system for air (Table VIII) have given results 
in reasonable agreement with experimental data (ref, 14) , 

The main sources for the dfita in Table VIII are cited in reference 
14. The reactions for the neutral species in high temperature air 
have been reviewed by Wray (ref. 17) , those for the charged spec- 
ies by Dunn and Treanor (rof , 18) , 


♦Equation (69) of Vo lime I. 
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In Table VIII, the reactions given in reference 14 hc've been 
rearranged to place the reactions having the same effect (e.g., 
dissociation of O 2 ) together, and to placa the reactions involv- 
ing only neutral species ahead of those involvinr charged species. 
Reference 14 gives both forward and badcward rates for each reac- 
tion. In each case, the rate coefficient in one direction is 
based directly on chemical kinetic data, while that in the ether direc- 
tion is a curvefit to the recults of calculations based on equation (62) 
of Volvune I. In Table VIII, the reactions have been written in such a 
form that t?.e rate based on experimental data is in th3 forward direc- 
tion. In NATA, the baclcward rate is computed internally using equation 
(62) of Volume. I. , 

The reaction rate data summarized in Tables VIII to XI are 
stored in an array RPRP(I,IR), contaii'.ed in common block /REAC/ , 

This array is dimensioned (29,92). The entries in this array 
are defined as follows, for the reaction with index IR in the 
master list of reactions: 


RPRP(I,IR) 1=1 

1 = 2 - 
1 = 3 

1 = 4 

*1 = 5 

* 1=6 
*I = 7-9 

*I = 10-12 


Coefficient A in eq, (7) 

(cm3 mole sec or cm^/mole2 sec) 

Exponent "Vj in eq, (7) 

Activation energy E^ in eq. (7), 
caV mol's 

QQ: 1,0 if a third-body list is 

provided in I = 20-29 
0,0 if not 

Number of reactant species { < 3 ) 

Number of product species f <, 3) 

Indices of reactant species in 
master list of species 

Indices of product species in 
master list of species 


*All values in the array are real. The values indicated by aster- 
isks are converted by the program into integers. To ensure round- 
ing down to the correct value, the stored values have been increased 
by 0,1 
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I 


I = 13-15 


Nvuiibers of molecules of reac- 
tants 


I = 16-18 


Numbers of molecules of pro- 
ducts 


*1 = 19 


Nuniber of third bodies ( <10) 


*I = 20-29 


Indices of third body species 
in mastv'ir list of species 


For convenience in adding to or altering the compiled- in data, 
RPRP is equivalenced to 92 singly dimensioned arrays of dimension 
(29), as follows: 

RPl(I) Equivalent to RPRP (1,1) 


• • 

RP92(I) Equivalent to RPRP (I, 92) 

4.4 Electronic Nonequilibrium Data 

When an ionized gas expands to low density, as in the di- 
verging nozzle section of an arc-heated wind tunnel, a condition 
of nonequilibrium betv/een the electron teiiperature and the heavy- 
particle temperature develops. Tliis phencatienon is a result of 
the smallness of the cross section for elastic energy transfer 
betv^een electrons and heavy particles. At high gas densities, 
there are enough collisions to keep the two temperatures approx- 
imately in equilibriixm in spite of the small cross section; but 
at low densities this is no longer true. In an e^qjanding plasma 
flow, electron-ion recoatrihination processes supply energy to the 
electron gas, so that the electron temperature normally falls 
more slowly than the gas temperature. 

The forward rate for a reaction which includes the electron, 
as either a reactant or a third body, usually depends upon the 
electron temperature rati er than the gas temperature. Thus, 
tliermal nonequilibrimi bet.\/een tlie electrons and the heavy par- 
ticles can affect the rates of production and destruction of 
species. Tliis phenomenon is not considered in the current NATA 
models for air. It has been studied, for the case of an atomic 
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nitrogen plasma, by Bowen and Park (ref, 19) . Nonequilibrium be- 
tween the electron cind gas temperatures is included e:q)licitly in 
I the NATA models for helixim and argon (Volume I, Sec. 7.1.2). 

These noble gas models require additional data for each reaction 
to specify the dependence of the forward and backward rate con- 
stants upon the electron and gas temperatures and to specify the 
partition of the energy of reaction between electron kinetic cn- 
ergy and radiative losses. Data are also required to provide the 
elastic collision cross section between the electrons and the neu- 
tral heavy particles (assumed to be the same for all neutral spec- 
ies) . These extra data for the thermal nonequilibrium models are 
stored in an array TNEP(I,INT), contained in common block /TNE/, 

This array is dimensioned (186,2). The index INT specifies the 
gas model in which the data are to be used. For helium, INT = 1, 
and for argon, INT = 2. The entries in TNEP(I,INT) are defined 
below in terms of the reaction index JR for whichever gas model 
is being used. The index JR runs from 1 up to ISR. The uxraen- 
sioning of TNEP allows the number of reactions ISR to be as high 
as 25 for gas models involving electronic nonequilibrixjm. For 
conventional gas models, ISR can be as high as 64. 

The entries in TNEP (I, INT) are defined as follows; 

*I s= JR KTF (JR) , indicator of type of formula for 

I forward rate constant for the JRth reaction: 

KTF =1 kf given by (7) as a fvinction 

of gas temperature T 

KTF =2 k£ given by (7) as a function 

of electron temperature Tg 

KTF = 3 kf = A (Te/lO"^)^ (1 - 

KTF = 4 kf = A (Te/10^) /max (1, ^ ) 

where T = b np R/Nq 

KTF = 5 kf = A/ /r 

Note; R denotes the local nozzle radius (or a corresponding ef- 
fective value in the case of a channel) . Also, n is the nvimber 
density of the atomic species appearing on the product side of 
the reaction, and b is a coefficient stored as BPAR (below) . The 
types KTF = 3, 4, and 5 are used only in the model for argon 
(IGAS - 3) ; see Appendix A. 

*All values in the array are real. The values indicated by aster- 
isks are converted by the program into integers. To ensure round- 
ing down to the correct value, the stored values have been increased 
by 0.1, 
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*I = JR + 25 


I 


*I = JR + 50 


KTR(JR) e equal to 0 if t. backward rate con- 
stant kj. = 0 for the JRth reaction; equal to 1 if 
kj. = kj-(T); equal to 2 if kj. = kj.(Te), where the 

fijnctional dependence xs given by equations (277-278) 
of Vol , I (ref. 1) . 

ITR(JR), indicator of rule for partitioning the 
reaction energy for the JRth reaction. In the 
definitions below, and - €j. denote the en- 
ergies gained by the electron gas in Nq reactions 
in the forward and reverse directions, respectively, 
and qf»“ qr <3enote the corresponding energies lost 
by radiation. Also Nq =Avogadro*s number. 


ITR =1 = - a Rq 

^ ^0 “ ^f 


c = q =0 


ITR = 2 


'•f 

‘If 


3 


- f ""O ""e 


^0 " ^f 


6 = 


= q^ = 0 


ITR = 3 = qf = " ^r " ° 


ITR = 4 = “ I ^0 "^e 

qf = qr = 0 


ITR = F 


= ^0 

q£ = q^ = 0 


ITR = 6 


qf - 


<^f = ^r = = 0 


*A11 values in the array are real. Tlie values indicated by aster- 
isks are converted by the program into integers. To ensure round- 
ing do\%’n to the correct value, the stored values have been increased 
by 0.1. 
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I 


I = JR + 75 

I = JR + 100 

I = 126-155 

I = 156-185 
I = 186 


EPAR(1,JR), parameter for the JRth reaction 
tal per Nq reactions) 

EPAR(2,JR), parameter "a" for the JRth reaction 
if ITR(IR) = 1 

TLIST(J), temperature values for table of elastic 
collision cross section, (see Appendix A) 

P0M(J) , values for table 

BPAR, parameter b for all reactions with KTF = 4 


For convenience in adding to or altering the compiled-in data for 
reactions involving electronic nonequilibrium, "mEP is equivalenced 
to two singly dimensioned arrays of dimension (186), as follows: 

TNl(K) Equivalent to TNEP(K,1) 

TN2(K) Equivalent to TNEP(K, 2) 

Tables XII through XV stammarize the precoded electronic non- 
equilibrium data for the helixara and argon models (IGAS = 4 and 3) . 
The data in these tables are docvuuented in Appendix A. 


TABLE XII 

ELECTRONIC NONEQUILIBRIUM DATA FOR HELIUM MODEL 


JR 

IR 

KTF (JR) 

KTR(JR) 

ITR(JR) 

EPAR(1,JR) 
ca 1/mole 

1 

35 

2 

2 

5 

109890 

2 

36 

2 

2 

5 

91540 

3 

37 

2 

0 

2 

109890 

4 

38 

2 

0 

2 

91540 

5 

39 

1 

1 

3 

0 

6 

40 

2 

2 

5 

98040 

7 

41 

2 

2 

5 

458270 

8 

42 

2 

1 

4 

0 

9 

43 

2 

1 

4 

0 

10 

44 

n 

/. 

2 

5 

18350 

11 

45 

2 

2 

5 

456730 

12 

46 

2 

2 

5 

475080 

13 

47 

1 

0 

6 

475080 

14 

48 

1 

2 

5 

346840 

15 

49 

1 

2 

5 

365190 

16 

50 

1 

2 

5 

383540 

17 

51 

1 

1 

3 

0 

18 

52 

1 

2 

5 

260350 

19 

53 

2 

2 

5 

413480 
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TABIE XIII 


ELECTRONIC NONEQUILIBRIUM DATA FOR ARGON MODEL 



KTF(JR) KTR(JR) ITR(JR) EPAR(1,JR) EPAR(2,JR) BPAR 

cal/niole 


96970 

95360 

96970 

95360 

363530 

266350 

267970 

267970 

1600 

0 

226000 

0 

226000 

0 

70680 

69070 

337040 


























TABLE XIV 

e" - He MOMENTUM TRANSFER CROSS SECTION 


Te 

°K 

5(1,1) 

h 

°K 

5(1,1) 

8" 

0 

0.00 

9000 

6.80 

100 

5.00 

10000 

6.77 

200 

5.59 

12000 

6.37 

400 

5.83 

14000 

6.55 

600 

5.99 

16000 

6.42 

800 

6.11 

18000 

6.29 

1000 

6.21 

20000 

6.15 

1500 

6.39 

25000 

5.80 

2000 

6.51 

30000 

5.46 

3000 

6.67 

35000 

5.14 

4000 

6.76 

40000 

4.85 

5000 

6.81 

45000 

4.60 

6000 

6.84 

50000 

4.32 

7000 

6.84 

70000 

3.4 

8000 

6.83 

100000 

2.4 
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TABLE XV 


e"*-Ar MOMENTUM TRANSFER CROSS SECTION 


Te 

°K 

82 

Te 

®K 

^(1,1) 

^°2 

PT 

0 

10. 

2500 

0.51 

25 

10. 

4000 

0.98 

50 

6. 

6000 

1.69 

100 

3.57 

8000 

2.40 

200 

2.12 

10000 

3.08 

300 

1.39 

12000 

3.73 

400 

0.97 

15000 

4.65 

600 

0.57 

20000 

6.02 

800 

0.40 

25000 

7.19 

900 

0.35 

30000 

8.11 

1000 

0.32 

35000 

8.81 

1200 

0.29 

40000 

9.30 

1400 

0.29 

45000 

9.63 

1600 

0.31 

50000 

9.82 

2000 

0.38 

100000 

10. 
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4.5 Standard Gas Models 

A "gas model" is the specification of a set of species with 
their thermochemical properties, a system of reactions among these 
species with their rate constants, and other data. The NATA code 
contains provisions for invoking certain standard gas models by 
input of a single index value, IGAS, The standard gas models 
available are summarized in Table XVI. The third and fourth col- 
umns in this table specify the pair of species whose binary dif- 
fusion coefficient is used in computing the Lewis number. The 
variable INT pertains to the treatment of electronic nonequili- 
briun in the model. For INT = 0, electronic nonequi librium is 
neglected. For INT > 0, electronic nonequilibrium is taken into 
account, and INT is the index for selecting the reaction para- 
meters from TNEP(I,INT). If the indicator LEWIS is equal to 1, 
the Fay-Riddell Lewis nvunber factor is used in calculating the 
stagnation point heat flux. For LEWIS =2, it is not. The reac- 
tion indices in Table XVI refer to the master list of reactions. 

AIR-1 is the general model ^or argon-free air, suitable for 
use in cases with reservoir temperatures up to about 15,000- 
20,000°K, Temperatures above 15000°K are beyond the range of 
validity of the thermo fits for the diatomic molecules. Above 
20,000°K, the specific heats for some of these species (as com- 
puted from the thermo fit) go negative, and the chemical poten- 
tials begin to decrease with increasing temperature. However, 
the tenperature capability of the AIR-1 model appears more than 
adequate for describing the flow in state-of-the-art arc heated 
wind tunnels. 

AIR-2 is a truncated air model obtained from AIR-1 by delet- 
ing all of the ion species except NO’*’. It is suitable for use in 
cases with reservoir temperatures up to about 6000®K. The ion 
NO"*" is retained because the ionization potential of NO p.5 ev) 
is much lower than those of the other neutral species in high temper- 
ature air (12.5-15.5 ev) . In problems with reservoir temperatures be- 
low 6000®K, use of the AIR-2 model in place of AIR-1 economizes on 
computer time without significantly affecting the results. 

HELIUM and ARGON are the electronic nonequilibrium models 
for helium and argon. Since the thermo fit is not used in these 
models, they are suitable for use up to temperatures at which 
the doubly ionized species He'*"*' and Ar’*^'*' become important. 
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CONAR* is a model for a planetary atmosphere containing 75 
mole percent CO 2 , 20 mole percent argon, and 5 mole percent N 2 . 
These mole fractions can be adjusted easily in the code input, 
so that CONAR is usable as a general C02~Ar-N2 model for the at- 
mospheres of Venus and Mars. 

C0NAR2 is a smaller version of CONAR with some of the ion 
species omitted, designed for use at temperatures up to 7000°K. 

In this temperature range, it gives practically the same results 
as CONAR with less expenditure of computer time. 

The data required for generating these standard gas models 
are stored in cin array GPRP (I, IGAS) , which is contained in common 
block /MIXT/. This array is dimensioned (124,6). Its entries are 
defined as follows, for the model with index IGAS: 


I = 
**I = 
**I = 
**I = 
**I = 

**I s= 

I = 
**I = 
**I = 
**I = 

**1 = 


1 Mixtvire name 

2 Number of elements in mixture (ISC) 

3 Nviinber of species in mixture (ISS) 

4 Number of reactions included (ISR) 

5 Niariber of ion species (IC) 

6-15 Indices (IE) of elements in master list of elements 

16-25 Mole fractions of cold species (QPJ) 

26-45 Indices (IS) of species in master list of species 

46-109 Indices (IR) of reactions in master list of reactions 

110-119 Indices (JCS) of cold species in master list of 
species 

120 Nvimber of cold species (NCS) 


*Acronym for Carbon-Oxygen-Nitrogen-Argon. 

**A11 values in the array are real. The values indicated by aster- 
isks are converted by the program into integers. To ensure round- 
ing do\7n to the correct valvie, the stored values have been increased 
by 0.1. 
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Atom index (ISAT0M) for Lewis number in master list 
of species 

Molecule index (ISI^L) for Lewis number in master 
list of species 

INT. If INT = 0, electron temperature equals gas 
temperature. If INT > 0, the model includes elec- 
tronic nonequilibrivim, and INT is the index of 
extra reaction properties required in TNEP(i,int) 

Indicator (LEWIS) for inclusion (1) or exclusion 
(2) of ?ay-Riddell Lewis number factor in stagna- 
tion point heat flux. 

Most of these definitions are self explanatory. However, the 
reference to "cold species" requires some discussion. In NATA, 
the overall composition of the gas in terms of the chemical ele- 
ments is specified by giving the composition of the cold gas mix- 
ture which is fed into the arc heater. The chemical species in 
this cold gas are called cold species. For example, the cold 
species in argon- free air are N 2 and O 2 , and their mole fractions 
are assumed to be 0.78823 and 0.21177, respectively. The weight 
fractions of the elements nitrogen and oxygen, which are deter- 
mined by these data, are invariant under all chemical changes in 
the system. 

For convenience in adding to or altering the compiled-in gas 
models, GPRP is equivalenced to 6 singly dimensioned arrays of 
dimension (124) as follows: 

GPI(I) Equivalent to GPRP(I,1) 


GP6(I) Equivaljcnt to GPRP (I, 6) 


4.6 Transport Cross Section Data 

The cross section data required for calculating the trans- 
port properties of arbitrary mixtures of the standard species 
(listed in Section 4.2) are compiled into NATA. The methods 
used in the transport property calculations have been explained 
in Section 3 of Volume 1 (ref. 1) . Briefly, the transport coef- 
ficients are computed from formulas involving the cross sections 


*I = 121 
*I = 122 

*I = 123 
*I = 124 
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pairs^if species ti/j). 


and B*. n^^'. ^^for collisions between 

These^cross^iections are calculated in 
a series of steps. First, the cross sections for _all pairs _are 
set to zero. Then, in each step, the values of X).t2,2)^ 

and B* are computed by a particular method (or “option") 

with a particular set of parameter values, and these values are 
added to the corresponding cross sections for each pair of spec- 
ies to which the step is applicable. The information concerning 
the applicability of steps to species pairs is stored in index 
arrays KKQ(M), NNQ(M), as explained below. If only one step of 
the cross section calculation is applicable to a particular species 
pair, then the cross sections f<^r the pair are the values computed 
during that step. If two or more steps are applicable to the pair, 
the cross sections for the pair are built up by adding contribu- 
tions from the different steps. If the cross sections are poorly 
known for several minor pairs of species, but ar«_ considered likely 
to be roughly the same for all pairs, then the cross sections for 
all of these pairs can be set in a single step. 


In the present section, the twelve options for calculating 
cross sections are defined, and the default methods used by the 
code to determine xanspecified cross sections are explained. The 
variables and arrays used to store the precoded transport cross 
section data are then defined. Finally, the precoded data for the 
standard species are tabulated and documented as to source. 


The options for calculating cross sections are selected by an 
index KKQ. For each option, tl re is an associated list of input 
parameters, W(J). Array dimensions limit the number of these 
parameters to five. Other nxomerical data required by some of the 
options are stored at specified locations in four arrays (TL, 
0MEGA1, ASTAR, BSTAR) , as explained below. Each of these arrays 
is dimensioned (1000). The cross section options available in 
NATA are as follows; 


KKQ = 2 Coulomb Cross Sections 


Here 

= 0.8 W(l) 

= W(2) (8) 

B* = 1,5625 
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where Q_ is defined by 


Qc = In (7 A ) 

(9a) 

A = | 

2 { (T n^) 

(9b) 

y =\l + ng^/^ ] ^ 

L 9 kT ^ 

(9c) 

(Section 3.2 of Volume I). In equations (9), e denotes the elec- 
tron charge, k Boltzmann's constant, T the absolute temperature 
and ng the electron density. 

KKQ = 3 Exponential Potential 


In this option, the cross sections are 
(ref, 20) tabulated collision integrals for 
tial 

obtained from Monchick 
the exponential poten- 

0 = Ae-^/P 

(10) 


which are compiled into NATA in the TL, 0MEGA1, ASTAR, and BSTAR 
arrays starting at location 1 in each array. The parameters are 

W(l) = A/k in°K 

W(2) = p n A (11) 

W{3) = 1,0 = position of first entry in tabulated 
collision integrals 


KKO = 4 Charge Exchange Gross Section 

In this option, -il and B* are calculated for a resonant 

charge exchange cross section of the form 


Q = (A - B log,^ V 
ex ' 10 


-13 9-- 



where v is the relative velocity in cm/sec. fi is not 

calculated in this option. Tlie required input parameters are 


W(l) = A in A 
W(2) = B in A 

W(3) = molecular weight of atom 


( 12 ) 


W(4) = control parameter 


For W(4) > 0., the computed cross sections il ^ ^ md B* 
replace those computed in earlier steps of the calculations, 
whd.le for W(4) ^ 0., they are adaed to the earlier values. 


KKO = 5 Tabulated Cross Section 

In this option, the cross section data are gxven in tabular 
form as a function oi temperature. The input parameters are 

W(l) = A = factor by which the tabulated values must 
be multiplied to give the collision integrals 
in 

W(2) = I - position of first entry in tabulated cross 

section data 

W(3) = N = mmber of entries in cross section table 


The cross section data themselves are stored in the TL, OMEGAl, 
ASTAR, and BSTAR arrays, starting at element I, as follows: 


TL(I) to TL(I - 1 + N) 


)Z5MEGA1(I to I - 1 + N) 
ASTAR(I to I - 1 + N) 


temperatures at which cross 
section data arc tabulated in 
^K. Values must be in order of 
increasing temperature, 

values of (^»^) at the tab- 
ulated temperatures 

values of -Ti- ' at the t£ib- 
ulatc^d temperatures 
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BSTAR(I to I - 1 + N) 


values of L* at the tabulated 
ten 5 >eratures 


KKP = 6 Power Law Potential 

This option calculates cross sections for an inverse power 
potential, 

^ = Ar”^ (13) 


based on the analysis of Kihara, Taylor, and Hirschfelder (ref. 
21) . The parameters are 

W(l) = ITL = index in 0^^EGA1, ASTAR, BSTAR arrays 
where data are stored 

w(2) = 'n 


For each value of used, the following additional data are 
stored; 


0MEGA1(ITL) = 
ASTAR (ITL) = 


TTr(3-4) (14a) 

(|-i) (v) )/A<^> (») ) (14b) 


BSTAR (ITL) 


(1 




(14c) 


where A/k is in ®K, A^ (>^ ) and A^ ^ (i^ ) are tabulated functions 
which are given for both attractive and repulsive potentials in 
reference 21, and T denotes the gamma function. 


KKQ = 8 Lennard-Jones (G-12) Potential 

This option calculates cross sections for the Lennard-Jones 
(6-12) potential, 

0(r)=4& [(^-)^^ - (.^)S] (15) 
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The parameters are 

W(l) = €/k in ^ 

W(2) = (Tin X 

W(3) = 501. 

Tabulated collision integrals for the Lennard-Jones potential are 
compiled into the code in the TL, 0MEGA1, ASTAR, and BSTAR arrays, 
starting at location 501. 


KKQ = 9 Scaling of Previously Commuted Cross Sections for Other 
Species 


This option allows cross sections calculated for one pair 
of species to be used also for other species, possibly with a 
constant multiplying factor. The cross sections are calculated 
from the formulas 




[L 


^( 2 . 2 ) 


= h 

12 ij 


B* = C B*. 

3 3-3 


where the Cy, are constant factors and the subscript ij indicates 
cross sections calculated previously for the pair (i,j)-. The 
parameters for the option are 

W(l) = i ■= first index of previously calculated cross 
section 

W(2) = j = second index of previously calculated cross 

section 

W(3) = 

W(4) = C 2 

W(5) = C 3 
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KKQ = 10 Empirical Mixing Rule 

This option calculates the cross sections for a pair of un- 
like species i,j (i ^ j) from the empirical mixing rule 



The values calculated from (17) are then added to the previously 
calculated cross sections for the pair. This option uses no W 
parameters. 


KKQ = 11 Fairing Option 

This option modifies the previously calculated cross value 
for a species pair according to the formula 


new old 

where f(T) is a linear fairing factor given by 

f(T) = max [o, min (1, )] 

^1 ^0 


(18a) 


(18b) 


Use of this option thus permits different forms to be used for 
the cross section in different parts of the temperature range, 
with a smooth transition between them. The parameters are 

W(l) = Tq = temperature at which the -Tl are to be set 
to zero 

W(2) = *^1 “ temperature at which the -f)- are to remain 

unchanged 


KKQ = 12 Generalized Mixing Rule 

This option is a generalization of the empirical mixing 
rule KKQ = 10 in which the cross sections are calculated from 
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the formula 




/ {X (^»s) 

^yxi(i,s)j 

(19) 



[y ij 

' mn 

where i, j, and 

m 

, n are any specified molecular pairs. 

The 

parameters are 






W(l) = 

i 



W(2) = 

j 



W(3) = 

m 



W(4) = 

n 



KKQ = 13 Scaling of Previously Computed Cross Section for 

the Same Species Pair 

This ^tion calculates one of the averaged collision cross 
sections for a ^air of species from previously calculated 

values of a different for the pair. In terms of the no- 

tation 

a 

a - Ai (20) 

ij ij 

a 


(1) 


(1,1) 

ij 

ij 

(2) 

^ n 

(2.2) 

ij 


ij 

(3) 

= B* 

a 


ID 


ID 


the option calculates a new value of the cross section n'”’. 
from the formula ^3 


= c 

13 ij 


( 21 ) 


where m and n are two specified integers in the range 1 ^ m 3, 
1 < n < 3 and C is a constant. The newly calculated value of 
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then IE places 
The parameters are 


the previous value of this cross section. 


w(l) 

= m 

W(2) 

= n 

W(3) 

= c 


KKQ = 14 Hultiplication bv a Constant 


This option multiplies previously calculated values of the 
collision cros'- sections for a pair of sp>ecies by a constant 
factor, accord.’ng to the formulas 


^( 1 . 1 ) , 

ij 


c, XI 

^ ij (old) 

-Q- 

ij (old) 


B*. 

^3 


^3®ij (old) 


( 22 ) 


The option is the same as KKQ = 9, except that here the cross 
sections for a species pair are obtained from previously calcu- 
lated values for the same pair, instead of from values for a dif- 
ferent pair as in KKQ =9. Parameters for the option are 

W(l) = 

W(2) = 

W(3) = 


NATA contains default provisions for estimating some cross 
sections if they are not specified explicitly in the preceded 
data or the input. If none of the specified steps in the cross 
section calcvilation is applicable to a particular pair, and if 
both of the species are ions, then the effective Coulomb cross 
sections (8) are used. If one species is neutreil and the otic r 
ionized, the formulas 
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( 23 ) 


^( 1 , 1 ) ^ ^( 1 ) ^- 0.4 

^( 2 , 2 ) ^^ 0.4 

B* T"°-^ 


(m^ 

give the default option. The constants a' * are canpiled into 
the program in the locations 0MEGA1 (996), ASTAR(996), and 
BSTAR(996) for m = 1, 2, 3, respectively. If both species are 
neutral and unlike (not the same species ) , the cross sections 
are estimated using the mixing rule (17). However, if cross 
section data are not specified for like- like collisions of a 
neutral species, the code does not attempt to provide estimates 
of the cross sections, but returns an error message and termi- 
nates the ca.'e. 


The variables and arrays used to store the precoded cross 
section data are as follows: 


Var iab le 

Name Dimension 


Definition 


NNKQ 1 

NNQ (M) 100 

NNQ(M) 100 

Ilm(K) 5 

JJm(K) 5 


Number of steps in the cross section cal- 
culation for which data are specified. 

Index specifying the option to be used 
in the Mth step of the cross section cal- 
culation (see above) . 

Number of species pairs to which the cross 
sections calculated in the Mth step are to 
be applied (NNQ(M) < 5), 

Indices of the species to which the cross 
sections calculated in the mth step are 
to be applied, referred to the master list 
of species (Section 4.2). In these vari- 
able names, m denotes an integer which is 
part of each name. Thus, for example, 
1123(2) and JJ23(2) are the indices de- 
fining the second pair of species to which 
the cross sections calculated in the 23rd 
step arc applied. There are 100 arrays 
of each type, c.g., Ili(K), II2(K),..., 
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IIIOO(K). There are NNQ(m) pairs of 
indices set for each step m. Only pairs 
with Ilm(K) < JJm(K) are used in the cal- 
culations.* 

Wm(K) 5 Parameter values for the mth step of the 

cross section calculation (see discussion 
of KKQ options above) . There are 100 of 
these arrays, Wl(K), WlOO(K).* 

ISEQ(L) 100 Sequencing array for specifying the order 

in which the defined steps are carried 
out during the cross section calculation. 
The index M or m in the preceding arrays 
is given by M ~ ISEQ(L), where L = 1, 2, 
3, ..., NNKQ. 


TL 

1000 

Additional 

array storage for cross sec- 

0MEGA1 

1000 

tion data. 

The data compiled into these 

ASTAR 

1000 

arrays are 

discussed below. 

BSTAR 

1000 




To prevent the data statements used in setting the TL, 0MEGA1, 
ASTAR and BSTAR arrays frc*n exceeding the 20-card limit in For- 
tran IV, those arrays are equivalenced to 40 ai^avs each of dim- 
ension (100), as follows: 

TLl(l) equivalent to TL(1) 

TL2(1) equivalent to TL(lOl) 


TLlO(l) 

TLll(l) 


equivalent to TL(901) 
equivalent to 0HEGA1(1) 


TL20(1) 


equivalent 


to 0MEGA1(9O1) 


*The arrays Ilm(K) , JJm(K) , Wn:(K) arc the same as the input ar- 
rays Im(K), Jm(K), Vm(K) discussed in Section 2.4. The shorter 
neiracs are used for input to keep the defining statement for name- 
list TINPUT within the 20-card limit allowed by Fortran IV. 
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TL21(1) equivalent to ASTAR(l) 


TL3Q(1) 

TL31(1) 


eqtiivalent to ASTAR(90l) 
equivalent to BSTAR(l) 


TL40(1) equivalent to BSTAR(901) 


The precoded data for cross section calculations will now 
be tabulated and docvimented . These data have been taken from 
previous transport property studies at Avco Systems Division and 
have not been revised during the present program. Thus, in scane 
cases the values used in NATA may not represent the latest avail- 
able data. Nevertheless, the data in the cods should be generally 
satisfactory for most engineering applications. 

For most of the important cross sections, with the excep- 
tions of those involving carbon-containing species, the precoded 
cross section data should be accurate to within 20 to 40 percent. 
For the carbon-containing species, very few data are currently 
available on the collision cross sections at high temperatures, 
and the values used in NATA are based for the most part on rough 
estimates. In general, it is believed that these estimates should 
be accurate to within about a factor of two. For some interactions 
involving minor species such as the raetastable states of He and 
Ar, nominal cross sections are used which may be in error by large 
factors. However, because of the low concentrations of the spec- 
ies in question, the effects of these cross section errors upon 
the calculated gas transport properties are small. 

The sources of the cross section data used in NATA arc indi- 
cated, for each pair of species in the master list, in Table XVII, 
In this table, the numbers preceded by A in the third column re- 
fer to the notes at the end of the table, while the numbers in 
the final colxunn indicate the steps in the cross section calcula- 
tion wliere computations for the given species pair are specified. 
The steps defined in the compiled-in data for the standard spec- 
ies are summarized in the cross section edit, figure 16. In the 
many cases whore the final column contains no entry for a species 
pair, one of the default options is used as indicated in the Notes. 
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TABLE XVII 


SOURCES OF CROSS SECTION DATA 


Species Species Notes Computation 

Indices Names Steps 


1-1 

e** 

- e" 

A1 


1-2 

e” 

- N 

A2 

43,50 

1-3 

e“ 

- 0 

A2 

46,50 

1-4 

e” 

- Ar 

A5,A6,A7,A8 

48 

1-5 

e“ 

- N2 

A2 

45,50 

1-6 

e“ 

" ®2 

A2 

47,50 

1-7 

e“ 

- NO 

A2 

44, 50 

1-8 

e“ 

- NO*^ 

A1 


1-9 

e“ 

- N+ 

A1 


1-10 

e" 

- 0+ 

A1 


1-11 

1-12 

e“ 

e" 

“ °2 

A1 

A1 


1-13 

e” 

- CO 2 

A5,A6,A9 

42,49 

1-14 

e” 

- CO 

A5,A6, A9 

41,49 

1-15 

e" 

- CN 

All 

44,49 

1-16 

e“ 

- He 

A8 

55 

1-17 

e“ 

- C 

AlO 

43,49 

1-18 

e 

- c+ 

A1 


1-19 

e“ 

- He+ 

A1 


1-20 

e“ 

- Ar+ 

A1 


1-21 

e“ 

- He(3s) 

A3 3 

55 

1-22 

e“ 

- He (Is) 

A33 

55 

1-23 

e*“ 

- He 2 '^ 

A1 


1-24 

e” 

- Hep 

A3 3 

55 

1-25 

*2“ 

- C0+ 

A1 


1-26 

e~ 

- Ar* (m) 

A14 


1-27 

c" 

- 7U:*(r) 

A14 


1-28 

e“ 

- Ar^"*" 

A1 


2-2 

N - 

■ N 

A2,A4 

12,60 

2-3 

N - 

0 

A2 

14,26 

2-4 

N - 

Ar 

A12 

4 

2-5 

N - 

^2 

A2 

13,26 

2-6 

N - 

O 2 

A2 

18,22,39 

2-7 

N - 

NO 

A2 

13,24,25 

2-8 


NO-''- 

A 14 


2-9 

N - 

N'‘ 

A2, A3 

10,60 

2-10 

N - 

O'*' 

A2 


2-11 

N - 

N.r 

A 14 


2-12 

N - 

°2^ 

A 14 
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TABLE XVII (Cont’d) 


Species Species Notes Computation 

Indices Names Steps 


2-13 

N - CO 2 

A13 


2-14 

N - CO 

A13 


2-15 

N - CN 

All 

13,24,25 

2-16 

N - He 

A13 


2-17 

N - C 

A18,A19 

5r7 

2-18 

N - C'*' 

A14 


2-19 

N - He+ 

A14 


2-20 

N - Ar+ 

A14 


2-21 

N - He(^S) 

A13 


2-22 

N - He (Is) 

A13 


2-23 

N - He2’^ 

A14 


2-24 

N - He, 

A13 


2-25 

N - CO^ 

A14 


2-26 

N - Ar*(m) 

A 13 


2-27 

N - Ar*(r) 

A13 


2-28 

N - Ar,’’’ 

A14 


3-3 

0-0 

A2,A4 

21,61 

3-4 

0 - Ar 

A13 


3-5 

0 - N2 

A2 

18,38 

3-6 

0 - O 2 

A2 

22 

3-7 

0 - NO 

A2 

18,22,39 

3-8 

0 - NO"*" 

A14 


3-9 

0 - N’*' 

A2 


3-10 

0-0"^ 

A2,A3 

11,61 

3-11 

0 - N„+ 

AM 


3-12 

0 - O 2 + 

AM 


3-13 

0 - CO 2 

A13 


3-14 

0 - CO 

A13 


3-15 

0 - CN 

All 

18,22,39 

3-16 

0 - He 

A13 


3-17 

0 - C 

A18,A20 

5,8,9 

3-18 

0 - C+ 

AM 


3-19 

0 - He+ 

1 

AM 


3-20 

0 - Ar 

AM 


3-21 

0 - Hc(^S) 

Aa.3 


3-22 

0 - He(^'-S) 

A13 


3-23 

0 - He 

AM 


3-24 

0 - Hcj 

A13 


3-25 

0 - CO'"’' 

AM 


3-26 

0 - Ar*(m) 

A13 
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TABLE XVII (Cont'd) 


Species 

Indices 

Species 

Ncimes 

Notes 

Computation 

Steps 

3-27 

0 - Ar* (r) 

A13 


3-28 

0 - AT 2 + 

A14 


4-4 

Ar - Ar 

A17,A4 

1,62 

4-5 

Ar - N 2 

A12 

3 

4-6 

Ar - O 2 

A13 


4-7 

Ar - NO 

. A13 


4-8 

Ar - NO'*' 

A14 


4-9 

Ar - N^ 

A14 


4-10 

Ar - O"^ 

A14 


4-11 

Ar - Np+ 

A14 


4-12 

Ar - 02 '*' 

A14 


4-13 

Ar - CO 2 

A13 


4-14 

Ar - CO 

A13 


4-15 

Ar - CN 

A13 


4-16 

Ar - He 

A13 


4-17 

Ar - C 

A13 


4-18 

Ar - C+ 

A14 


4-19 

Ar - He+ 

A14 


4-20 

Ar - Ar"*" 

A17 

2,62 

4-21 

Ar - He(3s) 

A13 


4-22 

Ar - He(-^S) 

A13 


4-23 

Ar - H 62 '*' 

A14 


4-24 

Ar - He 2 

A13 


4-25 

Ar - co"^ 

A14 


4-26 

Ar - Ar*(m) 

A13 


4-27 

Ar - Ar* (r) 

A13 


4-28 

Ar - AT 2 '*’ 

A14 


5-5 

N 2 - N 2 

A2 

17,38 

5-6 

N 2 - O 2 

A2 

20 

5-7 

N 2 - NO 

A2 

17,19,37 

5-8 

N„ - NO+ 

A14 


5-9 

N, - N+ 

A14 


5-10 

N 2 - 0 

A14 


5-11 

Nj - N 2 + 

A 14 


5-12 

N 2 - 0 + 

A 14 


5-13 

N2 - C 62 

A22,A24,A26 

16,34 

5-14 

^2 “ 

A22,A23 

15,29 

5-15 

N 2 - CN 

All 

17,19,37 

5-16 

N„ - lie 

A13 


5-17 

N, - C 

A13 


5-18 

N2 - c+ 

A14 
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TABTE XVII fCont'd) 


Species 

Species 

Notes 

Computation 

Indices 

Names 


Stens 

5-19 

N„ - He"^ 

A14 


5-20 


A14 


5-21 

N2 - He(-^S) 

A13 


5-22 

N - He(ls) 

A13 


5-23 

N2 - He2‘*' 

A14 


5-24 

N2 - He. 

A13 


5-25 

N2 - C0+ 

A14 


5-26 

N2 - Ar* (ni) 

A13 


5-27 

N. - Ar*(r) 

A13 


5-28 

N2 - Ar2' 

A 14 


6-6 

O2 - O2 

A2 

23 

6-7 

6-8 

O9 - NO 
0. - NO'*' 

A2 

A14 

19,23,39 

6-9 

6-10 

1 1 

<N CM 

0 0 

AM 

AM 


6-11 

0 - N2+ 

AM 


6-12 

°2 - 02-^ 

AM 


6-13 

0^ ^^2 

A22,A24 

17,35 

6-14 

O2 - CO 

A22,A24 

15,30 

6-15 

O2 - CN 

All 

19,23,39 

6-16 

Oj - He 

A13 


6-17 

^2 “ 

A13 


6-18 

6-19 

6-20 

0. - C'^ 

0. - He"^ 
^ + 

®2 ” 

AM 

AM 

AM 


6-21 

0 - He(^S) 

A13 


6-22 

O2 - He (Is) 

A13 


6-23 

O2 - He2'‘‘ 

AM 


6-24 

O2 - He 
°2 " 

A13 


6-25 

AM 


6-26 

O2 - Ar* (m) 

A13 


6-27 

0- - Ar*(r) 

A13 


6-28 

O2 - Ar/ 

AM 


7-7 

NO - NO 

A2 

17,19,20, 

7-8 

NO - NO'*' 

A2 

59,53 

7-9 

NO - N^ 

AM 


7-10 

NO - O’’' 

AM 


7-11 

NO - N2''' 

AM 


7-12 

NO - 02’'' 

AM 


7-13 

NO - CO2 

A30 

16,32 
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TABLE XVII fConfd) 


Species Species Notes 

Indices Names 


7-14 

NO - CO 

A29 

7-15 

NO - CN 

All 

7-16 

NO - He 

A13 

7-17 

NO - C 

A13 

7-18 

NO - C'*' 

A14 

7-19 

NO - He"^ 

A14 

7-20 

NO - Ar+ 

A14 

7-21 

NO - He(3s) 

A 13 

7-22 

NO - He (Is) 

A13 

7-23 

NO - He2‘‘‘ 

A14 

7-24 

NO - He, 

A13 

7-25 

NO - CO+ 

A14 

7-26 

NO - Ar* (m) 

A13 

7-27 

NO - Ar*(r) 

A13 

7-28 

NO - Ar2'^ 

A14 

8-8 

NO'*' - NO+ 

A1 

8-9 

NO+ - N+ 

A1 

8-10 

NO+ - 0+ 

A1 

8-11 

NO+ - N,'^ 

A1 

8-12 

NO"^ - 02’*' 

A1 

8-13 

NO+ - CO 

A 14 

8-14 

NO’*' - CO^ 

A14 

8-15 

NO"*" - CN 

A14 

8-16 

no’*’ - He 

AM 

8-17 

NO+ - C 

AM 

8-18 

NO’*’ - C+ 

A1 

8-19 

NO+ - He+ 

A1 

8-20 

NO’*’ - Ar+ 

A1 

8-21 

NO+ - He(3s) 

AM 

8-22 

NO+ - He(^S) 

AM 

8-23 

NO’*' - He2’*’ 

A1 

8-24 

NO’I' - He, 

AM 

8-25 

no’*' - CO+ 

A1 

8-26 

NO"*- - Ar*(m) 

AM 

8-27 

NO+ - Ar*(r) 

AM 

8-28 

no’*' - Ar,-*- 

A1 

9-9 

N’*’ - N+ 

A2 

9-10 

N^' - O’*’ 

A1 

9-11 

N-*- - N^'*' 

A1 

9-12 

N+ - 

A1 


-lb3- 


Computation 

Steps 

15,31 

40 



T ABLE XVII (Cont'd) 


Species 

pecies 

Indices 

Names 

9-13 

- CO 2 

9-14 

N+ - CO ’ 

9-15 

N"^ - CN 

9-16 

N'*' - He 

9-17 

N+ - C 

9-18 

N+ - C+ 

9-19 

N'*’ - He"^ 

9-20 

- Ar'*’ 

9-zl 

N+ - He (3s) 

9-22 

n’*' - He (Is) 

9-23 

N+ - He^”^ 

9-24 

9-25 

N+ - He, 
N+ - CO"^ 

9-26 

N+ - Ar*(m) 

9-27 

N+ - Ar*(r) 

9-28 

N+ - Ar,+ 

10-10 

0+ - 0+ 

10-11 


10-12 

10-13 

•f* “-I- 

0+ “ 02-" 
0^ - CO, 

10-14 

0*^ - CO 

10-15 

0+ - CN 

10-16 

0+ - He 

10-17 

0+ - C 

10-18 

0+ - C"^ 

10-19 

O'*’ - He+ 

10-20 

0+ - Ar+ 

10-21 

O"^ - He(^S) 

10-22 

0"^ - He (Is) 

10-23 

0 + - He/ 

10-24 

0+ - He, 

10-25 

0+ - C0+ 

10-26 

0 + - Ar*(ra) 

10-27 

0 ^ - Ar*(r) 

10-28 

11-11 

11-12 

o'*" - Ar,+ 

^ 2 : - 

n/ - O 2 

11-13 

Nj+ - CO 2 

11-14 

N2'^ - CO 

11-15 

IT 2 + - CN 

11-16 

N 2 + - Ho 


No.teri Computation 

Steps 

A 14 

A 14 

A 14 

A 14 

A14 

A1 

A1 

A1 

A 14 

A14 

A1 

A14 

A1 

A 14 

A 14 

A1 

A1 

A1 

Ai 

A 14 

A 14 

A 14 

A 14 

A14 

Al 

Al 

Al 

A 14 

A 14 

Al 

A 14 

Al 

A14 

A 14 

Al 

Al 

Al 

A 14 

A 14 

A 14 

A 14 
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TABLE XVII (Cont'd) 


Species Species Notes Computation 

Indices Names Steps 


11-17 


- c 

A14 


11-18 


- c+ 

A1 


11-19 


- He'*' 

A1 


11-20 


- Ar+ 

A1 


11-21 


- He(^S) 

A14 


11-22 


- He (Is) 

A14 


11-23 

N 2 + 

- He2-*- 

A1 


11-24 


- Hep 

A14 


11-25 

N 2 + 

- CO^ 

A1 


11-26 


- Ar* (m) 

A 14 


11-27 

N 2 + 

-Ar* (r) 

A14 


11-28 

12-12 

N 2 : 

O 2 

- Arp"^ 

- °2^ 

A1 

A1 


12-13 

O 2 : 

- CO 2 

A14 


12-14 


- CO 

A14 


12-15 

°2. 

- CN 

A14 


12-16 


- He 

A14 


12-17 

®2'^ 

- C 

A 14 


12-18 

°2 

- C+ 

A1 


12-19 


- He+ 

A1 


12 20 


Ar+ 

A1 


12-21 


- He(^S) 

A 14 


12-22 

0^: 

- He (Is) 

A14 


12-23 

02*' 

- He2+ 

A1 


12-24 

12-25 

O 2 : 

O 2 : 

- HCp 

- co"^ 

A14 

A1 


12-26 

°2 

- Ar* (m) 

A14 


12-27 

02'" 

- Ar*(r) 

A 14 


12-28 

02'^ 

- Arp"*" 

A1 


13-13 

co„ 

- CO 2 

A22,A25,A27 

16,33 

13-14 

C 02 

- CO 

A31 

15,28 

13-15 

C 02 

- CN 

All 

16,32 

13-16 

COo 

- He 

A13 


13-17 

C 02 

- C 

A 13 


13-18 

CO 

- C'*" 

A 14 


13-19 

coi 

- Hc+ 

• 4 


13-20 

13-21 

COy 

COp 

- Ar+ 

- Hc(^S) 

* -■ ** 

A 13 


13-22 

CO- 

- He (Is) 

A 13 


13-23 

C 02 

- IIC 2 + 

A14 
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TABLE XVII fConfd) 


Species 

Indices 

Species 

Names 

Notes 

Computation 

Steps 

13-24 

CO 2 - He^ 

A13 


13-25 

CO 2 - 

A14 


13-26 

CO 2 - Ar*(m) 

A13 


13-27 

CO 2 - Ar* (r) 

A13 


13-28 

CO 5 - PiTj* 

A14 


14-4 

CO - CO 

A22,A23,A28 

15,27 

14-15 

CO - CN 

All 

16,31 

14-16 

CO - He 

A13 


14-17 

CO - C 

A13 


14-18 

CO - 

A14 


14-19 

CO - He+ 

A14 


14-20 

CO - Ar+ 

A14 


14-21 

CO - He(3s) 

A13 


14-22 

CO - He (Is) 

A13 


14-23 

CO - He 2 **' 

A14 


14-24 

CO - He 2 

A13 


14-25 

CO - CO^ 

A 14 


14-26 

CO - Ar*(m) 

A13 


14-27 

CO - Ar*(r) 

A13 


14-28 

CO - Ar 2 + 

A14 


15-15 

CN - CN 

All 

40 

15-16 

CN - He 

A13 


15-17 

CN - C 

A13 


15-18 

CN - C+ 

A 14 


15-19 

CN - He"^ 

A14 


15-20 

CN - Ar'*' 

A 14 


15-21 

CN - He(3s) 

A 13 


15-22 

CN - He (Is) 

A13 


15-23 

CN - He 2 + 

A 14 


15-24 

CN - He 

A13 


15-25 

CN - CO"^ 

A 14 


15-26 

CN - Ar*(m) 

A 13 


15-27 

CN - Ar* (r) 

A 13 


15-28 

CN - Ai' 2 + 

A14 


16-16 

He - He 

A34,A4 

51,57 

16-17 

He - C 

A 13 


16-18 

He - C'*" 

A 14 


16-19 

He - Hc'^ 

A3 5 

53,57 

16-20 

He - Ar’’’ 

A 14 


16-21 

He - IIe(3s) 

A3 3 

51,57 
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TABLE XVII (Cont'd) 


Species 

Indices 

Species 

Names 

Notes 

Computation 

Steps 

16-22 

He - He (^S) 

A33 

51,57 

16-23 

He - He 2 '*' 

A33 

53 

16-24 

He - He^ 

A13 


16-25 

He - CO'T 

A14 


16-26 

He - Ar* (ra) 

A13 


16-27 

He - Ar*(r) 

A13 


16-28 

He - Aro'*' 

A14 


17-17 

C - C 

A18,A21,A4 

5,6,60 

17-18 

C - 

A32 

59,60 

17-19 

C - He'*' 

A 14 


17-20 

C - Ar"*" 

A 14 


17-21 

C - He(^S) 

A13 


17-22 

C - He (Is) 

A13 


17-23 

C - He2'*- 

A14 


17-24 

C - He 

A13 


17-25 

C - 

A 14 


17-26 

C - Ar*.{m) 

A13 


17-27 

C - Ar*.(r) 

A13 


17-28 

C - At2+ 

A 14 


18-18 

C+ - C+ 

A1 


18-19 

C+ - He+ 

A1 


18-20 

C’*' - Ar+ 

A1 


18-21 

C+ - HeC^s) 

A14 


18-22 

C+ - He (Is) 

A 14 


18-23 

C+ - He2'‘- 

A1 


18-24 

C"*" - He, 

AM 


18-25 

C'*' - CO^ 

A1 


18-26 

C+ - Ar*(m) 

AM 


18-27 

€■*■ - Ar* (r) 

AM 


18-28 

C+ - Ar,+ 

A1 


19-19 

He"*" - He+ 

A1 


19-20 

He"^ - Ar"*' 

A1 


19-21 

He"^ - He (3s) 

A3 3 

53,57 

19-22 

He+ - He(^S) 

A33 

53,58 

19-23 

Hc;j; - He 2 '*' 

A1 


19-24 

He - He 2 

A33 

53 

19-25 

He'*' - 00“^ 

A1 


19-26 

He"*" - Ar*(m) 

AM 


19-27 

He+ - Ar*(r) 

AM 


19-28 

He+ - Ar 2 ''' 

A1 
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TABUS XVII (Cont*d ) 


Species 

Species 

Notes 

Computation 

Indices 

Ncunes 


Stops 

20-20 

Ar'*' - Ar+ 

A1 


20-21 

Ar"*" - He(3s) 

A14 


20-22 

Ar‘*‘ - He (Is) 

A14 


20-23 

Ar"^ - He2'*' 

A1 


20-24 

Ar"*" - He_ 
Ar+ - CO^ 

A14 


20-25 

A1 


20-26 

Ar'*' - Ar*(m) 

A14 


20-27 

Ar'*^ - Ar*(r) 

A14 


20-28 

Ar"*" - Ar 2 '*' 

A1 


21-21 

He (3s) - He(^S) 

A3 3 

51,58 

21-22 

He (3s) - He (Is) 

A3 3 

51,58 

21-23 

He (3s) - He 2 ‘‘‘ 

A33 

54 

21-24 

He (3s) - He^ 
lie (3s) - CO+ 

A13 


21-25 

A14 


21-26 

He (3s) - Ar*(m) 

A13 


21-27 

He (3s) - Ar*(r) 

A13 


21-28 

He (3s) - AT2+ 

A14 


22-22 

He (Is) - He (Is) 

A33 

52,58 

22-23 

He (Is) - He 2 ‘*‘ 

A3 3 

54 

22-24 

He (Is) - He, 
He (Is) - CO^ 

A13 


22-25 

A14 


22-26 

He(^S) - Ar*(m) 

A13 


22-27 

Ile(-^S) - Ar*(r) 

A 13 


22-28 

He (Is) - Ar 2 ’^ 

AM 


23-23 

He2+ - He2+ 

A1 


23-24 

23-25 

He,+ - He, 
He,'*’ - CO^ 

A33 

A1 

54 

23-26 

He,'*’ - Ar* (m) 

AM 


23-27 

He 2 '*' - Ar* (r) 

AM 


23-28 

Ho 2 + - Ar,+ 

A1 


24-24 

Hc 2 - He, 

A3 3 

52 

24-25 

He, - CO+ 

AM 


24-26 

He, - Ar*(m) 

A13 


24-27 

He, - Ar*(r) 

A13 


24-28 

He^ ~ Ar,'*' 

AM 


25-25 

CO^ - co '^ 

A1 


25-26 

CO'*’ - Ar*(m) 

AM 


25-27 

CO'*’ - Ar*(r) 

AM 


25-28 

CO+ - Ar,+ 

A1 


26-26 

Ar* (m) - Ar* (m) 

A15 

1,62 
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TABLE XVII fCont'd) 


Species Species Notes Computation 


Indices 

Names 



Steps 

26-27 

Ar* (m) 

- Ar*(r) 

A13 


26-28 

Ar* (m) 

- Ar,"*"- 

A14 


27-27 

Ar* (r) 

-Ar*(r) 

A15 

1,62 

27-28 

Ar* (r) 

- Ar + 

A14 


28-28 

At2+ - 

Ar2^^ 

Al 


Notes to Table 

XVII 




Al. Default option ;uses 

effective 

Coulcxnb cross 

sections cal- 


culated from equations (100) of Volume I (ref. 1) . 

A2. Reference 22. 

A3. Reference 23. 

A4. The self-diffusion coefficient for atoms is set equal to 

the atom-ion charge exchange cross section in calculating 
the internal thermal conductivity, in order to account ap- 
proximately for the effects of resonant excitation energy 
exchange (see ref, 22) . 

A5. Effective cross sections are used, based on curvefit to 
mobility data. 

A6. Reference 24. 

A7. References 25, 26, 27. 

A8. Referenje 28. 

A9. Reference 29, 

AlO. For electron-car]?on atom collisions, we a^ume a constant 

collision cross section - tt il. 2) _ ^ ^ jLq 16 

cm^, in analogy to the case of e-N. This value appears to be 
consistent v.’ith available theoretical estimates (see ref. 30) . 


'r'-i 
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1 % 


All. 

A12. 

A13. 

A14. 

A15. 

A16. 

A17. 

A18. 

A19. 

A20. 

A21. 

A22. 

A23. 

A24. 

A25. 

A26. 

A27. 


I 


The collision cross sections for CN have arbitrarily been set 
equal to the corresponding cross sections for NO, 

Curvefit to data of reference 31. 

Default option; cross sections calculated from the empir- 
ical mixing rule, equation (17). 

Default option; cross sections arbitrarily set equal to 
the estimated values for N - 0+ collisions. See equation 
(23). 

Cross sections of excited argon arbitrarily assiimed equal 
to those for the ground state atoms. 

Curvefit to data of reference 32. 

Reference 33. 

Cross sections estimated from an approximate perfect pair- 
ing calculation, v;ith the para^meters determined from avail- 
able spectroscopic data and bv .analogy vi th the oxygen and 
nitrogen results (refs. 34-37) . 

Reference 38. 

Reference 39. 

References 40, 41. 

Cross sections obtained by fitting experimental tremsport 
property data below about 1000°K and e>:trapolating to higher 
temperatures assuming the same temi:)eraturc dependence as 
for collisions. 

Reference 42. 

References 43, 44. 

Reference 45. 

Reference 46. 

Reference 47. 
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A28. Reference 48, 

A29. Mean of CO-N 2 and CO-O 2 cross sections. 

A30. Mean of CO 2 -N 2 and CO 2 -O 2 cross sections. 

A31. Mean of CO-CO and CO 2 -CO 2 cross sections. 

A32. For C-C"^ collisions, the charge exchange cross section is 
arbitrarily set equal to the N-NT*' cross section, while the 
gas kinetic cross section is set equal to the N-O"^ value. 

A33. Cross sections of excited He arbitrarily assumed equal to 
those for ground -state He. 

A34. Reference 49. 

A35. Reference 50. 


The precoded data for NNKQ, NNQ, Ilm, JJm, Wm, and ISEQ 
can all be read or inferred from the cross section edit, figure 
16. The steps are performed in the order listed. Tlie first col- 
umn in figure 16 is a counter for the steps in this order. The 
second column gives the values of the sequencing array, ISEQ. 

For the precoded data, ISEQ(L)=L. The third colximn gives the 
values of the option index, KKQ. The columns headed V(l),..., 
W(5) list the parameter values for each step. Finally, the last 
column gives the pairs of species to which the step is applied. 

In some cases (e.g, , steps 51 and 52) , a step is repeated to cir- 
cumvent the limit of five species pairs per step. 

Table XVIII summarizes the precoded contents of the TL, 
0MEGA1, ASTAR, and BSTAR arrays. In the many cases in v/hich 
"cross section table" is entered under "Remarlcs" the data are 
talaulated cross sections for use with the option ICKQ=5. In these 
cases, the TL, 0MEGA1, ASTAR, and BSTAR arrays contain data as 
specified above in the discussion of this option. For the indi- 
ces 996-999 containing data for the power law interaction (KKQ-6) , 
no data arc stored in TL, and the data in the other arrays are as 
specified in cqviations (14). 
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TABLE XVIII (Concl'd) 


f 


1 




■ 











t=d 

see 
















H rj 

c 





*«— sM % 

0 






•p 





* 

■p 





i ic^ “ 

o 

to 

VD 

CD 

VD 


^ C tJi {Ji 

u 

Q) 

II 

II 

II 


•H c C 

-P 





II T3 -H H 

C 

ff 

a 



C XI TI 

•H 

1:^ 




* 0 c c 




PM 

(0 

H Ot 0 O 

c 




Ai 

W Oi D< 

0 



V 


4-1 (U U) w 

•H 

c 

c 

C 

to 

O S-( (U 0) 

1 

0 

0 

0 

B 

P P H 

rP 

•H 

•H 

•H 

0) 

W O P }-4 

to 

p 

p 

P 

Pi 

a) u 0 o 

P 

o 

u 

o 

u u 

4J 

to 

to 

(d 


iH tn 

3 

p 

p 

p 


lO C w 01 

(U 

<u 

(!) 

0) 


> -H C C 

C 

p 

P 

p 


(d 'r| T-I 


c 

C 

c 


m 4-> to (0 

P 


■H 



C C 4J +> 

0 





•H o C C 

4-1 • 





(0 U 0 0 

S 

to 

to 

to 


-P u o 

-P m 

fH 

fH 

I— i 


C H 

r-| (N 





0 < k: Di 

P 

p 

p 

p 


u o <c < 

to 

<u 

0) 

QJ 


W Eh E-I 

4h • 





1-3 S U) W 

0) C5^ 

0 

o 

0 


Eh is;. C W 

« 0) 

Oi 

ft 

ft 

0 






0 

1 

1 

1 — 1 

o 

in 

-p 

1 

1 

1— I 

iH 


U) 







1—1 




(0 


fd 




C 


•H 




0 


-p 




•H 






p 


<u 




u 


-p 




tO 


o 




p 






0) 





p 






c 


(N 




•H 


1— 1 





c 

1 




u 

0 

VD 

c 

c 

c 

I 

•H 

•w 

0 

0 

0 

o 

4J 


•H 

*H 

•H 


P 

m 

■p 

p 

p 

X! 

*H 

CD 

u 

o 

0 

C 


C 

fd 

fO 

to 

lO 

u 

0 

p 

p 

p 


U) 

b 

0) 

G) 

0) 

% 

o 

1 

p 

P 

p 

u 

Q 

X5 

c 

c 

c 

1 


P 

■H 

•H 

•H 

a 


to 






c 

+ 

+ 




c 

o 

O 

a 

u 


Q) 

1 

1 

1 

1 




o 

a 

u 


r- 






n 





X 

in 





0 

1 





V 

H 

IXI 


CO 

cn 

c; 

O 

cn 

CTi 

cn 

ov 

H 

in 

cyi 

cn 

cn 



cn 


-lCS-1- 











r 


! 


I 


4.7 Nozzle and Channel Geometries 

The geometric profiles for ten standard NASA Johnson Space 
Center nozzles and two rectangular channels are compiled into 
NATA. These data are indexed as explained in the definitions of 
NOZZLE and NPR0FL in Section 2.3 (Group 4). NATA users at other 
laboratories can advantageously replace these data with geometric 
descriptions applicable to their own facilities. 

The geometry of an axisymmetric nozzle is defined by a single 
profile. That of a rectangular channel requires two profiles for 
its description. As explained in Section 4.3 of Volume I (ref. 1) , 
each profile is represented by an analytical curvefit containing 
up to 12 sections. The sections are joined end to end with value 
and slope continuity. At least two sections must be upstream of 
the throat, and at least two must lie donwstream. The throat 
must be a section boundary. Each section in a profile fit may have 
one of three forms: 

(1) Straight Line (ISHAPE = 1) 

y (x) = P]^ + P 2 X (24a) 


(2) Circular Arc Convex Downward (ISHZU?E = 2) 

y(x) = -/pg^ - (x-P 2 )^ (24b) 

(3) Circular Arc Convex Upward (ISHAPE = 3) 

y(x) = Pj^ + /Tp _ (x-p 2 )^ (24c) 


In the second and third forms, Pg is the radius of the circular 
arc and (P 2 » !*]_) the x and y coordinates, respectively, of 

the circle center. The geometric summary in figure 4 gives an 
illustration of a NATA profile curvefit. The inlet position 
listed is the starting point for boundary layer calculations. 
The coluirui headed "ATPI(J)" contains the downstream boundaries 
of the sections. The parameters P,, P 2 , P 3 arc listed as 
PARAM(1,J), PARAM(2,J), PARAH(3,J), respectively 
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The precoded profile data are stored in an array ZPRP (I, NOZZLE ) , 
dimensioned (64,20). The dimensions allow as many as 20 compiled-in 
profiles. For convenience in adding to or altering the precoded 
data, ZPRP is equivalenced to 20 singly dimensioned arrays (ZPl(I), 
ZP2(I), etc., as follows: 

ZPl(l) equivalent to ZPRP(1,1) 


ZP20(1) 


equivalent to ZPRP(1,20) 


Thus, ZPl contains the precoded data for N0ZZLE = 1, ZP2 those for 
N0ZZLE = 2, and so forth. The data in each ZPn array are as fol- 
lows: 

Array Element Definition 


ZPn(l) 


Throat radius (cm) 


ZPn(2) 


Starting point for boundary layer calculations 
(negative value in cm upstream of the throat) , 


ZPn(3) 


Number of profile sections upstream of the 
throat. * 


ZPn (4) 


Number of profile sections downstream of the 
throat . * 


ZPn (4+ I) 


For I = 1 to 12, ISKAPE value for the Ith pro- 
file section.* 


SPn(16+I) For I = 1 to 11, the downstream boundary of 

the Ith profile section in centimeters from 
the throat (negative i 'istreeim) . 


*Thcse integer data arc stored as real values, increased by 0.1 
in each case to ensure reliable rounding down to the original 
integer values when the data arc used. 
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Array Element 
ZPn(?4+3I+K) 


Definition 


For K = 1 to 3 and .T = 1 to 12, the Kth par- 
ameter value (see eqs. 24) for the Ith 
profile section. The parameters having length 
dimension are given in centimeters. 

ZPn(64) Facility name (Hollerith data) 


The precoded data for standard channels are stored in an 
array CP(I,ICHMI), dimensioned (5,5). CP is equivalenced to five 
singly dimensioned arrays, CPl(I), CP2(I),..,, CP5(I). each of 
which contains or can contain the data for a channel. For example, 
CPI (I) contains the data for ICHA.J = 1. The contents of these 
arrays are defined as follows: 

Array Element Definition 

CPn(l) NPR0FL(1), the index specifying the precoded 

data for the first profile of the channel; 
these data are stored in ZPRP(I,NPR0FL(1) ) 
for I = 1 to 64. 

CPn(2) NPR0FL(2), the index specifying the precoded 

data for the second profile of the channel. 

CPn(3) Channel name (Hollerith data). 

CPn(4) Index (1 or 2) specifying the profile which 

diverges from the axis least rapidly downstream 
of the throat (corresponds to NBL in Section 
2.3, Group 4) . 

CPn(5) Facility name (Hollerith data). 


NATA includes pre-coded data for tv;o channels, as indicated in the 
definition of ICii^^N in Section 2.3 (Group 4). 

Figures 47 to 59 arc plots showing the throat regions of all 
of the iDrccodcd profiles. Each of these figures shov;s a 15.24-cm 
(6-inch) long portion of a profile. The profile actually continues 
indefinitely far to the right and left of the figure boundaries; 
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FIGURE 47 - PROFILE FOR DCA 1.90-cm THROAT (N0ZZLE=1) 


DCA 0.75-INCH THROAT 
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FIGURE 48 - PROFILE FOR DCA 3.81-cm THROAT (N0 ZZt,e=2) 


DC A 1.5- INCH THROAT 
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FIGURE 49 - PROFILE FOR MRA 5.72-cm THROAT (N{2fZZLE-3) 


MRA 2. 25- INCH NOZZLE 
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FIGURE 50 - PROFILE FOR MRA 2. 54 -cm THROAT (N02ZLE=4) 


MRA 1.0- INCH THROAT 

«.o 






FIGURE 51 - PROFILE FOR EOS 0.81-cm THROAT (N0ZZLE=5) 


EOS 0.32- INCH THROAT 
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FIGURE 52 - PROFILE FOR EOS 1.97-f'm THROAT (N0ZZLE=6) 


EOS 0.775- INCH THROAT 
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FIGURE 53 - PROFILE FOR MRA 1.9n-cm TimOAT (N0ZZLE=7) 


MRA 0.75- INCH THROAT 
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FIGURE 54 - PROFILE FOR MRA 3.81-cm THROAT (N0ZZLE-8) 


MRA 1.5- INCH THROAT 
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FIGURE 55 - PROFILE FOR 10 5.72-cm THROAT (N0ZZLE=9) 
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FIGURE 56 - PROFILE FOR EOS 2.77-cn. HIROAT (N/)'ZZLE-10) 
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FIGURE 57 - FIRST PROFILE FOR Tl2 AND T22 CHANNELS (NPR0FL=L1) 


DC A T12 CHANNEL. PROFILE 1 
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FIGURE 58 - SECOND PRCFIL.; FOR T12 CHANNEL (NPR0FL=12) 
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FIGURE 59 - SECOND PROFILE FOl 


DC A T22 CHANNEL. PROFILE 2 
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the program uses as much of the mathematically defined profile 
as it needs in each problem. The plots shown in figures 47—59 
were produced by the NOZFIT code, an a\uiiliary computer progrcim 
for setting up NATA-type profile curvefits from data provided by 
nozzle design drav/ings. A user's manual for NOZFIT is included 
in the present report (Appendix D) , 

The profiles as used in NATA, and as shown in the figures, 
differ in several respects from the profiles of the actual noz- 
zle hardware: 

(1) The NATA profiles expand conically to the left 
(in the upstream direction) , while the actual 
nozzles have finite-diameter plenum or arc- 
chaniloor radii. 

v2) Sharp corners in the actual nozzle profxles are 

rounded in the NATA fits, to provide continuity of 
the profile slope (dy/dx) as required by the code. 

A standard 0.127 cm (50-mil) rounding radius is 
used, except in cases where a larger radius ?' >.s 
proved necessary for code reliability. 

(3) In the fits for many ol the nozzles, sections of 
constant radius near the throat arc rexoresented 
as conical, usually \;ith a 3° con/ergcnce half 
angle. Tliis is done ro avoid instabilities in 
the nonequilibriuni solution (Section 4.3, Volume 
I). 
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APPENDIX A 


KEACTION DATA FOR THE HELIUM AND ARGON MODELS 


Tills appendix documents the reaction systene. and the elec- 
tronic nenequi librium parameters assumed in the standard gas 
models for helium (IG7\G = 4) and argon (IGAi- = 3) . In addition 
to chemical nonequi librium these models include effects of nen- 
equilibrium excitation of the gases by treating each of the im- 
portant excited states as a separate species. Approximate reac- 
tion parauTictors for the important reactions ameng these states 
arc then obtained from a survey of the available literature. 

Tlie species and parameter values used in the models arc 
given in Tables XIX, XX, X>CI, and XXII. The reasons for choice 
of the tabulated valx’es arc discussed belov/. 


A.l Helium Model 


E last i c collisions — The simplest type of collision process 
occurring .•>n a gas is the elastic collision in v;hich kinetic en- 
ergy is transferred from one particle to another without any 
change in the internal structure or excitation of the particles. 
Although such collisions ob' iously do not contribute to the spec- 
ies production term I'j in equation (32ia) of Volume I, the kinet- 
ic energy transferred bctv;een elections and heavy particles n 
elastic collisions can be important in detc .lining the net cn^''"f/y 
gain term q^ for the eleccron gas. Under the assumption that the 
electrons and heavy particles have Maxwellian v_!.ocity distribu- 
tions corresponding to tlic temperatures end T, respectively, 
it can bo shovm (ref. 51) that the contribution to the electron 
energy gain term q^ in equation (321c) (Voluiuo I) due to clastic 
collisions is given to a vein good approximation by tk.i formula* 
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= X 

3=2 
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j , cia£ 


l1 > c.las 
Nn 


(25) 


*It is assumed in ec aat.ion (25) and throughout this Appendi:: that 
the species j ' 1 rc ros^onts the oj.cctrons. 
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•For molecular species, the tabulated ground-state energy is the energy of the 
lowest vibrational level, v •> 0. 
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where the sum extends over all heavy particles j present in the 
gas. 


Sj.elas = f Rq “ ’'e) 


(26) 


represents the mean energy gained by the electrons in Nq elastic 
collisions with particles of the jth species and 


N 


ej,elas ~ 






= N^Tc . 

0 e] ,elas 


(/» 


(27) 


represents the n\imber of elastic collisions occurring between 
electrons and particles of the jth species per unit volvune per 
unit time. Here 

«j = Nq p Vj (28) 


is the number of particles of the jth species per unit volume, 
T^jsnjkgj elas/^0^® momentvuu transfer collision frequency 
for elastic collisions between electrons and particles of the 
jth species, and tlie reaction rate k^j^elas given by 


^ej,elas/^0 = 


3 V rrwe 


_ (1,1) 

®ej 


(29) 


where 

Q^. r O’ .'^(w)dv/ (30) 

2(kT^)3 Jo ^ 

is the Maxwe ll~averaged momcntvun transfer cross section for elas- 
tic collisions between electrons and particles of species j at the 
temperature and o ^(w) is the actual momentum transfer cross 
section as a function 5f electron energy w. 

The integral in equation (30) has been evaluated approximately 
for electron-ion collisions in reference 51 to obtain 
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for ions 


in the limit when the logarithm in (31) is much greater than 1, 
where is in °K and n^ in cm“^. This approximation should be 
adequate for most cases of interest in the present study. 

To obtain the cross section for ground-state helium 

atoms, we have evaluated the integ^’al in (30) numerically using 
literature data on the electron-helium momentum rransfer cross 
section ‘^e_fje('^) a function of electron energy. Below an 
electron energy of 5 ev, we have used the recent experimental mea- 
surements of Crompton, Elford, and Robertson (ref, 52) , who give 
values of the electron-helium moment\am transfer cross section 
e-He energy range from 0.008 ev to 6 ev with an esti- 

mated experimental error of about 2 percent. Since accurate ex- 
perimental data on the momentum transfer cross section are not 
presently available for electron energies above 6 ev,we have ob- 
tained the cross sections in this range by i ntegrating the theo- 
retical differential scattering cross sections of LaBahn and Call- 
away (ref, 53) over the scattering angle. On the basis of corapari- 
ijons v.’ith experimental data at both lov/ and high ( ;^ 100 ev) ener- 
gies, LaBahn and Callaway estimate that their cross sections should 
Le accurate to within about 5 % in this energy range. 


Figure 60 shows the electron-helium momentum transfer cross 
sections obtained from the data of Crompton, Elford, and Robert- 
son (ref. 52) and from LaBahn and Callaway (ref, 53) as a function 
of electron energy. The results of these two studies are in good 
agreement for electron energies near the upper limit of Crompton's 
measurements at 5 ev. For the present calculations, we have adopted 
the cross section values of Crompton, Elford, and Robertson below 
5 ev and those of LaBahn and Callaweiy above 5 ev, as shown by the 
solid curve in figure 60, Using this adopted cross section, the 
integral in equation (30) v/as evaluated numerically to obtain the 
Max\\/ell-avoraged electron-helium cross section shown in 

figure 61 as a function of electron temperature. Since numerical 
errors in the integration process should be negligible, the ac- 
curacy of the averaged cross section shown in figure 61 is 


108 - 



FIGU^ 60 - MOMENTUM TRANSFER CROSS SECTION FOR ELECTRON-HELIUM COLLISIONS 



ELECTRON ENERG' 


FIGURE 61 - MAXWELL AVERAGED ELECTRON-HELIUM MOMENTUM TRANSFER CROSS SECTION 
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determined by the accuracy of the original cross section data 
used in the computations, and should be within the 2 percent ex- 
perimental error of Crompton's measurements over most of the tem- 
perature range from 200 to 50,000°K shown in the figure. 

Sufficient data to determine accurate elastic collision cross 
sections for other neutral helium species (i.e,, molecules or 
excited atoms) do not appear to be available at present, and in 
our calculations we have simply taken the elastic cross sections 
for such species equal to the cross section for ground-state hel- 
ivim atoms sho\^?n in figure 61. Although not accurate in detail, 
this approximation should have a negligible effect on the final 
results of the calculations, since, for excited species, the energy 
transfer due to inelastic processes should almost always be much 
greater than the elastic losses. Thus whenever the concentration 
of excited species in the gas becomes large enough to significantly 
affect the average collision cross section for the gas in equa- 
tions (25) to (30) , the elastic energy loss term (25) will itself 
become negligible compared to inelastic loss processes in deter- 
mining the overall energy balance for the electron gas. 

In addition to the energy transferred to the heavy particles, 
elastic collisions also result in some energy loss from the elec- 
trons due to free-free radiative processes (bremstrahlung) . In 
the present model these losses are includea as part of the general 
collisional-radiative mechanism discussed belov/, so -that a sepa- 
rate radiative loss term to account for them is not required in 
our treatment of elastic collisions. For the usual experimental 
situation in which the electron thermal energy is small compared 
to the recombination energy, free-free radiation processes account 
for only a small part of the total radiative loss from the gas 
(ref, 54) . 

Collisiona 1-radiative recombination . - It now appears to be 
well established that the recombination of atomic ions ii'i a hel- 
ium plasma occurs primarily by the collisional-radiative process 
suggested by Bates, Kingston, and McWliirter (ref, 55), in which 
electrons rccomlDine first into highly excited atomic states and 
are then stabilized by collisional and radiative transitions to 
lower states. Detailed calculations of the electronic recombina- 
tion rates for this mechanism were carried out by Bates, Kingston, 
and Mcvniirter under the assumption tliat atomic excitation accovmts 
for a negligible fraction of the total gas energy, and using approx- 
imate theoretical values for the required collisional excitation and 
de-excitation cross sections between excited atoms ana electrons. 
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This calculation vas later extended by Bates and Khare (ref. 56) 
to account for the stabilization of excited atoms by collisions 
with ground state neutral atoms, and by Bates, Bell, and Kingston 
(ref. 57) to obtain a mere accurate determination of the popula- 
tion of atoms in metastable excited states. More recently, the 
original calculations of Bates, Kingston, and MeWhirter have beoi 
repeated by Johnson and Hinnov (ref. 58) over a limited range of 
gas conditions using excitation cross sections adjusted to fit 
their experimental spectroscopic data on the population of excited 
state atoms in helivim. The results of these calculations have 
been found to be in reasonably satisfactory agreement with avail- 
able experimental data on the recombination of electrons in helivim 
(refs. 58-60) , within the accuracy of the rather large uncertain- 
ties presently existing both in the experimental data and in the 
cross sections assumed in the theoretical calculations. Although 
these uncertainties have as yet precluded a detailed quantitative 
test of the accuracy of the theoretical predictions, the basic 
correctness of the collisional-radiative model appears to be v/ell 
substantiated by the general agreement between theory and experi- 
ment which has been obtained. 

Although the collisional-radiative model of Bates, Kingston 
and MeWhirter appears to be the most accurate theory presently 
available for treating the recombination of atomic ions in helium, 
there are several Jis^'‘'5*'^antagGs to the direct use of this model in 
the NATA code. First, of course, the model requires rather lengthy 
numerical calculations to determine the electron recombination 
rate for any given set of gas conditions, so that direct use of 
this model in a nonequilibrium flow program such as the NATA code, 
in which reaction rates must be determined at many points, would 
lead to excessively long execution times for the code. Further, 
the collisional radiative model requires input data on a large 
nvinber of excited state excitation and de-excitation cross sec- 
tions which appear to be hnown less accurately at present (ref. 

58) than are the overall electronic recombination rates (ref. 59). 
Thus, to obtain accurate results for the recombination rate from 
the model, it would probably be necessary to carry out a parametric 
study similar to that of Johnson emd Hinnov (ref. 59) in which the 
excited state cross sections were adjusted to obtain the best fit 
between the theoretical predictions and available experimental 
data, and these adjusted cross sections v;ere then used in the the- 
oretical model to predict the electronic recombination rate as a 
function of gas conditions. Such a study would go beyond the scope 
of the present effort; and furthermore, since sufficient data are 
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not presently available to uniquely determine the large nximber of 
adjustable cross sections in the model, the final accuracy which 
could be obtained in the recombination rates by this approach is 
somewhat uncertain. For the present study we have, therefore, 
adopted a simpler and less ambitious approach in which the elec- 
tronic recombination rates calculated by Bates, Kingston, and 
MeWhirter are curvefitted as a function of electron temperature 
and number density by a simple analytic formula and the parameters 
in this curvefit are then adjusted to give recombination rates in 
agreement with experiment. Following Bowen and Park (ref. 19) , we 
have taken this curvefit to be of the form 

kf = axTe" “‘^n^ + ^ 2 (32) 

where the a^ and adjustable constc.nts. This form has the 

correct theoretical dependence on electron density in the limits 
of high and low electron densities and, with the proper choice of 
constants, can be made to fit the calculations of Bates, Kingston, 
and MeWhirter (ref. 55) at intermediate electron densities v/ithin 
about a factor of three over the entire range of conditions cov- 
ered in their calculations (i.e., for electron temperatures from 
250 to 64000°K and electron number densities from lO^ to lO^^/cm^) . 
Although a more accurate approximation could no doubt be obtained, we 
feel that the accuracy of equation (32) is probably consistent with 
the accuracy of the other approximations made in the model, and should 
be adequate to give a good prediction of the overall heat balance and 
flow parameters for the arc tunnel. Details of the electron number .nd 

excited state distributions in the flo\^7 may be less accurately given, 
however . * 

The recombination of electrons .'.n helium plasmas has been 
extensively studied experimentally (ref. 60) , The status of 
tliese experimental studies has been summarized recently by Col- 
lins, et.al, (ref. 60). Initially, many of the experiments ap- 
peared to give discordant results, apparently because of uncer- 
tainties as to the exact ion involved in the recombination process 


*Note, however, that because of the steep dependence of the recom- 
bination coefficient on electron number density and temperature, 
the erx'or in these parameters at any point in the flow resulting 
from the curvefit (32) will be much less thoii t]ie error in the 
curvefit itself (ref, 61) . 
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and in the electron temperature. However, in the more recent 
experiments, in which care has been taken to identify the ex- 
perimental parameters more exactly, a more consistent picture 
\ of the recombination process has begun to emerge, although some 

points still remain unclear. For the present study, we have 
adopted the recombination coefficient for electrons and atomic 
helium ions at high electron densities which has been recommended 
by Collins, et. al. (ref. 60) on the basis of a fit to their own 
experimental data and earlier data in which the ions involved 
appeared to be clearly identified (ref. 59). This yields for the 
high density portion of the curvefit (32), 

= 7.1 X 10~^° (Tg/300°K)“'^*^ n^ cmVsec (33) 


for the recombination of electrons and atomic ions. 


In the low electron density region, electronic recombination 
rates are controlled by the direct radiative recombination of 
electrons and positive ions, as discussed in detail, for example, 
by Bates & Dalgarno (ref. 62) . A fit to the theoretical calcula- 
tions of Burgess and Seaton (ref. 63) gives the recombination rates 
for He"^ ions in this region as follows: 

kf = 6.3 X 10"^"^ (Tg/10^ cmVsec (34a) 

for recombination into excited singlet states of the He atom, 
kf = 2.10 X 10“^^ (Tg/io"^ Ok)“0-81 cmVsec (34b) 


for recombination into excited triplet states of the He atom, and 
kf = 1.59 X 10" (Tg/lO^ cmVsec (34c) 


for direct radiative recombination into the ground state He atom. 

To determine the remaining parameters reqrired for the re- 
combination of electrons and He'*' ions in our reaction rate model 
(equations (320) and (321) in Vol. I), it is necessary to specify 
the species formed in the recomlDination reaction and the fraction 
of the reconi)ination energy going into the electron gas and into 
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radiation. Although r ecomlDination orcurs initially into highly 
excited atomic states, Bates et, al (ref, 55) have shown that the 
net change in population of these states is negligible under condi- 
tions for which the collisional-radiative model is applicable, so 
that the state of the plasma can be described (.ompletely over times 
of interest for macroscopic flow problems by giving simply the net 
rates of recombination into the ground state lielium atom and the 
two motastable . excited states He (ls2s ^s) and He (ls2s ^S) . 
Further, it has been pointed out by Bates, Bell, ind Kingston (ref. 
57) that, under optically thick conditions, practically all recom- 
bining electrons in a helium plasma will pass througli one of tlae 
metastable excited states before reaching the ground state, so that 
it is not necessary to consider recombination directly into the 
ground state. This follows because the cross sections for colli- 
sional de-excitation directly into the ground-state are much smal- 
ler than those tor de-excitation into one of the metastable states 
when the electron energy is of the order of a few electron volts 
or less, while direct radiative transitions to the ground state 
are not effective in de-exciting the gas when the plasma is opti- 
cally thick, since the emitted radiation is re-absorbed by ground 
state atoms to produce new excitation before it can leave the plas- 
ma. 


For helium plasmas at temperatures of the order of a few ev or 
less, the mean free path for the line radiation emitted by radiative 
transitions to the ground state is of the order of lOl'^/no cm at the 
line center*, where ng is the nmiber of ground state atoms per cm^. 
Thus, for the conditions of int^'rest in the NATA code, essentially all 
of this radiation will be rc-exbsorbcd before it can escape from the 
plasma. Accordingly, I ' gas v/ill be optically thick to this radiation 
and direct radiative transitions from excited atoms to the ground state 
may be neglected to a good approximation in the code. 

The situation is less clear for the continuum radiation which 
results fr^m direct radiative rccornbir.ation of free electrons into 
the ground state according to the process indicated in equation 
(34c). For this radiation the mean free path is about 1.6 x lol^/n^ 
cm**, so tiiat, for example for helium at 20,00 Qok and atmospheric 
pressure t’oe mean free path would be about 0.5 cm. This is some- 
what smaller than tiqjical nozzle dimensions, so that under these 


*This estimate assuncs Doppler broadening of the line profiles; 
this should bo valid under the conditions of interest except for 
very high excited states. 

**Calculated from (34c), using detailed balance. 
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conditions one might e^^ct that the larger part of the continuxim 
radiation (34c) would be re-absorbed and re-ionize the gas, but a 
significcint fraction, especially near the edges of the nozzle, 
would escape. As the gas ejqpanded down the nozzle, the ratio of 
the radiative mean free path to the nozzle dimensions would in- 
crease, so that eventually a point would be reached at which most 
of tiie continu\im radiaticn (34c) escaped. For higher initial 
pressures or lovwjr initial temperatures, on the other hand, the 
mean free path of the radiation would be decreased, so that under 
some conditions it might be a good approximation to treat the flow 
as optically thick to the continuum radiation (34c) over the major 
portion of its expcinsion. 

Although the helium kinetic model developed in this appendix 
does not contain any provisions for treating radiative re-absorption 
in the gas explicitly, one can allow approximately for this effect 
by adjusting the radiative recombination rate for the process (34c) 
so as to match the net radiative recombination ejqpected in the flow 
as well as possible. For this purpose it is probably most impor- 
tant to match the net recombination rate in the high temperature 
region near the nozzle entrance, since the importance of radiative 
recOTibination is expected to decrease as the flow ejqpands (see 
equations (33) and (34)), and to become negligible far downstream. 

In many cases it should be an adequate approximation to assume that 
the flow is optically thick to the reccsribination radiation (34c) 
in the nozzle entrance region, so that direct radiative recombina- 
tion to the ground state according to the process (34c) may be 
neglected in tlie calculations; however, if this approximation is 
not adequate for a particular case a better estimate r:.ay be made 
on the basis of eq. (34c) and the particular nozzle geometry. 

The reaction rate parameters given in Table XX for the recom- 
Mnation of electrons and atomic He"^ ions are derived from equa- 
tions (32) through (34) on the assumption that the gas is opti- 
cally thick to all radiation arising from transitions to the He atom 
ground state, so that essentially all recombinations will produce 
a raetastablc helium atom in cither the 2s-^S or 2s^S state. Since 
direct information as to the relative nuiribers of electrons recom- 
bining into each of the two me testable states does not appear to 
be available at present for the higher e3ectron densities, we 
have assumed that the two states will be populated in proportion 
to their statistical weights (ref. 59), i.e., 3/4 of the recom- 
binations (33) will lead to atoms in the metastable state and 
1/4 to atoms in the state. For the lower electron densities, 
t:]ic nuuiber of electrons recombining into either of the two 
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metastable states is, of course, given directly by (34a) and (34b). 
For non-optica lly-thick conditions, a term based on equation (34c) 
may also be included in the model, as discussed above, to account 
for radiative reconbinations directly into the ground state. 

In addition to the rates for collisional-radiative recombi- 
nation discussed above. Bates, Kingston, and McWhirter (ref. 55) 
have also calculated rates for the reverse process of "collisional- 
radiative ionization" frcan their model. At the higher electron 
densities, their results show that the collisional-radiative ioni- 
zation rate is given to a good approximation by applying detailed 
balance arguments based on the electron temperature to the calcu- 
lated overall collisional-radiative recombination rate, as indicated 
in Table XX; however, at the lower electrcn densities the ioni- 
zation rates calculated from the collisional-radiative model may 
fall significantly below the values predicted from these simple 
detailed balance arguments. We have not attempted to fit the 
calculated ionization rate data in the present model, however, 
since for an ionized gas flow expanding through a nozzle, ioniza- 
tion will generally be negligible compaied to recombination in the 
region where the detailed balance estimates of the ionization rate 
become inadequate. This would not be true, however, for cases 
such as the ionization of a gas behind a shock wave in which addi- 
tional ionization is being produced in an initially cold gas, so 
that the reaction parameters given in Table XX would need to be 
revised to treat such cases. 

According to the collisional-radiative model, all of the re- 
combination energy of an electron-ion pair is transferred either 
into kinetic energy of the electrons or into radiation. For the 
lower electron densities, collisional processes are unimportant 
so that the entire recombination energy of the atom, together with 
the initial kinetic energy of the electron, will be emitted as 
radiation, as indicated ror reactions 3 to 5 in Table XX. For 
higher electron densities there is a close coupling between col- 
lisional and radiative de-excitation processes so that the exact 
fraction of the recombination energy which will be emitted as 
radiation can only be determined from a complete solution of the 
collisional-radiative equations as formulated by Bates, Kingston, 
and McT"Kiirter. However, as a rough approximation, experimental 
data indicate that the total radiant emission from a helium plas- 
ma at high electron densities does not differ from the predictions 
of the low density formula by more than about a factor of two over 
the range of conditions for which radiant emission makes a signi- 
ficant contribvition to the overall energy balance of the system. 


- 197 - 




r 


T 


f 


T 


In the present model we have, therefore, used the low electron 
density formula to calculate the ladiation due to recombination 
under all conditions and have assumed that the rest of the re- 
combination energy goes into the kinetic energy of the electron 
gas, as indicated for reactions 1 and 2 in Table XX. 

In addition to reactions 1 and 2 of Table XX in which the 
collisional processes contributing to recombination are assumed 
to be with an electron as the third body. Bates and Khare (ref. 

56) have also predicted recombination rates for processes stabil- 
ized by collisions with a ground-state helium atom. We have not 
included such processes in the present model, however, since their 
existence does not appear to be verified by the experimental data 
(ref. 60). 

De-excitation of metastable atoms .- Bates and Kingston (ref, 
61) have pointed out the importance of the metastable atom popu- 
lation in determining the overall electronic recombination rate 
in a decaying helixim plasma. This effect arises because, as we 
have seen above, the net collisional-radiative recombination rate 
is a strong function of electron temperature (see equation 33) , 
and the metastables serve as an energy source for the electrons, 
raising the electron temperature and thus impeding the recombina- 
tion process. A proper treatment of the processes determining 
the metastable population in the flow is thus important if one 
wishes to obtain an accurate prediction of electronic recombina- 
tion rates in an expanding gas. 

The processes which may lead to the destruction of metastable 
atoms in a decaying helium plasma have been surveyed by Bates, 

Bell, and Kingston (ref. 57) . For the present one-dimensional flow 
model we shall neglect the loss of metastables from the flow due 
to diffusion and de-excitation at the walls. This should be a 
good approximation in the region outside the boundary layer, where 
the one-dimensional model is expected to be applicable. Further, 
the de-excitation of metastables by direct radiative transitions 
to the ground-state (two photon emission) is completely negligible 
in helium for gas densities of interest in laboratory applications 
(ref. 57). Thus, the rate at which metastable atoms are removed 
from the flow will be determined entirely by collisional processes 
in the present model. 

Approximate rate constants for several of Uie processes lead- 
ing to the destruction of metastable atoms in a helium plasma have 
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been given by Bates, Bell, and Kingston (ref. 57) in their study 
of metastable atom populations in a decaying plasma. For the con- 
ditions of interest in the present study, it appears that the most 
important metastable destruction process will ordinarily be the 
de-excitation of metastable atoms by collisions with slow elec- 
trons, the excitation energy being transferred to the kinetic 
energy of the electron. Bates, Bell, and Kingston have calcu- 
lated the reaction rate for the de-excitation of a He(3s) meta- 
stablc atom by this process in the electron temperature range 
from 250 to 4000°K, using a cross section obtained by detailed 
balance from the measured He (3s) excitation cross section of Schulz 
and Fox (ref. 64) and averaging over a Maxwellian distribution of 
electron energies. The parameters given for this reaction in 
Table XX (reaction no. 6) were obtained from a curvefit to their 
calculations. Figure 62 compares this curvefit for the reaction 
rate with the original calculations of Bates, Bell, and Kingston 
(ref. 57). 

Bates, Bell, and Kingston do not give the rate for de- 
excitation of He (Is) metastable atoms to the ground-state; how- 
ever, since the excitation cross section for the He(ls) meta- 
stable state is about 1/3 that for the He(3s) state (refs, 65, ‘ 

66), it appears, taking account of the differing multiplicities 
of the two states, that the calculated rates given by Bates et, 
al. , (ref. 57) for He (3s) de-excitation should also be approxi- 
mately applicable to He(ls), as we have assumed for reaction 7 
in Table XX. 

In addition to de-excitation to the ground state, electron 
collisions with metastables can also produce transitions be- 
tween the He (Is) and He(3s) metastable states according to the 
reaction scheme 

He(ls) + e" ;^=iHe(3s) + e“ + 0.796 ev (35) 


Pheli'3 (ref. 67) has measured a reaction rate for this process 
of 3.5 X 10“^ cm^/sec at SOO^K, corresponding to a reaction cross 
section of 3 x lO"^^ cm“. To estimate the temperature dependence 
of the reaction rate, we make use of the work of Johnson and 
Hinnov (ref. 58) , v;ho have estimated a reaction cross section of 
about 10~15 cm2 foj- process (35) at electron temperatures of 
the order of 10,000°K, based on their spectroscopic studies of 
the population distribution of the helium excited states. Tlius 


- 199 - 



FIGURE 62 - REACTION RATE FOR THE COLLISIONAL DE-EXCITATION PROCESS 



ELECTRON TEMPERATURE, Tg 
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the reaction cross section is approximately proportional to 
1/Tg and the reaction rate to as shown in Table XX. 

It may be noted that the reaction rate for process (35) is 
two to three orders of magnitude greater than the rate for de- 
excitation of metastable atoms to the ground-state by electron 
collisions, so that «s the electron density in the gas decays 
one may expect the relative populations of the He(^S) and He(^S) 
metastable states to remain in approximate thermodynamic equi- 
librium with each other at the electron temperature over a con- 
siderable range of conditions. 

When the density of metastable atoms in helium becomes com- 
parable to the electron density, a significant nximber of meta- 
stable atoms may also be removed from the gas by the Penning 
ionization process 

He (3s or ^S) + He (3$ or ^S)— >He + He'*’ + e“ , (36) 

in which two metastable atoms collide and the excitation energy 
of one of them is transferred to ionize the other. This process 
has also been considered by Bates, Bell, and Kingston (ref. 57) 
who have shown that the reaction rate is given approximately by 

kf = 6.7 X 10-10 t 1/^ cm3/sec (37) 

3 

when both metastable atom.s are in the S state, where T is the 
heavy-particle temperature in°K. Their analysis ma^ also be 
applied to collisions in which one or both of the metastahle 
atoms is in the state by using the appropriate van der Waals 
force constant for the interaction (ref. 68) and noting that the 
spin conservation factor in the analysis of Bates, et. al. (ref. 
57) is 1 instead of 4/9 when either of the colliding atoms is in 
the Is state. This procedure yields the reaction rate constants 
given for the three Penning ionization processes (36) in Table 
XX (reactions 9, 10, and 11). 

For fractional ionizations less than about 0.01 percent, 
collisions with ground-state neutral atoms may also make a sig- 
nificant contribution to the de-excitation of metastahle atoms 
in helium. In the case of the He(ls), collisionally induced 
radiative transitions to the ground-state according to the 
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He (Is) + He — > He + He + hi/ 


(38) 


appear to be the primary de-excitation mechanism at low electron 
densities, where He denotes a helium atom in the ground elec- 
tronic state. The reaction rate for the process (38) has been 
studied both experimentally and theoretically (ref. 69) and all 
determinations appear to be in reasonably satisfactory agreement 
(i.e., within a factor of about 3 or 4) . For the present model, 
we have used an approximate curvefit to the theoretically calcu- 
lated temperature dependence (ref- 69) with the value normalized 
at 300°K to the reaction rate = 6 x 10“^^ cm^/sec measured 
by Phelps (ref. 67) . 

In the case of the He(3s) metastable atom, collisionally 
induced radiative transitions to the ground state of the type 
(38) arc forbidden by spin conservation, and the metastable 
atom is removed from the gas at low electron densities primarily 
by conversion into the metastable molecular state He 2 (a^ 2 "'") 
according to the three body reaction (ref. 67) 

He + He + He(3s)-^ He + Hc 2 (a^S^) (39) 


v/itli the reaction energy presumably going primarily into trans- 
lational Cind vibrational energy of the heavy particles. The 
reaction rate for this process measured by Phelps at 300°K is 
k = 2.b X 10“34 cm^/sec, so that the process should be negli- 
gible except at quite high gas densities. The temperature de- 
pendence of the reaction rate is unknown; however, it appears 
unlikely that the rate v;ould vary greatly with temperature over 
the range of conditions for which the process (39) might be 
important, and v/e have arbitrarily assumed a 7 T dependence in 
Table XX. 

Molecular species .- For pressures greater than about 1 mm 
Hg and low temperatures, it has been observed that atomic He’^ 
ions are reipidly converted into molecular ions. Although a num- 
ber of different molecular ions have been observed (ref. 69) 
only the ground- state diatomic ion HC 2 ''' will be considered in 
the present note, since it is the ion v;hich is formed initially 
by He'*' attaclunont and appears to be the only molecular ion which 
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could be present in significant concentrations under the rela- 
tively hi 0 i temperature conditions existing in an expanding 
plasma jet. The thermochemical properties of the Hc 2 '*' ion now 
appear to be fairly well established, and are given in Table XIX. 
The He 2 '*' dissociation energy had been rather uncertain until quite 
recently, but tne latest experimental and theoretical results 
(ref. 9) now appear to strongly support a value of 2.50 ev for 
the electronic dissociation energy of the ground-state He 2 ‘*‘ 
ion. Since tne energy differences between the various electronic 
states of He 2 end *Ie 2 ‘*' are accurately kno\-m from spectroscopic 
data (refs. 9, 2.2) hhe use of the above value for the He 2 "^ dis- 
sociation energy, together with the available spectroscopic data 
(refs. 12, 70) serves to completely determine the thermochemical 
properties of the He 2 molecule and the He 2 "^ molecular ion. The 
values of the lowest stable states of He 2 and He 2 ‘^ are summar- 
ized in Table XIX. 

The principal process leading to the formation of He 2 ‘^ 
molecular ions at pressures above about 1 torr appears to be 
the three-body attachment reaction 

He + He + He”^ ^ He 2 + + He (40) 


Several independent measurements (refs. 71-73) of the reaction 
rate for this process have yielded values for the rate constant 
which agree within about a factor of two at room temperature. The 
temperature dependence of the rate constant is somewhat unclear, 
with Beatty and Patterson (ref. 71) reporting a rate constant 
which is approximately independent of temperature, while Niles 
and Robertson (ref. 72) report a T“^ dependence over the tem- 
perature range from 77°K to 449°K; however, this difference is 
pe laps not too significant in view of the rather limited tem- 
p cat.' e range for which molecular ions may be expected to be 
importa.it in the gas. For the present model, we have adopted 
the reaction rate of Beatty and Patterson (ref. 71) for the pro- 
cess (40) , 3 ,e . , 

= 1.08 X 10"^ cm^/sec 
independent of gas temperature. 

Tlie production of molecular ions by the associative ioniza- 
tion (liornbcck-Molnar) process 

lie* + He HC 2 '^' + e“ (41) 
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has also been observed for excited He* atoms in 3p electronic 
states or above (ref. 74), but this reaction docs not appear 
to be a significant source of molecular ions under the conditions 
existing in an eicpanding gas flov;, and accordingly has not been 
included in the present model. 

The primary mechanism for the recombination of He 2 "*' mole- 
cular ions again appears to be the collisional-radiative pro- 
cess of Bates, Kingston, and McWhirter (ref. 55) in which elec- 
trons are initially captured into highly excited molecular states 
and are then subsequently stabilized by collisional and radiative 
transitions to lower states. In spite of repeated experimental 
studies, the dissociative recombination process 

He 2 ‘’‘ + e'^He + He* (42) 

has never been definitely observed in helium, and the reaction 
rate for the process appears to be almost certainly much less 
than 10“8 cm^/sec (ref, 75),* 

The best data presently available for the He 2 '^ recombination 
rate appears to be that of Berlande, et. al. (ref. 77) who find 
a rate constant of the form 

kf = 5 X 10“^*^ + 2 X 10“20 Hg + 2 X 10“2^ njjg cm /sec 

(43) 

for He 2 '*' recombination at an electron temperature of 300°K. 
Measurements at higher electron temperatures (ref. 78) indi- 
cate a temperature dependence at high electron densities similar 
to that found for atomic Ile"^ ions (equation 33), so that it is 
consistent with the available data to treat the reccmbination 
of He 2 '*' ions at the higher electron densities as a collisional- 
radiative process with the rate constant 

kf - 2 X 10"^° (Tg/300°K)"^*3 Hg cmVsec , (44) 


*This interpretation of the data has been recently questioned 
by Johnson and Gerardo (ref. 76), however. 
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or about one fourth the rate constant found for atomic He"*" ions 
under similar conditions (equation 33) , 

Since infomation on the final products of the recombination 
process (44) is not available, it again appears reasonable to 
assume that the singlet and triplet molecular states are popu- 
lated according to their statistical weights. The triplet 
states then presumably cascade dov/n by^collisional and radiative 
transitions to the metastable He2(a^2]!^) state at 17.937 ev 
above the He atom ground-state (see Table XIX) , while the singlet 
states cascade down to the lowest singlet state of the He 2 mole- 
cule, namely the unstable He 2 (X ground-state, which then 

imniediately dissociates into two ground-state He atoms. Since 
all electronic states of the Ho 2 molecule have approximately the 
same equilibrium internuclcar separation re (ref. 12) , it seems 
reasonable to assume that the He 2 (X '*') ground-state is formed 
with an internuclear separation equal t§ the separation r^ =1.08 A 
of the He 2 *^ ground-state, corresponding to a potential energy of 
about 2.31 ev (ref. 79). This potential energy is then converted 
into kinetic energy of the dissociating helium atoms, while the 
remainder of the recombination energy, equal to 22.190 - 2.31 = 
19.88 ev per molecule, is converted into electronic kinetic en- 
ergy and radiation by the collisional- radiative process. As with 
the atomic recom]Dination process, we have not attempted to dis- 
tinguish betv;een the energy going into electronic kinetic energy 
and into radiation in the present note, but have simply assigned 
all of the excess reconfoination energy to the electrons in Table 
XX (reactions 15 and 16) . Although we expect this to be a rea- 
sonable approximation for cases in which the He 2 '^ recombination 
energy is important, this has not been definitely verified. 

Since the first term in equation (43) appears to be much 
too large for a simple radiative recorolDination process, it has 
been tentatively ascribed (ref. 75) to the dissociative recom- 
bination process (42) , where He* may represent either a meta- 
stable or helium atom. In Table XX we have again assumed 
that the singlet and triplet states arc populated in accordance 
with their statistical weights, and have tajeen the rate constant 
to be independent of temperature, altho\igh there is some theore- 
tical evidence to indicate that it may actually be an increasing 
function of gas tcmperatuLc (ref. 80) . Although the rate constant 
for the process (42) is poorly kno\sm, it represents a minor 

correction to the calculated net reaction rate (43) under most 
conditions of interest, and should thus not contribute apprec- 
iably to the overall uncertainty in the calculated gas conditions. 
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The final te\m in the observed He 2 '*' recombination rate (43), 
which is proportional to the gas density, presumably represents 
the effect of stabilization by collisions with ground-state He 
atoms which was studied by Bates and Khare (ref. 56) . Since 
this term becomes important only when the fraction of ionization 
is of the order of lO”^ or less, we have not included it in the 
reaction rate model given in Table XX, 

To complete the present reaction rate model for helium, 
approximate rates for the destruction of the raetastable He^ 
molecule are included as the final two roagtions in Table XX. 
Since the experimental studies of Collins (ref. 81) and Phelps 
(refs. 67, 82) indicate that the processes 

He2(^2) + e”^2He + e~ (45) 

and 

+ He 2 (^Zl) ^ 3He + He"^ + e" (46) 


for the destruction of metastable He 2 (^ 2 ) molecules have the 
same rates at room temperature, within the experimental error, as 
do the corresponding processes for the He(3s) metastable atom, 
we have, for lack of any better information, simply used the 
reaction rates given previously for the He(3s) reactions (reac- 
tions 6 and 9 in Table XX) for the processes (45) and (46) as 
well. The destruction of metastable molecules by collisions with 
ground-state atoms appears to be negligible and is hence not in- 
cluded in the present model; Phelps' data (ref. 67) indicate a 
reaction rate for this process at least two orders of magnitude 
smaller than for the corresponding process (39) for metastable 
He(^S) atoms. 

Concluding- remarks - Since we have introduced a number of 
simplifying ap]proximations in constructing the reaction rate 
model for helium presented in Tables XIX and XX, and since, 
moreover, several of the important reaction rates for helium 
are still rather uncertain, especially at the higher tempera- 
tures, it would now be desirable to verify the model by compar- 
ing its iDrodictions with experimental data over as wide a range 
of conditions as possible, and, if neeessary, adjust the rate 
constants to obtain satisfactory agreement with experiment. This 
has not been possible within the scope of the present study, how- 
ever-; and accordingly the re^lction rates given in Table XX should 
be regarded as provisional until sudr time ns a more complete 
verification of tlie model can be obtained. 


- 206 - 


A, 2 Argon Model 


I 


( 


The nonequilibrium argon model used in NATA is basically 
similar to the helium model described above, but with the para- 
meters adjusted and a few minor modifications made to account 
for the difference in physical properties between helium and 
argon. Thus, much of the discvission given above for helivun is 
also applicable to argon, and only the differences between the 
two gases need be noted here. 

Although argon has been used extensively as a test gas in 
laboratory studies for various aerodynamic and arc tvinnel appli- 
cations, the reaction mechanisms in recombining argon have appar- 
ently not been studied in as much detail as they have for helium, 
and in consequence, as will be indicated in more detail in the 
discussion below, several of the important parameters in the 
argon recombination model appear to be significantly vmcertain 
at the present time. Accordingly, the errors in the nonequili- 
brivim model calculations for argon may be expected to be larger 
than for helium. 

The species and parameter values for the nonequilibrium 
argon model used in NATA are given in Tables XXI and XXII. 

Elastic collisions .- The electron energy loss due to elas- 
tic collisions in argon is again calculated from equations (25) 
to (30) with appropriate values of the momentum transfer cross 
sections for argon being used in equation (30) . As in the case 
of helium, the approximate Couloiiib cross section (31) is used 
in equation (30) for all electron-ion collisions. 

Data on the momentvim transfer cross section between electrons 
and ground-state argon atoms have been given by Frost and Phelps 
(ref, 83) and by Golden (ref, 84). For the present model, we 
have nvimerically integrated the data of Frost and Phelps over 
electron energy as indicated in equation (30) to obtain the 
Maxwell-averaged electron-argon atom cross section shown 

in figure 63 as a function of temperature. Use of the data of 
Golden in tnis computation would have given a noticeably lower 
cross section in the neigliborhood of the Ramsauer minimum at 
T = 1300°K; however, since the total cross section in this reg- 
ion is so small for either calculation, the effect of such a 
change on the overall electron energy balance for the gas would 
be negligible. 
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FIGURE 63 - MAXWELL AVERAGED MOMENTUM TRANSFER CROSS SECTION FOR 

COLLISIONS BETWEEN ELECTRONS AND GROUND STATE ARGON*”^OMS 



ELECTRON TEMPERATURE . T^ (^'K) 


As in the case of helium, wo have used the calculated 
electron-ground state argon atom cross section shown in figure 63 
for all elastic collisions between electrons and neutral argon 
species. As an indication of the error involved, figure 64 shows 
the Maxwell-averaged cross section for collisions between elec- 
trons and metastable argon atoms which we have estimated from 
the available cross section data (refs, 85 and 86), One sees 
that the estimated metastable cross section is about two to 
three orders of magnitude greater than the ground state cross 
section over the temperature range of interest. Thus, the ap- 
proximation employed in the code of using the ground-state cross 
sections for all neutral species should give an adequate repre- 
sentation cf the elastic energy losses as long as the concentra- 
tion ci metastable species in the gas remains <0,1 percent. As 
noted previously i our discussion of the helium model, the 
elastic energy losses themselves become negligible at higher met- 
astable concentrations, so that the error in the calculated elas- 
tic energy losses at high metastable concentrations ( ^ 1 percent) 
should not significantly affect the overall accuracy of the miodel 
predictions in this region, 

Collisional- radiative recombinatj ons .- Data on electronic 
recombination rates in argon have been reviewed recently by 
Diberman, et. al. (ref. 87) . Although the experimental uncer- 
tainty is larger than for helium, the available data for argon 
again appear to be generally consistent with the collisional- 
radiative recombination mechanism suggested by Bates, Kingston, 
and MeWhirter (ref. 55) . Accordingly, we have followed the ap- 
proximate model of collisional-ra'^iative recombination discussed 
above for helium in the present treatment of argon recombination 
also. 


As with helium, we have attempted to represent the observed 
recombination rate data for argon by a curvefit of the form (32) . 

For the higher electron densities, the observed recombination 
rates for argon (refs. 88, 89) are found to agree with the hel- 
ium data within the experimental scatter, so that the high den- 
sity portion of the reaction rate curvefit for helium (equation 
33) may also be used for argon. 

At low electron densities, radiative reconnbination becomes 
dominant and recombination rates may be determined from available data 
on the argon continuum radiation. From these data (refs. 90-92) 
one finds that the total rate for radiative recombination into 




the excited states of the argon atom agrees with the correspond- 
ing rate for helium within the experimental scatter, so that the 
net reaction rate for recombination into an excited 
at low electron densities becomes 

kf = 2.73 X 10“ 13 (Tg/104 ok)- 0.81 cm^/sec 

For radiative recombination directly into the argon 
the data of Samson (ref. 93) give the recombination 

ICf = 1.00 X 10“^^ (Tq/ 10^ cm3/sec 

which is somewhat lower than the corresponding rate 
helium. 

The data cited above indicate that the rates for electronic 
recatibination into excited atomic states are ccmiparable in helium 
and argon at both high cind low electron densities.' At intermed- 
iate electron densities, however, the data of Chen (ref, 94) 
give a reccmbination rate in argon which is several times higher 
than the corresponding helium rate and than the collisional- 
radiative predictions, and about an order of magnitude higher than 
the rate predicted from the simple curvefit (32) on the basis of 
the high and low electron density data,* The reason for this dis- 
crepancy is not clear, Chen suggests that the obseirved differ- 
ences between helium and argon in his work may be due to differ- 
ences in the electronic excitation and de--excitation cross sec- 
tions for excited states in the two gases; however, this c. lan- 
ation does not appear to be consistent with the close agreement 
between helium and argon recombination rates at high electron 
densities which has been observed in other studies (refs, 88,89), 
In view of this apparent inconsistency and the lack of other ex- 
perimental data to support the difference between helium and argon 
recombination rates found by Chen, the present model uses a re- 
combination rate based only on the high and low electron density 
data given by equations (33), (47) and (48), and ignores the 
data of Chen at intermediate densities. These data should be 
borne in mind, however, as a possible indication of significant 
uncertainty in the predicted recombination rates for argon. 

Because of the closed electronic p-shell in the heavier rare 
gas atoms, the relationship among the low-lying excited states 


argon atem 

( 47 ) 

ground state 
rate 

( 48 ) 

(34c) for 


*Ren^ember that the curvefit (32) gives recombination rates sev- 
eral times smaller than tlie correct collisional-radiative model 
at intermediate electron densities. 
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differs from that found in helium with the result that the decay 
of excited states during collisional-radiative recombination in 
argon follows a somewhat different pathway than that described 
previously for helium. Table XXIIl shows the four lowest lying 
excited states of the argon atom, together with the three com- 
monest designations by which tlxcy are referred to in the litera- 
ture. All higher excited states of the atom can decay rapidly 
into one of these four low-lying excited states, so that only the 
populations of these four states need be followed in the colli- 
sional-radiative model. Of the four states, tiie and ^Pq are 
metastable v/hile the and ^Pj^ can decay to the ground state 
by emission of resonance radiation. The decay is slow, however, 
because of the trapping of resonance radiation in the gas, and 
the two resonance states ^P^^ and ^P]^ thus behave somewhat like 
true metastable states under many conditions of interest. 

Within the accuracy of the present modal it seemed unneces- 
ary to distinguish between all four of the low-lying excited states 
in Table XXIII, and accordingly we have grouped the tv/o metastable 
states ^P 2 and into a single metastable state Ar* (ra) and the 
two resonant states, ^P^^ and ^Pj^, into a single resonant state 
Ar* (r) , as indicated in Table XXI. Assuming that these states 
are populated in proportion to their statistical weights by the 
recombination reactions (33) and (47) then leads to the reaction 
rates for electronic recombination into excited states given for 
reactions 1 through 4 in Table XXII. As in the case of helium, 
we find that an adequate fit to the total visible and infrared 
emission from the gas under all conditions (refs. 90-92) can be 
obtained by assviming that all of the recombination energy goes 
into the electron gas for the three-body recombinations 1 and 2 
while all of the energy* goes into radiation for the two-body re- 
combinations 3 and 4. 

For radiative recombination directly into the argon atom 
ground state, reabsorption of the emitted reconbination radiation 
by the gas can be effective in reducing the net reconbination rate 
even at moderately low gas densities. An estimate of this reduc- 
tion can be made by dividing the low density recoirbination rate 


*It may be noted that in computing the electron thermal energy for 
the reactions in Table XXII, we have taken account of the fact that 
the reaction rate depends on electron energy, so that the mean energy 
of: the electrons participating in the reaction is not the same as 
tthe mean thermal energy 3/2 k Tg of all the electrons in the gas. 

! .Under most conditions of interest, hov;ever, this difference will not 

bo significant for the final results of the calculation. 
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TABLE XXIII 


LOW- 

•LYING EXCITED 

STATES OF THE 

fVRGON ATOM 

LS 

Desianation 

Paschen 

Desianation 

j 1-Coupling 
Desianation 

Energy 

(ev) 

Statistical 

Weight 

- a 

4s ^p 2 

Iss 

4s 13 / 21 ° 

11.545 

5 

4s 

IS 4 

4s[3/2]i° 

11.620 

3 

4s \ 

IS3 

4s‘[V2](j° 

11.720 

1 

1 

4s Pi 

IS2 

4s'lV2]j^° 

11.825 

3 


(48) by the optical depth of the gas for the reccatibination radi- 
ation 


'r = 3.4 X 10” n, R, 

Ar 


(49) 


Where T is the optical depth, n^^j- is the number density of ground 
state argon atoms in the gas in cra"^, R is the channel radius in 
cm, and 3.4 x 10”^^ cm^ is the absorption cross section of ground 
state argon atoms for the recombination radiation, as measured by 
Samson (ref. 93). As indicated in Table XXII (reaction 5), this 
reduction must be applied to the rate constant (48) whenever the 
optical depth (49) of the gas becomes greater than one. 

Decay of excited atoms .- 'Xl^e de-excitation of lo^ 7 - lying atomic 
states in a decaying argon plasma "appears to be due primarily to 
collisions with electrons and ground-state atoms, and, in the case 
of the resonant states, the emission of radiation. For ionization 
fractions greater than about 10 ”^ to 10 " 5 ^ electron collisions and 
resonance radiation will be tlie most important de-excitation mech- 
anisms, while atomic collisions become important at lower frac- 
tions of ionization. 

De-excitation rates for electron collisions with excited 
argon atoms can be obtained from experimental data on the inverse 
process of excitation of ground-state argon atoms by electron 


213 - 


\ 


impact. Rates for the latter process have been obtained both 
from studies of the ionization rate behind argon shocks (refs, 

95-97) and frcan electron beam measurements of the excitation 
cross sections versus electron energy (refs. 98-102). 

In all of the shock tube studies, it has been found that 
the rate determining step in the ionization process behind the 
shock front is the initial excitation of the grovind state argon 
atom to a lov/-lying excited state, so that the excitation cross section 
can be detemined directly from the measured ionization rates. 
Experimentally, the ionization is observed uo proceed ir two 
stages; first an initial induction phase in v/hich the excitation 
is produced primarily by collisions with ground state atomp, 
followed by -i second, much more rapid stage, in which the elec- 
tron density has become sufficiently high for electronic col- 
lisions to contribute significantly to the observed excitation 
rates. Analysis of the ionization rate data for the second stage 
indicates that the measured rates are consistent with a linear 
dependence of the electronic excitation cross section on excess 
electron energy above the excitation threshold (~11,5 ev) , with 
a slope which varies in the different experiments (refs. 95-97) 
over a range from about 5 x 10~18 to 7 x 10 “ 18 cm^/ev. Since 
the e 2 q)eriments are not sensitive enough to determine exactly 
which state is being excited, this measured excitation cross 
section should probably be regarded as a sum over the four low- 
lying states of the argon atom indicated in Table XXIII, 

Total cross sections for excitation of the two argon meta- 
stable states ^P 2 and ^Pq by electronic collisions have been 
measured in electron bean studies over the energy range from 
threshold to -^200 ev (refs, 98-101). The data are in general 
agreement with the linear cross section dependence on electron 
energy assumed in the analysis of the shock tube data; hov;ever, 
considerable detailed structure in the cross section energy de- 
pendence is evident near threshold (refs. 98,99) so that the 
cross section slope deduced from the shock tube experiments 
w 'uld be expected to vary somewhat with electron temperature. 

Uv-rng Pichanick and Simpson's relative cross section measure- 
ments near threshold (ref. 98) normalized to the absolute cross 
ss^-tion values given by Borst (ref. 100), we find a mean cross 
section slope varying between about 4 x 10 “ ^8 cm^/ev and 8 x 10 ” ^8 
cm^/ev for the range of temperatures below about 40,000 °K, v/ith 
the value being ~ 8 x 10 “ ^8 cra^/ev for T < 1000°K, In the range 
of temperatures T ^ SOOO^K, it appears that a linear cross section 
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with the slope 


7 X 10*“^® cm^/ev (50) 

found by Petschek and Byron (ref. 95) should give a very satis- 
factory fit to the available electron beam data. 

The cross section for excitation of the resonant and 
IPl argon states has been measured by McConkey and Donaldson 
(ref. 102) in an electron beam aparatus over the energy range 
from thre.shold to about 2000 ev. . In the threshold region, their 
results indicate a linear dependence of the total excitation 
cross section on electron energy with a slope 

^ 1 X lO-^® cmVev. (51) 

Thus, the total excitation cross section for the argon resonant 
states is considerably smaller th^m for the metastable states 
at thermal energies, and the total excitation cross section for 
all lov/- lying states derived from the electron beam measurements 
is in excellent agreement with the shock tube results. 

The reaction rates for de-excitation of the low-lying meta- 
stable and resonant excited states of atomic argon by electron 
collisions V'/hich are used in the present argon kinetic model 
are derived from the corresponding excitation cross sections 
(50) and (51) by detailed balance. Tliese rates are listed as 
reactions 6 and 7 in Table XXII. 

In addition to causing de-excitation of the low-lying excited 
states, electronic collisions can also result in transition be- 
tween the metastable and resonant excited states according to the 
scheme 

Ar*(r) + -e" Ar* (m) + e“+ 0,07 ev (52) 

Although we have been unable to find any direct data, the rate 
for the process (52) is expected to be large (ref, 103) , in 
analogy to the corresponding processes in helium (ref. 67) and 
neon (ref, 104) . An approximate upper bound for the possible 


dependence 
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reaction rate is provided by the total cross section measure- 
ments of Celotta, et, al, (ref. 85) for the inverse process of 
electron-metastable collisions. 

For the present model we have taken the reaction rate for 
the process (52) equal to the corresponding rate estimated for 
neon at 300°K by Phelps (ref. 104), and have asstuned a 
temperature dependence, as in the helium model, to obtain the 
rate constant 

-Js 

kf = 10“® (— — ) cm^/sec (53) 

300°K 

for the process (52) . This rate lies about an order of magnitude 
below the upper bound provided by the total (elastic plus inelas- 
tic) cross section measurements of Celotta, et. al. 

It should be noted that the process (52) may play an impor- 
tant role in the decay of metastable argon atoms at low electron 
densities, since it can cause trcinsitions from the relatively 
long-lived metastable excited states to the much shorter lived 
resonant states. Thus, the uncertainty in the rate consteint 
(53) could result in significant errors in the model predictions 
of metastable decay rates for ar* argon plasma under some condi- 
tions. 

In helium, all of the excited atomic states which are cap- 
able of direct radiative transitions to the ground state can also 
decay rapidly into one of the two raetastable excited states, so 
that most of the excited atoms pass through one of the metastable 
states during the de-excitation process and only a small fraction 
are de-excited by direct radiative transitions to the ground 
state. Thus, the latter process could be neglected to a good 
approximation in the formulation of the helium kinetic model dis- 
cussed in the preceding section. For argon, on the other hand, 
the low-lying and resonant excited states cannot decay 
radiatively to any lower excited state, so that direct radiative 
transitions to the ground state become the dominant decay mech- 
anism for these states in many situations, and hence must be in- 
cluded in the kinetic model. 

Since the mean free path of resonance radiation in a gas is 
typically very much shorter than usual laboratory apparatus, 
dimensions, most of the photons emitted by the radiative decay 
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of excited atoms to the groxand state are reabsorbed by other 
atoms in the gas to excite them to the resonant state, and the 
net rate of decay of the resonant state population density in 
the gas is much slower than would be predicted solely on the 
basis of the spontaneous emission coefficient for the state. 

A model to predict the net rate of decay of the population den- 
sity under these conditions as a function of apparatus geometry 
and the absorption coefficient of the resonance line has been 
developed by Holstein (refs. 105,106) and has been verified 
experimentally for a number of gases (refs. 104-108) , including 
both argon (ref. 107) and neon (ref. 104) . 

Under the usual laboratory conditions, it is generally a 
good approximation to assume that the shape of the resonant 
radiation lines is determined by pressure broadening according 
to the dipole-dipole model of Fursov and Vlasov (see refs. 109, 
110). With this assumption, the net decay rate for the tran- 
sition to the ground state predicted by Holstein's model becomes 
independent of gas pressure and, for a cylindrical gas volume, 
reduces to the simple form 

kp = 0.205 A„j(-^)^ (54) 


where Ajjj is the probability for spontaneous emission of a reso- 
nant photon from the excited state per unit time, X is the 
wave-length at the center of the resonant line, and R is the 
cylinder radius. When the appropriate constants for the argon 
resonant states (ref. Ill) are inserted into equation (54) , one 
finds the reaction rates 

kf = a-.2#-12i (35a' 

V R 

for the radiative decay of the Ar (^Pj^) state, and 

kf = sec”^ (55b) 

>/r 

for the Ar (^Pj^) state, where R is the cylinder radius in centi- 
meters. Equation (55a) has been verified experimentally by 
Ellis and Twiddy (ref. 107). Since we expect the ^Pj^ state to 
be more highly populated than the ^Pj^ state under most conditions 
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of interest, equation (55a) has been used for the radiative 
decay rate of tiie Ar resonant state in the present kinetic 
model (reaction 9 of Table XXII) . 

The reaction rates (55) should be valid in argon for noz- 
zle radii R^O.l cm and ground state argon atom number densi- 
ties n^j.^ lOi^ cm“3. For lower nxiiriber densities than this, 
pressure-broadening becomes small compared to the natural line 
width of the resonance lines, and the effective rate constant 
begins to increase with decreasing number density. Eventually, 
at very low number densities, reabsorption of the resonance 
radiation of course becomes negligible, and the effective rate 
constcint approaches the spontaneous emission coefficient kf=Ajj^, 
Comparing this limit with equation (54) , one sees that the net 
radiative decay rate for the resonant argon excited states is 
reduced by about three to four orders of magnitude by the trap- 
ping of resonance radiation under typical laboratory conditions. 

For loii 7 fractional ionizations, the de-excitation of excited 
atoms by collisions with neutral atoms can become significant. 
Experimental data on the reaction rates are available both for 
the direct de-excitation process (refs. 107,112,113) and for 
the inverse excitation process (refs. 95-97, 114-117), 


Several studies of the decay of excited state population den- 
sities in low temperature argon gas (refs. 107,112,113) have 
indicated that both two- and three-body de-excitation processes 
are significant, and have given reasonably consistent values 
(within about a factor of two) for the de-excitation rates. The 
most detailed study was that of Ellis and Twiddy (ref. 107), who 
gave two- and three-body reaction rates at 300°K of 


lc 2 = (1 i 0.3) X 10“^^ cm^/sec 


IC 3 = (1.7 i 0.2) X 10“^^ cm^/sec 


(56) 


for de-excitation of the Ar (^P-) metastable state at 11.55 ev 
and 


Tc 2 = (5.7 i 0.7) X 10~15 cmVsec 
= (1.14 i 0.15) X 10“^^ cm^/sec 


(57) 


for de-excitation of the Ar (^Pq) metastable ,:ate at 11.72 ev. 
Futch and Grant (ref. 113) measured the Ar (^P 2 ) de-excitation 
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rates at 77°K and 300°K and their data yield a temperature de- 
pendence over this range of k2'^T2/3 fpr the 2-body de-excitation 
process (56a) and ]C3~T“*^*56 for the three-body process (56b). 

In their work, Ellis and Twiddy attribute the two-body de- 
excitation rate (56a) to a two-step process in which the meta- 
stable Ar (2P2) atcsn is first excited to the resonant Ar (^Pi) 
state by collision with a ground state argon atom, and then de- 
cays radiatively to the ground state. No evidence is given to 
support this assignment, however, and it appears to be inconsis- 
tent with the temperature dependence observed for the two-body 
de-excitation process at 77°K by Futch and Grant, which suggests 
that the direct de-excitation process 

Ar + Ar (^P 2 ) — >2 Ar + 11.55 ev (58) 

is responsible for the observed two-body de-excitation rate. 

Further evidence for the reaction (58) is provided by sever- 
al shock tube studies in which excitation cross sections for 
collisions between ground state argon atoms were derived from 
the initial ionization rates behind an argon shock (refs. 95-97, 
114-117). In all cases, it was found that the experimentally 
observed ionization rates for the initial ionization stage dir- 
ectly behind the shock front could be accounted for by an atom- 
atom excitation cross section which was (approximately) a lin- 
ear function of the excess atom energy above the excitation 
threshold. The slope of the derived cross section energy depen- 
dence varied somewhat between the different ejqieriments ; however, 
all results were consistent within about an order of magnitude, 
with the reported values of the cross section varying from about 
2.5 X 10"20 cm2/ev to 2.5 x 10 " cm^/ev. Using these data to 
derive a cross section for the inverse de-excitation reaction 
(58) by detailed balance, one obtains a two-body de-excitation 
rate, k 2 , which is proportional to /t and has the value k2 -10-15 
to 10-1^ cmVsec at 300°K, in good agreement with the measured 
low temperature de-excitation rates. Thus, the assignment of 
the observed two-body de-excitation rate (56) to the reaction (58) 
is verified and an approximate /t dependence of the de-excitation 
rate over the temperature range from T = 77*^K to — 10,000®K is 
substantiated . 
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The three-body de-excitation process in equations (56) and 
(57) is attributed to the two-step reaction process 


Ar + Ar + Ar* — ^Ar 2 * + Ar 

AT 2 * - — ► Ar + Ar + hv 


(59) 


in which the excited argon atom combines with a ground state atom 
to fo'^m an excited At 2 * molecule, which then immediately radiates 
to form an unstable ground-state AT 2 molecule which dissociates into 
two ground state atoms, Tlie radiation emitted during the process (59) 
is observed to peak at a wavelength of about 1265 8 (ref, 118), cor- 
responding to a mean energy of -^9,8 ev for the emitted photon, with 
the remaining '-'1,8 ev of the excitation energy of course going into 
kinetic energy of the three argon atoms participating in the reaction. 


For the present kinetic model, the two- and three-body de- 
excitation rates for the metastable Ar*.(m) atom given in reac- 
tions 10 and 11 of Table XXII have been set equal to the Ar (^^ 2 ^ 
de-excitation rates from equation (56), since we ej^ect that most of 
the metastable atoms will be in the 3p^ state under the conditions 
encountered in an expanding argon plasma jet, because of this 
state's higher statistical weight and lower excitation energy. As 
noted above, the T'^ temperature dependence of the two-body de- 
excitation process appears to be well established, while we have 
provisionally assumed the temperature dependence observed 

by Futch and Grant at lov/ temperatures for the three-body pro- 
cess. Since collisional de-excitation rates have not been ob- 
served for the resonant excited states ^Pp and ^Pp, we have also 
used the ^P 2 de-excitation rates for the resonant excited states 
in reactions 12 and 13 of Teible XXII. This approximation should 
have a negligible effect on the computed population densities for 
the model under the conditions of interest in an expanding plasma 
jet, since the radiative decay rate (55) for the resonant excited 
states is much greater than the collisional de-excitation rate 
(56), except at very high gas densities. 

Molecular species .- The only molecular species included in 
the argon kinetic model is the ground-state molecular ion Ar 2 "^. 

In contrast to helium, there appear to be no metastable molecular 
states in argon (ref. 121) so that any excited molecular states 
formed will decay very rapidly to the ground state, in times of 
the order of 0.1 ^sec or less (ref. 118) and need not be consid- 
ered in the kinetic model. The Ar 2 ground state molecule itself 
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is very weakly bound (dissociation energy Dq^'O.OI ev (ref, 123)), 
and thus should be rapidly dissociated into ground-state atoms un- 
der the conditions e:;isting in an expanding argon plasma jet. 

The thermochemical properties of the argon molecular ion 
Ar 2 ’*‘ are still very uncertain at the present time; the values 
used in the present model have been taken from the recent equi- 
librium drift tube measurements of Teng and Conway (ref. 13) and 
are shown in Table XXI. The Ar 2 + dissociation energy Dq — 1.14 
ev found by Teng and Conway is in reasonable agreement v/ith prev- 
ious experimental (refs. 119, 120) and theoretical (ref. 121) 
estimates obtained by less direct means. The vibrational fre- 
quency cUq ~ 80 cm“l derived by Teng and Conway from the thermo- 
chemical data, however, is about a factor of five lower than 
the value estimated by O'Malley (ref. 122) on the basis of an 
interpretation of observed dissociative recombination rates in 
argon. We have accepted the value of Teng and Conway in the 
present work since it appears to be more directly determined 
from the experimental data; however, the uncertainty could be 
large . 

As in the case of helium, the production of molecular ions 
in a decaying argon plasma at pressures near 1 torr or higher 
appears to be due primarily to the three-body conversion process 

Ar+ + 2Ar Ar 2 ’’' + Ar (60) 

A number of measurements of the reaction rate for the process 
(60) at room temperature have given results in reasonably sat- 
isfactory agreement (refs. 124,125), with the more precise det- 
erminations yielding a value 

k£ (2 ± 1) X 10”^^ cn^/sec (61) 

for the reaction rate at 300°K. The temperature dependence of 
the reaction rate (61) has apparently not been accurately det- 
ermined; however, available theoretical and experimental results 
(refs. 124,126,127) indicate tha . the rate should probably de- 
crease somewhat with increasing temperature. As a reasonable 
compromise among the various estimates, we have assumed a 
temperature dependence for the rate constant (61) in the isresent 
analysis, as predicted by the simple Thomson theory (ref. 124). 
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The drift tube measurements of Liu and Conway (ref. 124) 
indicate that the rate constant (61) applies only to reactions 
of the ground-state Ar"^ (^^3/2^ ion, and that the excited 
fine-structure component of the Ar"*" ion (at 0,18 ev above the 
grovmd-state) does not participate in the reaction (60) under 
the essentially electron-free conditions of their experiment. 

For the experiments carried out in decaying plasmas (n^ ~ loH 
cm“^) (refs. 125, 128), however, the rate of conversion between 
the tv\70 Ar^ fine-structure states was apparently sufficiently 
rapid that no distinction between the ’-.wo states was observed 
in the data, and the reaction rate (61) applies to the total 
atomic Ar"^ population density in the plasma. Since the latter 
experiments appear to be closer to the conditions to be expected 
in an argon plasma jet, we have grouped the two Ar+ fine-structure 
states together into a single species for the present kinetic 
model (see Table XXI) and have used the experimental reaction 
rate (61) for the combined species in obtaining the reaction 
rate data given for reaction 14 in Table XXII, 


Tlie associative ionization process 

Ar + Ar* ^ Ar 2 '^ + e~ (62) 

is known to be rapid in argon for excited atoms having ex- 
citation energies above a threshold of ~ 14.7 ev, with the rate 
constant being kf 10”^ cm^/sec near room temperature (refs. 
119,129). Under conditions of collisional-radiative recombina- 
tion, however, the population of excited Ar* atoms above the 
14.7 ev threshold is very small compared to the ion population, 
so that the associative ionization reaction (62) should not con- 
tribute appreciably to the molecular ion production rate in an 
expanding argon plasma at the relatively high pressures (p > 10“^ 
atm) of interest here. Tlius, a detailed treatment of associative 
ioni tion harr not been included in the present kinetic model. 


In contrast to the situation in helium, the recombination 
of molecular ions in argon under most conditions occurs primar- 
ily through the dissociative recombination reaction (ref. 75) 

Ar 2 '^ + e” Ar + Ar* , (63) 

and this process appears to be sufficiently rapid to make an 
important contribution to the overall recoitibination rates in an 
argon plasma even at the higher temperatures, where the mole- 
cular ion Ar 2 '*' represents a relatively minor constituent of the 


222 - 


plasma. Dissociative recombination rates for argon have been 
measured in a number of independent experiments, both at room 
temperature and at elevated temperatures (refs. 75,\22,127 130). 
The low temperature measurements (ref. 75) give a value k£--7xl0‘ 
cm^/sec for the reaction rate of the dissociative recombination 
process (63) at room temperature. The measurements at elevated 
temperature indicate that the reaction rate for dissociative re- 
combination depends on both the electron temperature T^ and the 
gas tempera urc T. A good fit to the data of Mehr and Biondi 
(ref. 130) and Cunningham and Hobson (ref. 122) has been given 
by O'Malley (ref. 122) in the form 

kf = 9.6 X 10-7 (Te/300°K)”°*^7 (p .. g-630OiVTj 

(64) 

O'Malley's fit (64) does not appear to be consistent with the 
high temperature data of Chen (ref. 127) ; however, the validity 
of the latter data has been questioned (ref. 75). 

In his original derivation of equation (64), O'Malley in- 
terpreted the observed gas temperature dependence of the dissoci- 
ative recombination coefficient in terms of a model in which 
only the lowest vibrational state of the Ar 2 '*' molecule partici- 
pates significantly in the dissociative recombination pro'-'. 

(63) and, on the basis of this model, derived a value of 
co^ ~ 630°K ~ 440 cm”^ for the Ar 2 *^ vibrational frequency. 

Thio interpretation of the data appears to be inconsistent with 
the value of the vibrational frequency 80 cra“^ recently 

derived by Teng emd Conv/ay (ref. 13) from thermodynamic data; 
however, other interpretations of the observed temperature depen- 
dence (64) which do not lead to this inconsistency appear to be 
possible. Thus, for the purposes of the present analysis, it is 
unnecessary to inquire into the validity of O'Malley's model in 
detail and we may regard equation (64) simply as a curvefit to 
the available experimental data on the dissociative rccoii:ibination 
rate in argon. 

The exact state of the excited Ar* atoms prodv^ccd in the 
dissociative recombination reaction (63) is not known at present; 
hov/ever, it apjocars likely (ref. 75) to be one or more of the 
many excited atomic states lying in the range between about 14 ev 
and the molecular ion ground-state at 14.61 ev (see ref. 6). The 
excited atom produced in the original reaction (63) will then cas' 
cade do\'m rapidly by the collisional-radiativo process to one of 



the four low-lying excited atomic states listed in Table XXIII, 
with the excess excitation energy being lost either by radiation 
or tra;->sfer to the electron gas. In setting up the reaction 
scliemo given in Table XXII (reactions 15 and 16) we have as- 
sumed ratlier arbitrily that electronic collisions will be tlie 
dcMninant de-excitation mechanism under the conditions existing 
in ein argon plasma jet, so that the entire excess excitation 
energy is transferred to the electron gas. Lacking any de- 
tailed knowledge of the final state of the reaction prodi'cts, 
we have simply assumed in reactions 15 and 16 that the resonant 
and metastable excited states will be populated in accordance 
with their statistical v/eights. Further, the rate of the inverse 
associative ionization process is calculated from detailed bal- 
ance based on the electron temperature. This is consi stent v.*ith 
t’-if assumption that de-excitation of excited atoms occurs pri- 
marily by electronic collisions, since under this assumption the 
population of the excited Ar* atoms which can initiate the assoc- 
iative ionization process will be controlled primarily by the 
electron temperature; however, the approximation is evidently 
very crude and serves only to give an indication of the condi- 
tions under which associative ionization is likely to be impor- 
tant. To obtain a more quantitative prediction of associative 
ionization rates, it would be necessary to couple the associa- 
tive ionization-dissociative recombination model, as suggested 
by Biberrcan, et al. (ref, 87) , in order to obtain a prediction 
of the population densities in the various excited atomic states. 
Such a detailed model could perhaps explain some of the apparent 
anomalies which hav been observed in the measured argon recom- 
bination rates (refs. 94,127) but would be far beyond the scope 
of the present study. 

Because of the very large dissociative recombination rate 
for Ar 2 ’*‘ molecular ions, the three-body collisional- radiative 
recombination process has apparently never been observed experi- 
mentally for Ar 2 ’’' ions. Nevertheless, this process would be ex- 
pected to be the predominant recombination mechanism for AT 2 '*’ 
ions at sufficiently high electron densities and has been in- 
cluded as the final reaction process in Table XXII, using the 
corresponding KC 2 '’' reaction rate from equation (44) . Although 
the final state of the reaction products is not knovm, we have 
assimied in Table XXII that the recoiriloining electron will cas- 
cade down through successive Rydberg states of the molecule to 
reach the unstable Ar 2 ground state, v;hich will then dissociate 
into two ground state argon atoms. As in the case of helium, 
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it should be a good approximation, at the high electron densi- 
ties for which the three— body recombination process is important, 
to assign all of the recombination energy to the electron gas. 


- 225 - 


APPENDIX B 


DIAGNOSTIC MESSAGES 


This appendix lists all of the diagnostic messages produced 
by NATA when error conditions are detected. For each message 
the following information is given: 

(1) The name of the subroutine in which the message 
is produced. 

(2) A brief description of any dumps or additional 
messages which follow the given message. 

(3) A description of the error condition indicated. 

(4) A summary of the subsequent action taken by the 
code. 

The messages are listed in alphanumeric order for the entire 
code. Lovrer case letters appearing within the messages indicate 
numerical values. 

Most of these diagnostics occur only very rarely. Many of 
them have never been observed at all with the present version 
of NATA, and are included in the program only to allow for un- 
usual input errors or for programming errors which might occur 
during future modifications. 

The Fortran variables appearing in the various dumps are 
to be defined in the NATA Programmer's Manual, Volume III of 
this final report. 

BACKSTEPPING OF PERTURBATION SOLUTION TERMINATED AFTER n STEPS. 

DIAGNOSTIC DATA FOLLOW. 

Produced in subroutine N0NEQ. 

Followed by a dump with the namelist name NEQDMP, followed 
in turn by "ERROR EXIT NO. 5 FROM N0NEQ." 

This message is produced in the nonequilibrium solution by 
the perturbation method when ISX.I is greater than 1.2 C-^ 
and the temperature has been increased n times without finding 
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a flow point at which the right-hand inequality in equation 
(381) of Volume I (ref. 1) is satisfied. 

The DUMP subroutine is called to produce a dun^ of common 
data and terminate the case. 

BETA MATRIX OF INSUFFICIENT RANK 

Produced in subroutine N0NEQ. 

Followed by "ERROR EXIT NO. 1 FROM NCfcnEQ," 

Indicates that the rank of the . matrix is less than the 
number of dependent species (n-c) ; see Section 7.3.4 of Volume I 
(ref. 1) . The diagnostic is encountered only when a new gas 
model is being used. It normally indicates that too few linearly 
independent chemical reactions are included in the gas model, but 
can also occur when an error has been made in the stoichiometric 
coefficients "^'ij ^ reaction. 

The DUMP subroutine is called to produce a dump of common 
data and terminate the case. 

BOUNDARY LAYER ITERATION NOT CONVERGED 

Produced in subroutine DERIVS. 

Followed by a dump with the namelist name DRVDMP. 

Indicates that the self-consistent solution for the dis- 
placement thickness <S* and the derivatives of the flow vari- 
ables (Sect. 7.6 of Volume I) has not converged after three 
iterations . 

The nonccn verged 5* value from the final iteration is 
used and the solution proceeds. 

CONVERGENCE FAILUPsE IN AGSpIN AE = a DEL = d^ X = x 

Produced in subroutine AGS0LN. 

Followed by a dump with the namelist name AGDMP. 
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Indicates failure of the iterative solution of equation 
(135) in Volume I, In the message, AE denotes the input value 
of the effective area ratio Ag, the two values following DEL 
are the displacement thicknesses at the current flow point, 
and X is the estimate of the axial coordinate from the final 
iteration, expressed in centimeters. 

The DUMP subroutine is called to produce a dun^) of common 
data and terminate the case. 

CONVERGENCE FAILURE IN RESTMP 

Produced in subroutine RESTMP. 

Indicates failure of the iterative solution for the condi- 
tions in the upstream reservoir, based upon either the pressure 
and mass flow (for ISW2B = 0) or the enthalpy and mass flow 
(for ISW2B < 0) , Such failures can be caused by errors in the 
input data for FI^W and PRESAI or HSTAG. 

The DUMP routine is called to dump common data and temi- 
nate the case. 


DLOGR IS POSITIVE 

Produced in subroutine N0NEQ. 

Indicates that a positive value of d ^n />/dx has been en- 
countered in the nonequilibrium solution. 

Tlie subsequent action depends upon the value of the ratio 
(AFNX-DATEST)/DATEST, where AFNX is the effective area ratio and 
DATEST is a control parameter which is preset to 1.01. If this 
ratio is greater than or equal to 0.05, the message "ERROR EXIT 
NO, 3 FROM NONEQ" is written and the DUMP routine is called to 
terminate the case. If the ratio is loss than 0.05, the positive 
d /fn/3/dx value is taken to indicate that the nonequilibrium sol- 
ution is on the subsonic branch of the dovmstroam solution. An 
attempt is made to recover the desired supersonic branch of the 
solution by the inverse method at the previous switch point; see 
Section 7.4 of Volume I. 
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DUMP ROUTINE CALLED BY name 

Produced by subroutine DUMP. 

Indicates the name of the subroutine from which DUMP was 
called. Besides printing this message, subroutine DUMP sets a 
logical indicator ERR to the value .TRUE.. Tests on ERR in 
higher level routines then causes immediate return of control 
to the main program. After calling subroutine DUMPEX to print 
dximps (DMPl, DMP2, DMP3, DMP4) containing most of the variables 
in common, the main program terminates the case. Input data axe 
then read for the next case, if any. 

ERROR EXIT NO. 1 FROM NONEQ 

Produced by subroutine NONEQ, following the message "BETA 
MATRIX OF INSUFFICIENT RANK" (see above) . 

ERROR EXIT NO. 2 FROM NONEQ 

Produced by sxabroutine NJZJneq. 

Indicates step failure (flunking of a validity check) in 
a lower level routine (DERIVS or C0MM)., either in the perturba- 
tion solution or upon restart at the switch point after detec- 
tion of a positive d ^n ^/dx value in the downstream solution. 

The DUMP routine is called to dump common data and terminate 
the case. 


ERROR EXIT NO. 3 PROM NONEQ 

Produced by subroutine N^NEQ. 

Indicates that a positive d in//dx value has been encoun- 
tered in the non equilibrium solution beyond the switch point 
from the inverse method to direct integration. (See message 
"DLOGR IS POSITIVE," above.) The present error exit message 
is printed if (AFNX - DATEST)/DATEST is greater than or equal to 
0.05. 

The DUMP routine is called to dump coimmon data and terminate 
the case. 
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ERROR EXIT NO. 4 FROM NONEQ 

Produced by subroutine N0NEQ, 

This message is printed if the d p/6x> 0 condition 
is encountered when the solution has already been restarted 
at the switch point four times, in attempts to find the super- 
sonic downstream solution. 

The DUMP routine is called to dump common data and termi- 
nate the case. 


ERROR EXIT NO, 5 FROM NONEQ 
Produced by subroutine N0NEQ. 

This message is printed following the "BACKSTEPPING OF 
PERTURBATION SOLUTION" message, discussed above. 

The DUMP routine is called to dump common data and termi- 
nate the case. 


ERROR EXIT NO. 6 FROM NONEQ 

Produced by subroutine N0NEQ. 

Indicates step failure (flunking of a validity check) in 
subroutine DERIVS or C0MM at the beginning of the first step 
of the numerical integration. 

The DUMP routine is called to dump common data and termi- 
nate the case. 


ERROR EXIT NO. 7 FROM NONEQ 

Produced in subroutine N0NEQ. 

Indicates that a temperature greater than the reservoir 
temperature has been computed at a point in the nonequilibrium 
numerical integration. 

Tne DUMP routine is called to dump coimnon data and termi- 
nate the case. 
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ERROR EXIT NO. 8 FROM NONEQ 


Produced in subroutine Nj^NEQ. 

Indicates step failure in DERIVS or C0MM on restart of step 
after a step size reduction. 

The DUMP routine is called to dump common data and termi- 
nate the case. 


ERROR EXIT NO. 9 FROM NONEQ 

Produced in subroutine N0NEQ. 

Indicates that the square of the concentration of an inde- 
pendent species is zero in the element conservation calculation 
for the nonequilibrium solution. 

The DUMP routine is called to dump common data and termi- 
nate the case. 

ERROR EXIT NO. 10 FROM NONEQ 

Produced in subroutine N0NEQ. 

Indicates that the step size in the nonequilibrium integra- 
tion has become vanishingly small. This message is printed if 
Ax falls below 10”^^ cm, or if the step size has been reduced 
more than 30 successive times without completion of a valid step.* 


*Note - if this diagnostic occurs upstream of the throat (x < 0) , 
and if the terminal dump shows that d in Ae/dx = DL0GA is posi- 
tive, Ag = AFNX==1, and dT/DX = DT is near zero or positive, a 
likely cause of the failure is insufficiently rapid convergence 
of the nozzle profile curvefit upstream of the throat. In such 
cases, a successful solution can be obtained using a modified 
profile curvefit with a somewhat larger convergence angle in the 
region just upstream of the throat. The N0ZFIT program (Appen- 
dix D) can be used to prepare profile curvefits. 
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The DUMP routine is called to dump common data and termi- 
nate the case. 

ERROR IN INPUT DATA FOR NOZZLE GEOMETRY 
X = X ARATI^ = a DERIVA = d 

Produced in subroutine GE0MAR. 

Indicates that, at the axial coordinate X, the subroutine 
calculated either a geometric area ratio A„ = ARATI0 less than 1 
or a derivative d Ag/dx = DERIVA for which X • DERIVA < 0. 

The DUMP routine is called to dump coifimon data and termi- 
nate the case. 


EXCEEDED 50 ITERATIONS IN SH0CK 

\ 

Produced in subroutine SH0CK. 

Indicates that the .terative solution of the cubic equation 
for a classical oblique shock has not converged after 50 itera- 
tions. 

An error indicator in the argument list is set to inform 
the calling routine (WEDGE) of the failure. Wedge calculations 
are omitted for the angle for which the failure occurred and for 
all larger angles. The failure occurs because the assumed angle 
of attack is too high to allow an attached shock. 

FINDX CALLED WITH AN AREA RATIO LESS THAN UNITY, A = a 

Produced in subroutine FINDX. 

Indicates that FINDX was called to determine the axial co- 
ordinate corresponding to a geometric area ratio A less than 
unity. Since the geometric area ratio cannot be less than unity, 
there is no solution. Tlie error indicated by this diagnostic 
originates in the calling routine or in higher-level routines. 

The DUMP routine is called to dump common data and termi- 
nate the case. 



+++++ FIX REQUIRED IN AGSj^LN 


Produced in subroutine AGSjiLN. 

Subroutine AGS0LN solves for the geometric area ratio Ag 
and axial coordinate x corresponding to a specified value of 
the effective area ratio Ag. In channel flow problems, this 
requires an iterative solution of equation (135) in Volume I 
(ref. 1) . The above message is printed if, during this iter- 
ative solution, an x value is obtained whose sign is incon- 
sistent with that of the argument UPD0WN (which specifies 
whether an upstream or downstream solution is desired) . 

When this condition is first encountered, the programming 
assumes that it is a result of an unusually large separation be- 
tween the throat and the sonic point (caused by rapid change of 
the boundary layer displacement thicknesses in the throat region) . 
On this assumption, an attempt is made to fix the problem by re- 
setting the assumed displacement thicknesses to their values at 
the throat. If the problem recurs, the message "CONVERGENCE 
FAILURE IN AGS0LN" (above) is printed and the case is terminated 
in the usual way. 


GJ(i)** 2 UNDERFLOWED 
Produced in subroutine N0NEQ. 

Followed by message "ERROR EXIT NO, 9 FROM N0NEQ" (above) . 

INDEXING OR STORAGE FAILURE IN MATINV 
Produced in subroutine MATINV. 

Indicates that the dimension statements in arrays LPIJ and 
BTA are inconsistent; can occur only if someone tinkers with the 
programming . 

The DUMP routine is called to dump common data and termi- 
nate the case. 


IN NEWRAP, CAPX (k) = 0 
Produced in subroutine NEWRAP. 
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Indicates that the mole fraction of the kth species was 
•':ound to be zero. 

The equilibrivim solution is continued. 

IN NEWRAP, P a 0 

Produced in subroutine NEWRAP. 

Indicates that a zeio value of pressure was obtained during 
the iteration to determine the conditions at a point in the equi- 
librium solution. 

The solution is continued. 

INVALID INPUT DATA... NPRFLS = n 

Produced in subroutine READ. 

Indicates that NPRFLS was specified in the input as a value 
other than 1 or 2. 

The job is terminated. 

ITERATION TO FIND STAGNATION CONDITIONS DID NOT CONVERGE 

Produced in subroutine M0DEL, 

Followed by a dump with the namelist name M0DDMP. 

Indicates that the iterative solution for the conditions 
behind the normal shock, or the solution for the stagnation con- 
ditions, did not converge within the maximum number of allowed 
iterations. 

Conditions at the model point where the convergence failure 
occurred are not calculated. The solution continues. 

ITERATION TO OBTAIN FREE STREAM SOLUTION AT MODEL POINT DID NOT 
CONVERGE X = X X140DEL = y CM 

Produced in subroutine FR0ZEQ. 

Follov;cd by dump with the namelist name FRDMP2. 
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Indicates that the iteration to determine the free-stream 
conditions at a model point in the frozen or equilibrium solu- 
tion has not satisfied the convergence test after the maximum 
allov/able nixmber of steps. 

The conditions from the final iteration are used in the 
model condition calculations and printed out, 

MATINV, MATRIX SINGULAR 

Produced in subroutine MATIl^^^. 

Indicates that the square submatrix of defining the 

elemental composition of the independent species is singular. 

This diagnostic occurs only when a user-specified gas model is 
being used and indicates an error in the inputs defining the in- 
dependent species. It can occur if the ISC = c species listed 
first are not independent, where c denotes the number of chemi- 
cal elements present in the gas model. 

The DUMP routine is called to dump common data and termi- 
nate the case. 

MATRIX OF COEFFICIENTS IS SINGULAR 

Produced in subroutine DSMS0L. 

This subroutine is used throughout the program to solve sys- 
tems of linear equations. The diagnostic indicates that the mat- 
rix of coefficients for such a system is singular. The preceding 
output, and the data in the terminal dump, should indicate the 
kind of calculation in which the failure occurred, i.e., reser- 
voir conditions, frozen or equilibrium solution, nonequi librium 
perturbation calculation or numerical integration. 

The DUMP routine is called to dump common data and termi- 
nate the case. 

MODEL PARAMETER ROUTINE CALLED FOR A 
MACH NUMBER LESS THAN 1.5 

Produced in subroutine M0DEL. 

Indicates that M0DEL was called at a point where the frcc- 
streara Mach number was less than 1.5. The method used for the 
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normal shock solution in M0DEL does not work reliably for 
M < 1.5. 


The model condition calculations are skipped. The sol- 
ution continues. 

MORE TH;\N 50 ITERATIONS IN FINDX, A = a 
UPDj?'WN = u lENTRY = i MBL = m 

Produced in subroutine FINDX. 

Followed by a dump with the namelist name DMP. 

Indicates a convergence failure in subroutine FINDX, which 
solves for the axial coordinate x at which the geometric area 
ratio has the value A. According as the argument UPD0WN has 
the value - 1. or +1., the upstream or downstream solution is 
sought, lENTRY indicates whether the subroutine was entered at 
its beginning or through the entry point FINDXC, which is called 
to determine the x at \;hich the MBLth profile has a half-width of A, 

The DUMP routine is called to dvimp common data and termi- 
nate the case. 

NEGATIVE CONCENTRATION ENCOUNTERED IN COMM 

Produced in subroutine C0MM. 

Indicates thr C0MM found a species concentration to be 
negative during the nonequilibrium solution. 

A step-failure indicator is set, and control is returned 
to the calling routine (DERIVE) . Upon return to subroutine N0!lEQ, 
the integration stop size is reduced. The solution continixes. 

NEGATIVE OR ZERO VALUES OF :jyi0DPl NOT ALLOWED. DATA IGNORED. 


Produced in subroutine READ. 

Indicates that the input value of XM0DP1 (the location 
of the first model point in a sequence) is negative, 

XM0DP1 is reset to 10^*^ and the calculations continue. 
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NEGATIVE RHO IN GEOM 

Produced by subroutine GEOM. 

Indicates that a negative density value has been encountered 
during the iteration to deterndne the density corresponding to 
the geometric area ratio at the current ^low point during the 
nonequilibrium solution by the inverse method. 

The DUMP routine is called to dump common data and termi- 
nate the case, 

NO THERMAL PROPERTY DATA DEFINED FOR SPECIES NUMBER i IN . HE 
CURRENT GAS MODEL (NO. m IN THE M: TTER LIST OF SPECIES) 

Produced in subroutine READ, 

Indicates that ETAJ(i) = 0. and IGJ(i) = 0 for the ith spec- 
ies in the gas model, because of an input error in user specifi- 
cation of a nonstandard species. 

The jcb is terminated. 

SU2 LESS THAN 1 IN COMM 

Produced in subroutine C0MM. 

Indicates that C0MM computed a total enthalpy larger than 
the reservoir enthalpy. 

A step failure indicator is set and control is returned to 
the calling routine (DERIVS) . After the return tc N0NEQ, the 
step size is reduced and the calculation is continued, 

TEMPERATURE G1H3ATER TIIAI^ RESERVOIR VALUE 

Produced in subroutine N0NEQ, 

Followed by the message "ERROR EXIT NO. 7 FROM NOl'JEQ" (see 
above) . 

TOO MANY ITERATIONS IN VJES0LN ZETA = z CAPGAM = c 

Produced in subroutine WE;s0LK. 
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Followed by a dump with the namelist name VIEDMP. 

Indicates that the Newton-Raphson solution of equation 
(482b) in Volume I (ref. 1) for X as a function of C has 
not converged after 20 iterations. 

The final, vinconverged value of X is accepted and used. 

The solution continues. 

TCK) MANY NEWTON-RAPHSCa? ITERATIONS 

Produced in subroutine EQCALC or in subroutine NEWRAP. 

Irdicates that the Newton-Raph*: m solution for the equi- 
libriuni mole fractions (Volume I, Sections 6.1 and 6.2) failed to 
converge. If the failure occurs in EQCALC during the reservoir 
calculations (EQCALC called by the entry INTA of subroutine INGAS) , 
it ir probably caused by an error in the input specifications of 
the reservoir condition (input variables ISW2B, PRESAI, CTAPI, 

FL0W, HSTAG) . 

The DUMP routine is called to dun^ common data and termi- 
nate the case. 

TRANSPORT PROPERTIES OF DESIRED MIXTURE t.’^NNOT BE CALCULATED 
FROM AVAILABLE DATA. REVISE CROSS SECTION INPUT DATA. 

Produced in subroutine XSECT. 

Followed by a dump with the namelist name XSDMP. 

Indicates that cross section data have not been specified 
for the like- like interactions of a neutral atom or molecule. 

The DUMP routine is not called, but the case is terminated. 

X DECREASED IN FROZEQ 

Produced in subroutine FR0ZEQ. 

Indicates that the axial coordinate x decreased during a 
stop of the frozen or equilibrium solution. Since x is calcu- 
lated from the effective area ratio in those types of solution, 
a decrease in x could result from improper specification of the 
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< ‘ nozzle profile (for example, a discontinuity in the profile at 

^ - the point where two profile sections aie joined). If the sol- 

ution includes the boundary layer, rapid growth of the displace- 
ment thickness can cause the effective area ratio to decrease. 

A decrease in the effective area ratio downstream of the throat 
would lead to a decrease in x. This error condition is usually 
encountered only in the high Mach number region far downstream 
of the throat. The most common cause is instability of the coup- 
led inviscid flow and boundary layer. 

The subsequent action depends upon the circumstances. If 
the boundary layer is being neglected, or if the effective area 
ratio is less than k of the geometric area ratio, a dump with 
the namelist ncime FRDMP is written and the current equilibrium 
or frozen flow solution is terminated. The solution of the cur- 
rent case continues, however. If the boundary layer is included 
and the effective area ratio is greater than % of the geometric 
area ratio, the program assumes that the failure is a result of 
instability, and attempts to generate a valid solution by cutting 
the stability parameter w in equation (218) of Volume I in half, 
cind restarting the solution in the upstream reservoir. However, 
if the same error occurs after three successive restarts, the 
I dump FRDMP is written and the current (equil5,brium or frozen) 

solution is terminated. 

ZERO OR NEGATIVE STEP IN BLAYER, X = x XP = y 

Produced in subroutine BLAYER. 

Followed by a dump with the namelist name DMP. 

Indicates that the axial coordinate X at the current call 
to subroutine BLAYER is less than or equal to the value XP at 
the last previous call. 

The RETURN is executed and the solution proceeds. 
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APPENDIX C 

ILLUSTRATIVE TEST PROBLEMS 

This appendix presents the inputs and selected portions 
of the output for five test problems, chosen to illustrate var- 
ious features of the NATA code. Two additional test problems 
(nos. 1 and lA) have already been discussed in Sections 2 and 

3. 

The images of the input cards for the test problems are 
shown in figure 65. There are four groups of data cards, each 
group comprising the data for a NATA run. The third run includes 
two cases (4A and 4B) . Tlie se runs were executed on an IBM 360/75. 

Test problem no. 2 illustrates NATA flow calculations for a 
rectangular channel. In a channel case, the problem summary 
includes geometric data for two profiles (figure 66). Figure 
67 shows the output of reservoir conditions for this case. Fig- 
ure 68 illustrates the output of the flow solution. In channel 
flow problems including the boundary layer, NATA prints two com- 
plete sets of boundary layer data at each flow point. The 
first set . ppears in the fourth and fifth lines of output for 
the flow point, the second set in the sixth and seventh lines. 

Each set begins with the lateral dimension ("WIDTH" or "HEIGHT") 
of one of the channel walls. The boundary layer data in each 
set pertain to the wall whose lateral dimension is included in 
the set. For example, at the last flow point sho\%?n in figure 
68, the lateral dimensions of the channel are 2 by 18 inches. 

The first set of boundary layer data (including THETA = 0.100 
and STANTN = 1.697D-3) refers to the boundary layer on the walls 
which are 18 inches wide. The second set (including THETA = 0.340 
and STANTN = 4.777D-4) refers to the layer on the walls which are 
2 inches wide. The flow points at X = 16.000, 17.000, and 18.000 
in figure 68 are special points requested by the TSDIAM inputs 
in figure 65. 

Tes . problem no. 3 illustrates a flov/ solvition based on 
the larger of the two standard planetary atmosphere models. 

In this case, the reservoir conditions were specified by direct 
input of the reservoir pressure and temperature. Figures 69 
through 73 sh'>'w’ the problem summary for this case. Figure 74 
shov;s the reservoir conditions, and figures 75 to 78 the first 
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5 FIGURE 65 - INPUT DATA FOR TEST PROBLEMS NO. 2, 3, 4A, 4B, AND 5 

TEST PROBLEM NO. 2 > FROZEN AIR FLOW IN A CHANNEL 
♦INPUT 

1SW2A=0. tSW3A=0* CXMAXI=57* PRESAI=.762. FL0W=. 1 • ICHANst. 

TSDIAHslS. 16. 17. 18, 19. 20 
♦END 


TEST PROBLEM NO. 3 - PLAf<CTARY ATMOSPHERE MODEL 
♦INPUT 

ISW1A=0. ISW3A=0. CXMAXI=50. ISW2B=I. PRESAI=1. CTAPI=10000. NOZZLE=2. 
IGAS=5 
♦END 


TEST PROBLEM NO. 4A - WEGENER EXPERIMENT C - NONSTANDARD GAS AND GEOMETRY 
♦INPUT 

ISW2B=1. NOTRAN=T, CXM AXI=:5.906. TPRNTI=.001, READG=T, PRESAI=2. CTAPI=4D0. 
N02ZLE=0. DIAM=.5472, JOIM=0. NSECTS=2.2. I SHAPE = 1 .2 .2 • t , 
ARAMI=7.60fe602E“2, -.5773502. 0. 4,695442. 0. 4.0005. 4,695442. 0. 

4,0005. .6941124. 2.0368255-2. 0. A TPI=-2, 000249, 0, .08154619. 

XZEROI=-10. IGAS=0. NCi>=2, JCS=5.30, QPJ=. 995025. .004975, 1SCI=2. 

ISSI=3, ISRI=1. ICI=0, IE=5.6. 15=5.30.29. 1R=76, ISAT0M=29, ISMOL=30, 
ISW4A=1 . 

♦END 
♦E I NPUT 

SP29=0 ,2.5,6. C, 1, 2,0, 4.003. -3.75E-4, 2.45E-6. 2*0, 5.945. 3586. 4*0. 

1,2340. SP30=0.2,5.6.0.2.4.0 .3. 55 3. .Oil 625.-4.55E-6 .2 *0 . 10 . 028 . 4473. 

440. 1, 2340. RP76=3,E14. 340. 2. 2 , 5 . 29 ,0.5 . 30 • 0 . 1 .2 » 0 . 1 ,1 , 0 .0 . 10 40 . 

♦END 

TEST PROBLEM 4B - WEGENER EXPERIMENT F STACKING OF CASES 

♦INPUT 

PRESAI = 2.16. CTAPI=402, QPJ=.97S61. .02439, 1SW4A=0. READG=F 
♦ END 


TEST PROBLEM NO. 5 - ELECTRONIC NONEQUl L IBR lUM MODEL FOR ARGON 
♦INPUT 

ISW1A=0. ISW3A=0. ISW3B=0, CXMAXI=40, ISW6B=-1, ISW2B=1. PRFSAI=1, 
CTAPI=10000. NOZZLE=l, IGAS=3. XMOOPl=20, NMODPT=3, TPRNTI=,001, 
ISW5B=-1COOOOO 
•MD 


OP POOR QUALITY 
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PieSURB 68 - FLOW-SOLUTIOM p^*?OT FOR TEST PROBLEM NO 
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four pages of the nonequilibrium flow solution. In this case, 
the switch from the inverse method to direct integration was 
unsuccessful on the first try, as indicated by the diagnostic 
"DLOGR IS POSITIVE" in figure 77. After an underflow message 
produced by the IBM operating system, MATA printed the conditions 
at the flow point where din p/dx was positive (top of figure 78) 
and restarted the solution at the previous switch point. After 
the restart, the swi-cdiover from the inverse method to direct 
integration occurred sufficiently far downstream of the throat 
to give the desired supersonic branch of the downstream solution, 
as shown in figure 78. 

Test problems no. 4A and 4B simulate two of Wegener's ex- 
periments on MO 2 recombination. The MATA solutions of these 
cases are shown in figures 29 and 30 of Volume I (ref. 1) . 

These problems illustrate the use of nonstandard nozzle geom- 
etry, gas species, and reactions. In the input (figure 65), the 
geometric data describe Wegener's wind tunnel (shown in figure 
28 of Volxime I ). Species number 29 is MO 2 and nunber 30 is 
N 2 O 4 . The recombination reaction 2 MO 2 + M 2 M 2 O 4 + M 2 is de- 
fined as reaction number 76. Figures 79 and 80 show the problem 
summciry for case 4A. Mote that the standard properties are 
used for M 2 . The theannal properties of MO 2 and M 2 O 4 aure defined 
by means of thermo fits; no physical model data are given for 
these species because MATA is not programmed to treat nonlinear 
triatomic molecules such as MO 2 or polyatomic species such as 
M 2 O 4 . All transport property calculations were suppressed in 
this run by input cf M0TRAM = T. Figure 81 shows the calcu- 
lated reservoir conditxons, and figure 82 the first page of the 
equilibrium solution. 

Test problem no. 5 illustrates the use of MATA with a sta. 
dard electronic nonequilibrium model (argon, IGAS = 3), together 
with some of the output controls. Figures 83 to 85 show the 
problem summary. Note that the electronically excited species 
AR*M and AR*R have the same alpha matrix (elemental composition) 
as the ground state species AR, but have different enthalpies of 
formation and appear in different reactions. The electron ther- 
mal nonequilibrium parameters are tabulated in figure 85. Fig- 
ure 86 shows the reservoir conditions. The boundary layer was 
neglected in the solution. Figures 87 and S 8 show two pages of 
the output from the nonequilibrium solution. As shewn in figure 
87, the code begins the solution by taking three steps using the 
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perturbation echnique. The numerical integration then begins. 
Stability requirements force the use '>f a very small step size. 

The input ISW6B = -1 in test problem no. 5 (figure 65) 
gives the output of reaction rate data (PI, CHI, PICHI, DIG) in 
each step (see Section 3.5). Tbe input ISW5B = -1000000 gives 
a one- line message (X, DELTAX, T, etc.) for every completed 
integration step. The data included in this message are as fol- 
lows: 

Axial coordinate (cm) 

Integration step size (cm) 

Heavy-particle tenperature (°K) 

Electron temperature (°K) 

Nondimensional sta*~ 'ation enthalpy 

Energy transfer to thv, electron gas (cal/cm^sec) 

Nronber of step— size reductions required to 
atdiieve a successful integration step 

Thu messages "CAUSES OF STEP FAILURE" in figures 87 and 88 
are also triggered by the negative ISW5B value, and indicate the 
location and nature of the flurhed validity check responsible for 
each step size reduv-tion. The nuiaarical code used in this mes- 
sage is documented in Section 4.55 of Volume III of this report. 


X 

DELTAX 

T 

TEP 

CHA 

QDPE 

IC0UNT 
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APPENDIX D 


USER'S MANUAL FOR THE NOZFIT C(X)E 


D. 1 Introduction 

The geometries of nozzles and channels are specified, in 
NATA, by analytical curvefits to the profiles of the surfaces 
confining the gan flow. In these curvefits, each profile is 
represented by means of a sequence of straight line segments and 
circular arcs. The details of this system of geoioetric specifi- 
cations have been documented in Section 4 of Volume I (ref. 1) . 

In order for the curvefits to be usable in NATA, the various 
sections must join together continuously, and with continuous 
slopes, to high accuracy. The NOZFIT program has been developed 
to facilitate the preparation of such curvefits. The main inputs 
to NOZFIT are geometric data, most of which can be read directly 
from nozzle design drawings. The outputs include 

(1) Printout of the parameters in the profile curvefit; 

(2) Pvinched cards, containing a Fort -an DATA statement, 
which can be incorporated directly into NATA, to add 
the curvefit to the set of precoded standard profiles; 

(3) Printout of the throat region of the profile in tabu- 
lar form; and 

(4) A coit 5 >uter-gene rated plot of the profile in the throat 
region. 

Section D.2 c this appendix defines the inputs to N02FIT, 
and Section D.3 discusses the various types of output. Section 
D,4 presents a sample rvin of the code with its printed and plotted 
output and a listing of the punched output. 

D.2 Inputs to NOZFIT 

Tlie profile curvefits produced by NOZFIT consist of sec- 
tions joined end to end with continuity of slope, each section 
being a straight line, a circular arc concave upward, or a cir- 
cular arc concave do\mward. A,.ray dimensions in NATA limit the 
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total number of sections to 12. In the analysis to detemine 
the profile curvefit parameters, the sections are separated into 
two groups, those lying upstrecim of the throat, and those lying 
dovmstream. The analysis for each group starts at the throat 
and determines the sections sequentially, proceeding away from 
the throat. In each group, the section having a boundary at the 
throat (x = 0) is always assumed tc be a circular arc concave 
upward with zero slope at the throat. Thus, the profile ordi- 
nate is a minimum at tlie throat. For the profile to be usable 
in NATA, it is necessary that the ordinate be a monotonically 
decreasing function of x upstre m of the throat and a monotoni- 
cally increasing function of x downstream.* 

The input to a NOZFIT case begins with a card containing 
alphanumeric information describing and identifying the nozzle 
or channel whose profile is to be fitted. This card is read 
with an A format. The information is printed at the head of the 
output and reproduced as a conraent card preceding the data cards 
produced. The first 48 characters are also used as a title for 
the nozzle profile plot. Finally, the first 4 characters on 
this header card are incorporated into the punched DATA state- 
ment and are used, by NATA, as a facility name (DC.., etc.). 

The remaining inputs are all read in under the namelist 
name INPUT, using the usual namelist format (see Section 2.1). 
They are defined below: 


Variable Preset 

Name Dimensions Values Definition 


NSECTS(I) 2 


2*0 NSECTS(l) = 
in nozzle 
throat. 

NSECTS(2) = 
in nozzle 
of throat 


number of sections 
profile upstream of 

number of sections 
profile downstream 


*More precisely, che area ratio A must have dAg/dx < 0 for 
X < 0, dJ\g/dx > 0 for x > 0. In a^channol, it is not necessary 
for both of the profiles t ecrease and increase as described 
above, so long as the area ratio has the required behavior.. 
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Variable Preset 

Name Dimensions Values Definition 

ISHAPE(J) 12 12*0 Index specifying shape of the Jth 

section (counting frcxa the up- 
stream inlet) 

1 Straight section 

2 Circular arc convex toward 
axis 

3 Circular arc ccncave toward 
cixis 

PAR(I,J) 2,12 24*0 Parameter values for the Jth sec- 

tion: 

For ISHAPE(J) = 1, PAR(1,J) = 
angle of inclination to nozzle 
axis in degrees (positi. e value) 
For ISHAPE(J) = 2 or 3, PAR(1,J)= 
circle radius in inches 
(See ICOND for PAR(2,J)) 

- Throat diameter in inches. 

12*0 Index specifying condition de- 

fining the Jth section. 

IC0ND(J) = 1 Throat condition* 
ICjeiND(J) = 2 Straight section (J) 
is tangent to adjacent circular 
section nearer the throat. 
IC;6 nd(J) = 3 Circular section (J) 
is tangent to adjacent straight 
section (nearer the throat) at 
an axial distance of PAR(2,J) 
inches from the throat** 

IC/)ND(J) = 4 Circular section (J) 
is used to break a sharp angle 
between two straight sections 
which intersect at an axial dis- 
tance of PMl(2,J) inches from 
the throat** 


*Noto - there are always two throat sections, one upstream and 
i < ne do\imstream of the throat. 

**Note - PAR(2,J) is negative if it represents a point upstream 
of the throat. 


DTK 1 

J.C0ND(J) 12 
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Variable 

Name 

Dimensions 

Preset 

Values 

Definition 

XSTART 

1 

* 

Upstream limit on V. for calcula- 
tion of nozzle profile (negative 
value, inches). The profile is 
calculated for a nozzle section 
6 inches long. 

XZER0I 

1 


Inlet position in inches above 
the throat (negative) ; for use 
in NATA boundary layer calcula- 
tions . 

N0ZZLE 

1 

- 

Nozzle index for use in NATA 
(integer value between 1 and 20) . 

CARDS 

1 

.TRUE. 

Set to .FALSE, to suppress card 
outgput. 

PI^S 

1 

.TRUE. 

Set to .FALSE, to suppress plot 
output. 

ENDJ03 

1 

.TRUE. 

Set to .FALSE, if there is another 
case in the job follov/ing the cur- 
rent case. Set to .TRUE, in last 
case. 


D.3 Outputs of NOZFIT 


The outputs of N0ZFIT are illustrated by the results for 
the test problem in the next section. Tlie printed output con- 
sists of the folic ’ing: 

(1) A listing of the values of all the input variables 
in Namelist format. 

(2) The image of the "header” card containing alphanum- 
eric identifying infcrmation. 

(3) A table headed "nozzle profile parameters," This 
tctblc contains the following five columns: 
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ATP 


PARAM(1,J) 
PARAM(2,J) 
PARAMO, J) 


Nozzle section index, starting 
fr<xn the. upstream end. 

Position coordinate of the down- 
stream boxindary of the profile 
section (cm, positive downstreeua) , 

See below. 

See below. 

See below. 


The parameter values are coefficients in the emalytical 
e 3 q>ressions for the profile sections. Let Pq^= PARAM(1, J) 
P 2 sPARAM( 2 ,J), and P 3 = PARAM(3,J) . Then: 

For ISHAPE(J) = 1, 

y = Pi + P 2 X (D-1) 


For ISHAPE(J) = 2, 

Y = P]L " - (X - P2)2 (D-2) 

For ISHAPE(J) = 3 

Y = Pj^ +/ P 32 - (X - P 2 )^ (D-3) 


v/here Y is the profile ordinate and X the axial co- 
ordinate, The units of P^^, P 2 and P 3 are such as to 
yield Y in centimeters when X is ejq>ressed in centi- 
meters. For ISHAPE(J) = 2 or 3 , the profile section 
is a circular arc, P 3 is the circle radius, and the 
circle center is at X = P 2 f Y = Pj^, 

(4) Printed images of the DATA cards produced. 

( 5 ) A two-column table headed "nozzle profile," giving the 
X- and y-coordinates (in inches) of points on a section 
of the profile 6 inches long, beginning at XSTART. Tlieso 
arc the semio data that arc represented in the plot output. 


- 2 - 71 - 



1 


The punched output consists of 12 punched cards for each 
case. Their fomat is illustrated by the printed card images 
in the output fo" the sample problem. They consist of a comment 
card followed by an eleven-card DATA statement defining an ar- 
ray ZPn, where n is the profile index (equ-'l to the NOZFIT input 
N^ZZItB) . The newly fitted profile can be incorporated into 
NATA by inserting tliis DATA statement into the block data rou- 
tine BLKDl and recompiling the routine. Note that n must be 
less than or equal to 20, and must be different from the indices 
of all the other profiles defined in BLKDl. 

All of the entries in ZPn are floating-point numbers. Some 
of these values represent integers. To ensure rounding-down to 
the correct integer values, N0ZFIT increases such values by 0.1 
in the DATA statement. The entries in ZPn are all defined in 
Section 4.7. 

The plot output consists of one plot per case, showing a 
section of the ■ mputed profile 6 inches long in the axial dir- 
ection. The CO ordinate scales along the x and y axes are the 
same so that the shape is not distorted. The plot is a little 
larger than full scale. 

D.4 NOZFIT Test Problem 

The following pages (figures 88-93) present -e printed 
oulpi’t from a N0ZFIT run on the IBM 360/75 at Avco Systems Div- 
ision. The values of the input variables for this -^est problem 
are listed on the first page of the output. Figure 55 shows the 
plot produced. 
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FIGURE 90 - OUTPUT OF N0ZFIT TEST PROBLEM (Second Page) 
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FIGURE 91 - OUTPUT OF N0ZFIT TEST PROBIiEM (Third Pige) 
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FIGURE 92 - OUTPUT OF N0''""’T TEST PROBLEM (Fourth Page) 
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FIGURE U3 - OUTPUT OF N0ZFXT TEST PROI3LEM (Fifth Pago) 
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